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(16514 SHEAR LAG IN BOX GIRDERS 
WORDS: Bending; Box ox girders; Concrete (reinforced); Load wa 
distribution; Loads (forces); Reinforced concretes; Shear lag; Slabs a a 

ABSTRACT: An initial approximate method is presented to evaluate the shear lag _ 
_ phenomenon in reinforced concrete single-box girder bridges. Assuming an ideal cross- 
section and considering a symmetrical deformation of the top and bottom slabs, but q 
_ preserving continuity conditions at the web-slab intersection points, a fourth order non- 
homogeneous differential equation with constant coefficients was obtained. Making 
_ further simplifications, quick results are obtained for uniformly distributed loads. In 
_ spite of all the simplifications these results were in good agreement with the German 
_ specifications for shear lag provided that midspan reduction factors are desired. the — 
| aan are the most important anyway, because they are valid practically almost for : 
REFERENCE: Kusmenovic, Bogdan O. (Prof., Dept. of Civ. Univ. of 
_ Lawrence, Kans.), and Graham, H. James, “Shear Lag in Box Girders,” Journal of the 

4 _ Structural Division, ASCE, Vol. 107, No. ST9, b aang Paper 16514, September, 1981, 
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' KEY WORDS: Analytical techniques; Cyclic loads; Damage estimation; oe 
Earthquake resistant structures; Inelastic action; Reinforced 
_concretes; Seismic design; Structural deformation analysis; 
| ABSTRACT: Seismic damage in reinforced concrete frames is predicted by onctetinal 
methods. Analytical models are developed which can reproduce the behavior of 
| reinforced concrete members under cyclic loads. Accuracy of the analytical models is 
1 checked by considering a set of quasi-static cyclic load tests, and the models are found 
2 to be sufficiently accurate. Several damage indicators such as ductility, flexural damage 
_ ratio, dissipated energy, and cumulative plastic rotation are investigated for each 
_ experiment. Based on the test results and physical grounds, flexural damage ratio and Ay o 
dissipated energy are chosen as damage state parameters. The results of the present — 
_ study are later used to develop a probabilistic model of member —— in pat mo | 
REFERENCE: —Hooshang( Structural Mechanics Associates Beach, 
Calif.), Biggs, John M., and Irvine, H. Max, * “Seismic Damage in Reinforced Concrete _ 
Journal of the Structural Division, ASCE, Vol. 107, No. ST9, Proc. _— _ 


16496, September, 1981, pp. 1713-1729 


it 16496 | SEISMIC DAMAGE IN REINFORCED CONCRETE FRAMES _ ia 


16520 COLLAPSE OF AIR-SUPPORTED SPHERICALCAPS 
WORDS: Air ‘structures; Collapse; Laboratory tests; 
- Ponding; Ponding tests; Structural behavior; Structural design; Structural 
ABSTRACT: A theory is presented to predict the critical central weight which, in the © 
| presence a ponding medium, will lead to collapse of an inextensible air supported 
a spherical membrane. Results are developed for a range of radii, pressures and ponding 3. 
| medium densities. Also explored is the manner in which these results will be effected 
- by meridional action only, such as would be caused by a system of cables converyag 
7 at the apex. A description is given of a series of tests carried out on two laboratory — 
models to confirm the theoretical results. The experimental results indicate that the 
theory may lead to under design and that this instability is particularly sensitive to 
surface irregularities near the apex, 
~ REFERENCE: Malcolm, David J. (sr. Engr., DAF Indal Ltd., 3570 Hawkstone Road, 
Mississauga, Ontario, LSC 2V8 Canada), and Glockner, Peter G., “Collapse by | 
Ponding of Air-Supported Spherical Caps,” Journal of the Structural Division, ASCE, 
Vol. 107, No. ST9, Proc. Paper 16520, September, 1981, pp. 1731-1742 eqadl 
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16521 PROCEDURE FOR PULSE- SIMULATED RESPONSE 
KEY WORDS: Algorithms; Dynamic response; Earthquake simulation 
Impulses; Pulse frequency; Pulses; Structural dynamics; 
Stability; Vibration response 


- ABSTRACT: A computationally efficient technique is presented for optimizing the - 
| selection of pulse train characteristics, to be used for simulating the response of general 7 


| types of structural systems to arbitrary dynamic environments. The optimization 
_ procedure uses an adaptive random search algorithm that incorporates a periodic 
exploratory search for the optimal step-size variance, which signifianctly improves the 
cenvergence characteristics. Pulse trains (to induce specified or criterion structural 
motion) are not unique and are largely dependent upon constraints. The method * 
REFERENCE: Masri, Semi F. (Prof., Cw. Engrg. Dept., of Southern ag! Los 
ngeles, Calif. 90007), and Safford, Frederick B., “Optimization Procedure for Pulse- _ 
imulated Response,” Journal of the Structural Division, ASCE, Vol. 107, No. ST9, 
roc. Paper 16521, September, 1981, pp. 1743-1761 
fort 


16486 APPROPRIATE DESIGN BY GEOMETRIC SERIES 


KEY WORDS: Approximation method; Frome design; ‘Frames; — 
distribution; Optimization; _ Rectangles; Series (mathematics); Structural a 
ABSTRACT: An approximate optimization method is presented that is applicable to- 
rectangular, two-dimensional multistory frames. Using special moment distribution - 
techniques, the optimization of linearized frame “weight” is achieved for an infinite 4 
number of iterations. A closed-form exact formula makes it possible to avoid these — 
iterations. Thus, when this program is used, a minicomputer with 4K memory can 
handle a four-story, four-bay frame subject to coepnenetnest loading. The method is 

REFERENCE: Kuzmanovic, Bogdan O. (Prof. of Civ. _ Engrg., Univ. of Kansas, 

_ Lawrence, Kans. 66045), and Willems, Nicholas, “Appropriate Design by Geometric — 
Series,” Journal of the Structural Division, ASCE, Vol. 107, No. ST9, Proc. = 
16486, September, 1981, PP. 1763-1775 


ABSTRACT: After reviewing general methods available fo: determining seismic design ; 

, forces for structures which can tolerate limited amounts of inelastic deformations, the 

a reliability of two representative procedures is evaluated. In the methods evaluated, 
inelastic design response spectra are obtained by modifying a linear elastic design 
response spectrum in terms of a specified ductility factor as suggested by Newmark 

| and Hall and the Applied Technology council. The effect of different accelerograms, as ; 

well as of different system damping and hysteretic characteristics, on the inelastic — 

- response of single degree-of-freedom systems designed using these methods is 
thoroughly investigated, considering maximum displacement ductilities, maximum and 
permanent drifts, number of yield events, and hysteretic energy a 
Application of these design methods to multistory buildings is also briefly considered 
Results obtained for ideal inelastic systems indicate that these methods do not Er ad : 


REFERENCE: Mahin, Stephen A. (Asst. Prof. of ‘Civ. Engrg., Univ. of California, 
i Berkeley, Calif. 94720), and Bertero, Vitelmo V., “An Evaluation of Inelastic Seismic , 
| Design Spectra,” Journal of the Structural Division, ASCE, Vol. 107, No. ST9, Proc. ; ; 
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46495 CELLULAR AND VOIDED SLAB BRIDGES 
KEY WORDS: Bridges (structures); Cellular structures; Characteristics; 
Onthotrople (structural); Longitudinal force; Longitudinal stability; 


ABSTRACT: Bridges of the voided ‘slab type a are ‘explored. These types of 
_ bridges have either rectangular or circular holes running in the direction of the bridge 
= ‘span. In these structures the transverse distortion of the cross section may significantly _ 
| affect the load distribution properties of a bridge to the extent that it can no longer be 
| analyzed by conventional orthotropic plate theory. A simplified method is developed 
7 | which can account for the increases in longitudinal moment and shear intensities which 
_ result from distortion of the cross section. Existing simplified methods, which are 
- ~ applicable only when transverse cell-distortion is absent, may be used in conjunction — 
certain factors even when such such effects are 
Bakht, (Sr. Research Policy Planning and Research 
- - Div., Ontario Ministry of Transportation and Communications, Downsview, Ontario, 
ae Jaeger, Leslie G., and Cheung, Mo S., “Cellular and Voided Slab Bridges,” 


eS Journal of the Structural Division, ASCE, Vol. 107, No. ST9, Proc. c. Paper 16495, 


16512 BEAM-TO-COLUMN CONNECTIONS tes 
KEY WORDS: Beams _ (structural); Connections (joints); Fabrication; 
_ Flanges; High strength nad Stiffeners; Structural — Structural steel; 
Tests; Welds 
ABSTRACT: review is of tests on beam-to- ‘column moment ‘flange 
connections. The main objective of these tests was to study the behavior and to 
“te develop the design method for flange moment connections. A test program of flange — 
Bt connections consisted of 12 full-size beam-to-column connections. This test series 
7 included specimens representing connections tor the lower stories, middle stories, and 
upper stories of a multi-story frame. Primary attention is focused on strength, 
‘ deformation capacity and overall stiffness of the connection. The load-deflection plot 
2 of each connection is compared with its control test behavior, in the adequacy of AISC © 


' specifications is studied. 
| REFERENCE: Chea, W. F. (Prof. of Struc. Engrg., School of Civ. Engrg., Purdue — \. 
1 Univ., | West Lafayette, Ind. 47907), and Patel, K. V., “Static Behavior of Beam-to- i. 
£ Column Moment Connections,” Journal of the Structural Division, ASCE, Vol. 107, 
No. ST9, Proc. Paper 16512, September, 1981, pp. 1815- 


16504 LIMIT DESIGN OF SLABS FOR CONCENTRATED LOADS | 
KEY WORDS: Design; Limit design; Loads (forces); Orthotropism; 
Plasticity; Plates (structural members); Reinforced concretes; Reinforcement 


ABSTRACT: Yield Line Theory is used to derive expressions for the flexura2 limit 
loads of reinforced concrete slabs subjected to simultaneously acting concentrated and 
distributed loads. The “concentrated” loads are assumed to have a rectangular 
footprint of any reasonable size and aspect ratio. The slab reinforcement can be 
orthotropic, different over each supported boundary, and the positive and negative 
moment reinforcements can also be different. The slab boundaries may be rigidly : 
supported or consist of beams whose bending and torsional strengths, and weight, can _ 
be incorporated into the analysis. One of the boundaries can also be unsupported. The _ 
geometrically admissible collapse mechanisms for the combined loadings are studies 
and results obtained by their use are compared with some presented in the literature 


REFERENCE: Gesund, Hans (Prof. of ‘Struct. Engrg., , Univ. y. of Kentucky Lexington, 
Ky. 40506), “Limit Design of Slabs for Concentrated Loads,” Journal of the Structural — 
_ Division, ASCE, Vol. 107, No. ST9, Proc. Paper 16504, September, 1981, pp. 1839- 
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“WORDS: ‘Buckling; Buildings; Cyclic loads; E: 

Engineering mechanics; Plastic properties; Steels; Stretching; Structural 
ABSTRACT: An experimental study is presented of the inelastic behavior of axially 
loaded steel members subjected to repeated buckling and stretching. Tests on 24 
commercially available steel struts commonly used as bracing members are reported. A 
large variety of shapes were tested including wide flanges, structural tees, double- 
angles, a double-channel, and thick and thin-walled square and round tubes. The 
boundary conditions were of two types, pinned-pinned and fixed-pinned, while the 
effective slenderness ratios were either 40, 80, or 120. The effects were investigated of 
loading patterns, end conditions, cross-sectional shapes, and slenderness ratios on the 
hysteresis response of members. An explanation is given regarding the fundamental 
_mechanisms responsible for the observed degradation in the buckling load capacity 


REFERENCE: Egor P. (Prof., Dept. of Civ. Univ. of California, 

- Berkeley, Calif. 94720), and Black, R. Gary, “Steel Struts under Severe Cyclic 

_ Loadings,” Journal of the Structural Division, ASCE, Vol. 107, No. ST9, Proc. Paper | a 

(16497, Se tember, 1981, 1857-1831 
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SHEAR LAG IN Box GIRDERS 


By Bogdan O. Kuzmanovic,' F. ASCE and and H. James Graham,” ASCE 


_ _ Shear lag phenomena, or the nonuniform distribution of a stresses nan 
wide flanges of a beam cross section, have long been recognized (3,5,6,12). 
First studies were ‘mainly | concerned | with steel beams of J of T shape without Y 
composite action (1,3,5 ,6) while later. investigations dealt with composite action 
— 2, 4,9). When box girders made their comeback in bridge construction, most 
of the studies either considered only steel box sections (7,16,18,19) or sections 

with a wide concrete slab acting with the steel part of the box 


reinf orced concrete box sections of the dimensions usually e encountered in — 
bridges were discussed very rarely (13,17). The American Association of State 
Highway and Transportation Officials, AASHTO, specifications (3) draw attention | 

s the shear-lag problem, but do not give any guidance on how to ‘proceed. 4 


The German Industrial Standards, DIN, issued in August, 1973 (8), give the 
effective width of reinforced concrete box sections for uniform load and various : 
Support conditions. These widths for uniform loads are in agreement with test - 
results (21), and they are on the conservative side. Presently, most studies oe 
on shear lag have been performed using computerized finite element or folded : 
plate analysis (16,17,18,23). A quick evaluation of the possible shear lag effect 
is of importance to a practicing engineer at the early stage of the design of | 
a box girder bridge. A computer run at this stage is . neither feasible nor economical — 
as even the cross section itself might be changed in further studies. Therefore, — 
_ there is a real need for an approximate but quick method method of reasonable accuracy ‘ 
entailing a limited amount of numerical work. 
In this paper, the writers | develop a simplified method for bridge sections 
of the type shown in Fig. \(a) using the principle of minimum potential energy 
and adapting the method ony Ss by E. Reissner (19) for the cross — 


The effective width | is a function of geometrical factors, such as slab width, 
_effective span, and the type of stress distribution. Boundary conditions at bridge 


_'Prof., Dept. of Civ. Engrg., Univ. of Kansas, Lawrence, Kans. __ 


3 Assoc., Beiswenger, Hoch, and Associates, North Miami Beach, Fla. n OF wih. 
_ Note.—Discussion open until February 1, 1982. To extend the closing date one month, 
7 . written request must be filed with the Manager of Technical and Professional Publications, 3 
ASCE. Manuscript was submitted for review for possible publication on October 8, 1980. 
This paper is part of the Journal of the Structural Division, Proceedings of the American q 
7 Society of Civil Engineers, ©ASCE, Vol. 107, No. ST9, September, 1981. ISSN 0044- 
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"supports and the type of loading, e.g., point load, load, and prestressing 
are also. affecting this width. A few codes that do consider shear lag do not 

m take into account the load parameters. Similarly, the mutual dependence ce 
the effective width and the distribution of moments along the bridge axis in 

practical bridge design is either completely disregarded, or the total effect of 
he moment included, e.g., as a fixed reduction or increase 


2 


moments, or both. Again, the reason lies in considerable amount 


numerical work in solving this problem. 


To reduce the amount of work in developing the method ideal 


| 
. 

Hl 

onUx) 

—  —_ 

1 
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will ‘be applied to an . actual bridge section, Fig. 3. _The top slab h has a width 
of 4w and a eave thickness t,. The bottom slab is 2w in width and has ‘ 
a thickness tye The web has a thickness t,, and depth h. A | | symmetrical load a 
distribution | q(x) is specified along the span L _ and it acts | symmetrically on pe 
the webs, 1 normal to the slab plane. This load causes a known distribution of : 
_ bending moments along the length of the span, M(x), neglecting any influence 
7 7 of the effective widths. Both slabs are assumed to be in a two-dimensional 
__ stress state, usually identified as a plane stress state. The model is further 
_ simplified by assuming the slab rigidity in the chordwide (y- ~axis) direction to 
be large while the -Spanwise (x-axis) slab displacement u(x, y) is assumed 0 
be pon symmetrical for all four segments of w as shown in Fig. Ney: 


‘+ (1-4) 


A 
_ ‘The va variation of u(x,y) of course for a cantilever slab is different from that — ; 
for a slab between wets, | However, this assumed symmetry cuts down numerical __ 


; 


work considerably and, as examples demonstrate, is still 
_ The continuity of displacements at the intersection between web and slab — 
(in which y = w) is preserved in Eqs. | and 2. At points where y = 0, the | 
spanwise displacement u u(x, 0) is expressed in terms of an unknown function 
U(x). This function is a measure of the magnitude of the shear lag effect. 
E For the points between y = 0 and y = w a parabolic variation of u(x,y) is 
_ assumed. By applying variational calculus to the total potential energy, two. 
_ differential equations are obtained from which the unknown function U(x) and _ 
_ its derivatives can be eliminated. Therefore, the final differential equation is § 
in | terms of derivatives of ‘the deflection: z z only. Although its solution 7 
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= -—The total potential nergy | ‘Tho of the load load system, slabs and webs 


mis divine + I, is the moment of inertia of the whole cross section, i.e., 


4 


t 


From the minimizing condition, = 0, three differential equations are 
_ which enable the elimination of the unknown function U(x), its” — 
and the satisfying of its boundary and continuity conditions. 


Using Reissner’s two parameters k and n, in which 


'2"(x) = 


This is the same of R s ‘equation b but be ‘now 
different values. 


r fixed end (slab attached to support) z #©) = 


slabs not attached te to "(x) = 

By expressing the stress in the slab at the a with the webs in the | 

form o,, = Ee,, and noticing that ¢, = cz”, where c is the distance of the 


point wi under consideration from the y- axis, ‘the stress including aoe lag e effect 


The use of the method is illustrated in the following example. vied re * 
oe Example.—The shear lag of the first and second span of a three-span continuous 
bridge 106’ + 188’ + 106’ (32.3 + 57.3 + 32.3 m) is considered. Only waiform 
load is applied. The somewhat idealized cross-section is in Fig. 
= 84 ft* 
6728; n 


q 
1 
tk t 2(x 
is 
| 
| 
k = 0.1983; kl/2 = 10.511; G/E = 0417, 
— 


SHEAR 


ment 7 
(alseady simplified by disregarding the shear effect) is replaced 
“4 by that for a propped cantilever (Fig. 4). In the beam on elastic foundation — } 
analogy method used for the design of box-girders (22), the solution is approxi- — 
mated is a mean value between the solution for free and fixed supports. He Here 


P Soll" bd Bua | 


for uniform load and not too unbalanced spans, the assu ption ofa rigid ee: 


Support is close ‘enough. The bending moment at any x xis 
‘The differential “equation is is Eq. 6 
solution is straight- -forward, only the results are given. . The complemen- 
taryfunctionis 


+C,e" 


= 7 4 


at xsl; 2=0; 2’ =0; 2 
constants of of — are 
a 
the values from Eqs. 16 are substituted into into Eq. 1 12 and Eq. 13 
added together, the following expressions for z and z” are obtained: iy, Meeps 
16EI 2k? ED 


| 
| 


oarow’ the value for z’ ’ from Eq. “18 into Eq. 9, the stress in the mid-slab 


at the left support (x = 1 
¢t ‘This ‘means that no shear lag can be assessed where this stress is zero. ‘This 


is the deficiency of this method. 
At mi 
| =ql’/16 16 we 


1) 


using the values from Eq. 5 


a" 


_ DIN 1075 (8) gives for a top fanater slab a value of 0.920, a value that is 
“ ‘slightly more conservative. The width reduction fact_r at the first interior opens 
will be obtained from consideration of the second span. — 


Ps 


My, Second | Span. —Again, in order to minimize the amount of numerical work, = q 


this span is taken as fixed at both ends (Fig. Baa 


sas + 


ECO 
i 4 
2.724 x - and M,,.__ 
> 
Mp . Proceeding in the same way as for the first span, the expression for the 7 A 
m the midspan). 
d 


The coments 


At the left support 


The width reduction factor now pe Oe 762; — 1075 (9) gives a value of 
0. 740 <0. 
At midspan (x = 


factor is is B = 


is very close. iy 


1 The proposed approximate method developed for an idealized cross- “section 
“yield initial values about expected shear lag in a box girder design. = 

minimum potential energy principle Provides a a suitable approach to 7 
_ 3. In case of uniform loading and using “several simplifications, it appears 
_ that manual calculstion is still practical, yielding width ame factors close 
to given by DIN 1075, although less conservative. 

_ 4, The same method can be applied for other types of aie by integrating © 7 

-within several domains where M(x) is a continuous function. ~The 


“now, however, are much longer due to added continuity conditions betwe 


The results of this paper were obtained as part of work done by the Bridge — 7 


_ Section of the firm Beiswenger, Hoch, and Associates, Consulting Engineers, 
North Miami Beach, Florida. For permission to ane this paper the authors — 
are indebted to Mr. F. A. Norofia. 
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Al 
0 the 1075 (8) value is — which 
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s 


seer 


2 


At (Ee + Gy?) dx dy 


ou, 


E and G are the elastic and shear modulus, respectively. Made 
By substituting Eqs. 30 and 31 into Eq. 29, and a a aa with Tespect 
to the following expression is obtained. > tee: 


in which J, = 4wt, ‘moment inertia of both slabs. By 
Egs. 26, 27, and 32, the total potentialenergy Mis 


U 


+ — (x) dx + 


The potential energy of the load system, oily 1A 

The strain energy of the sidewebs, I1,, is 

iA 

shearing strains ure, respecively, 
| 
| 
4 
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a. E 3w? Val 
If we denote i x, and x,, the ends of the interval of integration and apply 


in whic = moment of inertia of the whole cross section. vex er Bi 4 


ut 
to Eq. 33, the is 


Elz” (x) + M(x) + — EI,U’ (x) 


EI,| — U' (x) +—2" (x) 
To satisfy Eq. 34 for any arbitrary values of 8z” and 8U (positive or negative), 


both terms under the integral must vanish. This condition produces the following 


‘This equation defines the boundary and continuity conditions of U(x ) .) for adjacent — 

=; of different \ ieee At a secti a section where the slab is fixed to the © 


«= is arbitrary. at a section yn where. the slab is not fixed, therefore the bracketed - 
value in Eq. 38, must vanish 


‘EL, U' (x) +— 2”(x) 


expression between brackets in Eq. an ge 


_ These conditions, as in Eqs. 39 and | 40 and the a Alaa are 
in addition to those imposed on 2(x) and M(x) by the simple technical bending — 


7 
> 
pa / 
2 
( 
| 
_ Phe transition conditions Tequire continuity of the function U(x) and of the 
| 


17 710 

By elimination of U(x) and its derivatives from Eqs. 36 and 37, a ‘differential 

equation involving only the derivative of the vertical displacement . of zas an 


The last bracketed term in Eq. 41 represents the s shear lag effect. If shear 
deformability were zero, i.e., G = 0, this equation would reduce to the known 
equation of simple bending. 
ba Eq. 41 has the same form as | obtained by E. Reissner (8) in in which he developed — 
7 for a rectangular doubly symmetrical box section (Fig. 2). Only I, and J nll 
now different values. As by Reissner, this equation can be 


By the follow ing two parameter kandn, where 


(x) = 


constant cocfficients. Using the same paramnetets n and k, the boundary condi- 


for fixed end (slab attached tos support) z 
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following s symbols are used in this paper: q 


Mw = distance of neutral axis from top and bottom of CTOSS 


E, G= Young’ s and shear ‘modulus, respectively; 
= = moment of inertia of whole section, ‘slabs and we 


respectively; 


span length; 


| 
= 


u(x,y), U(x) = spanwise displacement, measure of shear lag, respectively; 
he: = linear strain in in span direction of top and bottom ab, 
5 


re. x,y,z = coordinates; we 


= shear strain of top and bottom slab, respectively; eaSA, ae 
_ = stress in slab at the contact with web at bottom and top; 
= total potential energy, potential energy of load system, 


of the slabs and webs, respectively. 


a 


P= 
om 
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am SEISMIC DAMAGE IN REINFORCED CONCRETE F ened = 
ae y Hooshang Banon,' John M. Biggs,’ M. ASCE and H. Max Irvine* — 


a Damage prediction for reinforced concrete (RC) fi frames is usually measured 2 


of demands for beams columns. In order t to develo op 


accuracy of both mechanical models and damage parameters. 
al be Reinforced concrete members, whee subjected to cyclic inelastic deformations, , 
bor both ‘stiffness and strength degradation. Development and testing of 
§ mechanical models s with satisfactorily | reproduce the behavior of RC members 
_ under cyclic loads i is an important step in the prediction of damage by calculation. — 
Once the accuracy of the model is tested, damage indicators such as damage 


and rotation may be used in a model 


elements e employing. the results of. the is ‘developed 
(2, and will be published elsewhere). == 
When a reinforced concrete frame is subjected to severe ‘snthgsiin’, some > 
of the beams and columns are expected to undergo substantial inelastic actions. 
A mechanical model used to analyze RC frames should be able to reproduce : 
inelastic cyclic deformations of the structure under inelastic cyclic loads. The 
: choice of a mechanical model depends on factors such as expected accuracy © 


§ results and cost of the analysis. Since the e goal of the present study is. to ¥ 
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are introducec and their eifectiveness in prediction Of damage is teste< 
&§ 
&g 
| 

jj 


be able to predict damage v via a set of cia response parameters, special _ 
care has been taken in choosing a mechanical model. The following is a list 
of some of the more popular models which are available but were not a 


“ 
a dual model, each member is replaced by an elastic element 


dual component model i into a multicomponent model by: assuming several r members 
in parallel. Since the end rotations of both components in a dual component — 
model are the same, the element stiffness matrix may be easily formulated. 
_ One advantage of the dual component model is that the yield condition at one — 
av end of the member depends on end rotations at both ends. On the other hand, | a 
since the dual component model can reproduce only bilinear behavior and poe 


any stiffness degradation, it is not suitable for analyzing reinforced concrete 
In a fiber model, each section is divided into many layers or fibers (8,13), 
and the moment-curvature relationship for the section is determined from steel © 
and concrete constitutive laws. _ Member stiffness r may then be determined | by 
— integrating along the member length. Since setting up the stiffness matrix of | 
lt) each member requires many computations at each time step, the fiber model 
- becomes very expensive. Finally, various finite element models have been used 


- to analyze RC structures, but they are also uneconomical for inelastic dynamic 


this ‘study, a component model (6) is to analyze inelastic 


behavior of RC frames. In this model, each member is represented by an elastic — 


beam element with inelastic springs (hinges) at its two ends [Fig - \(a)}. These — 4 
two inelastic springs represent inelastic flexural deformations within the length 


of a member, i.e., inelastic flexural deformations are lumped at the two ends. 


In order to find the characteristics of the two hinges, it is assumed that the 

member deforms in antisymmetric bending [Fig. 1(b)] ‘Since the point of 

contraflexure is in the middle, each half of the member may be viewed as 

a cantilever. By setting the tip displacements of the half-length cantilever and | 

the model equal [Fig. I(c)], stiffness characteristics of each hinge is found. 
In reality, plastic rotation at one end of the member depends o on rotation at al 


characteristics o of the two hinges become dependent. The ; antisymmetric bending 
= assumption is fairly accurate for girders if the effect of dead load is negligible, — 


and it is also accurate for columns in low-rise buildings, = = —— 
Flexural deformations of each hinge are assumed to follow the Takeda model 
@, A modified version of the Takeda model (10) is employed in the present — 

‘study. Fig. 2 depicts the moment-rotation hysteresis curve for the modified 

 _ model. The basic modifications to the original Takeda model are as 
: follows: (1) The hysteresis curve has a bilinear (instead of a trilinear envelope); 


and ag stiffness degradation is built into the model. A complete ——— 


a "There are “two odie sources s of deformations, other on fexerel 
- ‘in a reinforced concrete member subjected to high flexural, shear, and axial 
loads. These two sources are inelastic shear deformations and slippage of main 
_ longitudinal reinforcement. shear loads in a reinforced concrete member 


il 
| 
q | 
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hand, slippage of the main longitudinal reinforcement causes a crack opening 
at the joint interface; i.e., the member starts to rotate at the joint interface. 7 
_ The main effect of the aforementioned two inelastic deformations is a pinching 
of hysteresis loops which is often observed in laboratory experiments: (18). 
If the effects of shear and slippage are minimized in the laboratory, test of 
a specimen, the Takeda model can accurately reproduce the hysteretic behavior — 
of the specimen. On the other hand, it is important to take shear and slippage 4 


deformations into account when they are present. In the present study, a further 
pair of inelastic springs is added at each end of the member for this Purpose. - 
a The characteristics of these springs s may t be found by consi a 
of opening and closing of cracks due to shear and slippage deformations (2). 
Since the mechanism is the same for shear and slippage, both deformations | 
are found to have the same effect on the hysteretic behavior of a member g 


— (15). The proposed hysteresis model for shear and _ slippage is shown in Fig. 


: FIG. 1.—Single Cor Component Model ‘ AG. 2- 
/* This model is completely defined by a set of seven tules, which are shown ; 


Moment-rotation relationship is elastic up to the yield point. 
2. Once the yield point is exceeded, loading proceeds on beh — slope §f 
4 Unloading from rule 2 is parallel to o the elastic stiffness. re 
4. Once the unloading from 3 has ended, there is a marked reduction in 
- gtiffeces; i.e., little load may be applied in the oppostite direction until the 
— is closed. In the present study, the stiffness of this part is assumed 
to be 50% of th the second slope of the bilinear e1 envelope [a result in conformity — 
: g 5,6. If the direction of loading changes while closing the crack, loading or 
unloading will be parallel to the elastic stiffness. «ts 
7. Once the crack is closed, loading will be towards the previous maximum 
rotation point (Point A in in Fig. a a strength-de degrading feature 


; 


also been built t int ot 


A value of = 0.8 is ‘suggested i in this study. i 


a 8. If the direction of moment changes while unloading “ 


It may be the inelastic springs at ‘the end of a member act 
: independently; therefore there is no interaction between flexural, ‘shear, om 
slippage deformations. In reality, both shear and flexural deformations in a 
— 


IG. 3.—Moment-Rotation Hysteresis Model for and 


member are primarily concentrated near the joint; thus | there is interaction between 
these two inelastic deformations. Slippage of the main longitudinal reinforcement _ 
would make the picture even more complicated. Although the simplifying — 

assumption of independent springs is not quite correct, it is both convenient 

-\ and (itis believed) sufficiently accurate. Other studies (9,14), which have employed 


a procedure, have also found the method to be satisfactory. 
? is practice, most RC beam sections are designed to have different amounts" 
of steel at the top and bottom. This means that most sections have different | 
= and yield moments in the positive and negative directions (Fig. 4). 


A simple method of analyzing such sections is to use an average stiffness in 


both directions, which would result in an overestimate of stiffness in the positive 
"directions and an underestimate in the negative direction. Instead, a modification 
to the single component model is made to take into account this effect. Fig. 
5 shows a member which is acted upon by two moments of opposite sign. 


the model. Thus, anew point of is 


4 


‘Ther point of divides ‘the member cantilevers, 
cantilever has different stiffness properties. Let us assume again that, as with 
the single component model, all the inelastic rotation is lumped at the two 
end hinges, and the properties © of these hinges which follow the Takeda model 
are set a priori. The difference is that the elastic element in this case is ‘made 
up of two segments with different stiffness properties. The stiffness matrix _ 
_ of the elastic element is constructed by static condensation of degrees of freedom : 
i at their connecting point (Point C in Fig. 5). A more consistent model is the 
- connected two-cantilever model suggested by Otani (14), where the Takeda model — 
is s applied to > the load d deflection of each cantilever. Unfortunately, the connected - 
_ two- cantilever model results in a a nonsymmetric ‘element stiffness matrix. Although 
‘ ‘the present model fixes the point of contraflexure in the middle to find the 
. properties and then later takes the location of the point of contraflexure | 
into consideration to construct the member stiffness, it has the advantage over 
the connected two-cantilever model of producing a symmetric stiffness matrix. 
it may also be noted that since the length of each segment (Fig. 5) changes 


_ FIG. 4.—Moment-Curvature Reiation- FIG. 5.—(a) Moment Distribution; and 

“ship for Nonsymmetric Section (6) Stiffness Properties for 
metric Section 


at each time step, the element stiffness 1 matrix : hes to be set up at each time | 
_ step. Therefore, this element is not attractive for for dynamic analysis of RC frames _ a 


in 
: terms of ductility demands for individual members. Although much effort has" 
been spent on the computation of ductility in RC frames, there is little me 
available which correlates ductility demand with actual observed damage in 
laboratory experiments. One major problem is that ductility alone may not be 
: a good indicator of damage in RC members. Low- -cycle fatigue-type damage, 
which is caused by a number of inelastic cycles, is also important in failing 
_areinforced concrete member. Therefore, a combination of sudden high deforma- 
tions (e.g., denoted by ductility demand) and cumulative fatigue damage may 


of prediction o of damage in ae frames | through the use of damage indicator: 
| = which are computed from analytical models. It is felt that any mechanical model 
__ which is sufficiently accurate in reproducing the member behavior may be used 
= for this purpose. The emphasis here i is mainly on computing the damage indicators ae : 
4 


of various s damage indicators for reinforced concrete members is s given in the 


3 The most widely used indicator of damage is rotation ductility, which is “a 


‘in wi which 6 in . Fig. 6(a) ri represents total def deformation ; at the end of of a m member. “a 


: 


er. 


» Ductility; 


In a single component. model, rotations (6,) occur at the end hinges. 


| In order ” estimate 6,, it is assumed that the member yields in antisymmetric 


in which L = the member length; and EJ = the member stiffness. The same 
definition of ductility may also be applied to the moment-curvature relationship 


a section. Curvature ductility is simply [see Fig. 6(b)}: 


"FIG. 6.—Detintions of: (a) Rotation Ductility; (b) Curvatur id (c) 
— 
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a bilinear moment-curvature relationship, curvature may be 7 


written as Pigs GOs re Coe 


in 1 which 7 p= = the anf the porary slope t to initial elastic si stiffness an and M... 

; = the maximum moment. This definition of ductility is intended to eliminate 

the need for assuming antisymmetric bending of an element. On the other hand, 

curvature ductility applies only to the most damaged section along a member, 

and it ¢ does n not reflect the overall state of damage for a | member. ‘This definition 

: is even less reliable for the single component model because it does not take 
into account any stiffness or strength degradation, 
_ Another interesting idea for prediction of damage was first proposed by Sozen 
an. Damage ratio (DR) is defined to be the ratio of initial na “en 


to a reduced secant stiffness at maximum ‘displacement [Fig. hae 
mits 


: For a ‘single component model, the definition is applied to each half of a 
q member, in which | b= the tip displacement of each cantilever; and P= the | 
“shear load (see Fig. 1). Although the DR is independent of the yield displacement, 
‘it is directly dependent on estimating the initial tangent stiffness (K,). It is” 
usually difficult to compute the post cracking stiffness of a RC member accurately, ; 
Z shear and slippage deformations greatly reduce the stiffness. This is 


examined in more detail later when comparing experimental and analytical results. 
A modified version of damage ratio | is “Suggested in this this study. Since - flexural 
stiffness « of a cracked member may be ‘computed separately, ee is replaced 7 
_ by this quantity. The Flexural Damage Ratio (FDR) is defined Pic’ satire 
a member which is in K, = 24 The 
ratio in general is expected to be a better predictor ‘of damage than ductility, : 
because the reduced secant stiffness also reflects strength degradation. 77) Lae 
_ All the damage parameters introduced up to this point lack one important — 
feature: they do not reflect the cumulative damage effect on failure of a member. 
It must be realized that low- -cycle fatigue- type damage i is possible under earthquake — 
excitation. Two other "damage parameters are introduced for this purpose. 
_ Normalized cumulative rotation (NCR) is defined as the ratio of the sum of 
al plastic rotations | in inelastic springs, except for | unloading | parts, hae yield 
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‘that is in springs to that time to ‘the maximum elastic 


_ energy that would be stored in a member if it is ete to } antisymmetric: 


F which ¢ = time elapsed since the beginning of loading; and 0(d_) is the 


cy increment at one end of the member during the time interval from | 

™1 + dz. For the present mc model, hinge moments are the same as the member © 
end moment , and the summation “of hinge rc rotations plus | the rotation of the — 
elastic element make up the member end rotation (8). 

7 _ From the foregoing review on damage parameters it may be noted that they 
en. to two general categories. Rotation ductility, | curvature ductility, and 


and dissipated energy ‘reflect ‘cumulative fatigue- ype ‘damage which is “due to 
a significant number of inelastic cycles. Any method of predicting damage in 
members should include one parameter from each category. The 
- damage indicators would then depend on whether one is more informative : 
tha another, and on the accuracy of the computed comnge indicators from 
— Choice of Biiectienitil Specimens.—In or: order to test t the accuracy « of mechanical 
| models and determine the effectiveness of damage indicators in predicting damage, © 
it was decided to analyze a sample of laboratory tests of RC frames. Although 
the results of many dy namic (shaking table) tests of specimens | have become 
available in recent years, they do not directly provide the type of information 
_ which would lead to a model of damage in reinforced concrete members. 
Therefore, for the purpose of this study, it is found appropriate to limit our 
- to available static test results in the literature. Changes in stiffness 


and strength properties of a specimen may be directly observed in a quasi- “static: 


ode load test. Once some kind of damage (or failure) criteria is established, 
it is also possible to identify the time of failure in a static cyclic load test. 
Another advantage is that loading in a static cyclic test may be closely monitored. 
The sample in this study includes laboratory tests on 32 specimens from 
eight different sets of experiments (1,3,5,7,12,13,17,20). All of these tests are 
on large-scale members or subassemblages. Special care was taken in choosing \ 
the > experiments to ensure that they r | Tepresent ina reasonably realistic manner 
4 the behavior of RC buildings under earthquake excitations. It should be kept 
in mind that inertia and damping forces are absent in quasi-static tests, and 
the rate of loading is also different from a real seismic loading. Since there 
is no conclusive evidence on the effect of rate of loading on the failure of 
a member, the results are expected to apply to dynamic cases as well. Poy 
There are three main objectives in comparing analytical and experimental 
results of the sample of cyclic load experiments. First, it is determined whether | 


the mechanical model, which is the single i tcianes model in this this — is 


| 
| 
| 
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sufficiently accurate. Second, by monitoring the computed damage owl 


fo the analytical model, a relationship to damage may be established. At 
: - this point, it is also possible to compare the computed values of damage indicators 


with those _ measured in the test. Third, by choosing a criteria for the failure 

of a member, experimental failure points are determined. The failure point marks _ 

the end of the analysis, at which point all of the damage indicators at failure 


measure re of in 1 buildings. There a are for 
post earthquake damage into different categories (21), but these schemes do 
not apply in this case. Instead, it is decided to define only the point of excessive _ 
damage or failure for each specimen. Therefore, the objective is to predict 4 
; failure of each member (local — in a frame, given t that it has undergone. 
certain inelastic deformations. pate! | the 
Failure of a member in the seed study is assumed to be an instantaneous _ 
_ phenomenon. This is also true of many laboratory tests of specimens that are 
- considered in this study, i.e., the load-carrying capacity of some members starts 
: to drop off very fast. In such instances, it is fairly easy to define the point © ; 
of failure. In other cases, failure of the member is gradual, e. g., the member 
can go through a few more ‘cycles of loading although the capacity of the member 
to sustain loads drops off with each additional cycle. For members with gradual | 
failure, load-carrying capacity is used as the failure criteria. If there is more 
than 20% difference between analytical and experimental load at maximum 
displacement in a cycle, the member is assumed to have failed. It is assumed 
_ here that the analytical model can reproduce the behavior of a member until — 
a event, such as buckling of main longitudinal reinforcement, causes failure 
of the member—at which point the model becomes ineffective. In practice, 
although the member may be able to carry some load beyond its theoretical 
failure point, _ it is best to assume that the member has failed. = | 
- Accumulation of damage in RC members may be modeled through the use 
4 of damage indicators. Values of damage indicators are monitored throughout 
‘the test for all the specimens considered in this study. The failure point for 
a member, as defined previously, is considered the terminal point of a test, 
and the objective is to predict the failure point via the knowledge of damage © 
indicators. One may note ‘that any other logical definition of failure could also 
be used for this purpose. However, it is felt that other definitions - of ‘failure 
_ will not change the final values of damage indicators appreciably. he: ey {on 
Comparison of Theory and Experiments. —The sample of laboratory experiments 
which is gathered in the present study covers a wide variety of loading conditions _ 
on RC members. As mentioned before, the sample consists of 32 experiments — 
‘ - eight different sets of tests. Although the sample size is rather small, 
it is felt that general conclusions and recommendations would not be ‘mage 


by including more data in the sample. On the other hand, for a . model of damage 
in RC frames, a large sample size would result in a more accurate and reliable - 
_ model. The following is a comparison of the experimental and analytical ae ; 
some individual cases and the sample asa whole. 
Fig. 7 is a comparison of the hysteretic behavior of specimen 8 (A8) tested 
by Atalay and Penzien (1). These sets of tests were designed to study | the — 


Gatenier of reninguet concrete columns under high axial and flexural loads. 
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‘Thet test set- -up ) for the | experiment is also in Fig. 7. The tests. specimen 
= 


represents two columns of a high-rise building between their inflection points, 
a which are assumed to be at their midheights. It is seen from Fig. 7 that the 
test set-up does not simulate double curvature in columns; however, this type | 
of test set-up is used for simplicity. The test procedure consisted of yoy 
an axial load on the specimen, and then applying predetermined cyclic displace- 
ments on the joint. All specimens in this set of tests had cross sections of — 
12 in. x 12 in. and a length of 11 ft, and for specimen 8 in Fig. 7 the axial 
load was 120 —. Fig. 7(a) and 7(b) are the experimental and analytical 


da bispuacement 
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FIG. 1 —te) and ( (b) Ang Analytical Le Load-De' Deflection Curves for 
A8Tested by AtalayandPenzien 


load-deflection curves for the specimen. The latter was computed using the 
model described previously. The secondary effect of axial load is clearly seen 
in these two figures, e.g., the load starts decreasing as the displacements increase. _ 
_ Hysteresis ra for this specimen are very stable up to the failure point. Failure” 


There are two aspects of the inelastic cyclic behavior which are important _ 
in this type of test, namely, peak loads and the energy dissipation. A visual 
inspection of Fig. 7(a) and 7(b) reveals that the model is fairly accurate in 
both loads | the energy | for this” experiment. Accurate 
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secant stiffness i in DR and FDR is 


ts ment. In fact, the displacement ductility, which is the same as rotation éuctility q 
in this case, is measured to be 5.4 and 3.6 from the analytical model and the 
experiment respectively. On the other hand, since the experiment is carried * 

_ out well beyond the yield displacement, overestimating the initial tangent s stiffness 

a of the member has little effect on the overall inelastic behavior of the specimen. bale 
It may also be noted that the analytical model is fairly accurate in this case, ' 7 


because shear and slippage deformations are not present in this set of experiments. 
In fact, the main longitudinal bars were welded to steel plates inside the joint 
to minimize the effect of slippage, and due | to rather long span-t to-depth ratio, 


study the single component model employing the Takeda model 
accurate in predicting the inelastic behavior of RC members under flexure. bend 


Analytical and experimental load-deflection curves for specimen 6 (S6) in 


the experiments by Scribner and Wight (17) are shown in Fig. § 8. The test set- = 


uw 


for the ‘specimen is also depicted on ‘the ‘same figure. The specimen "represents 
an: exterior joint, with two columns and a beam in between their inflection 
points. This specimen is a half-scale model of a joint in a tall building. ‘After 
application of axial load to the columns, predetermined displacements ‘were | 
applied at the tip of the beam, and the necessary loads were measured. The 
loading scheme consisted of six cycles of displacement ductility six in the positive 
direction and ductility four in the negative direction. If the specimen survived 
these cycles, more cycles of higher ductilities were also applied. Since the 
inelastic action is intended to be confined t to the oem, all columns were designed + 
_ Pinching of hysteresis loops due to the effects of shear and slippage i is clearly i 
seen in Fig. 8. Although the computed yield displacement is much lower than 
the actual yield displacement, the overall match between the experiment and vs ; 
is very good, €.8., normalized dissipated energy is computed to ) be 


"experiment. gradual failure i: is observed for this spe specimen. Strength degradation 
. during cycles of equal ductility may be attributed to the presence of shear _ 
and slippage deformations in this specimen. It may = - noted that, using 


prediction of load also implies that the 
accurately computed. On the other hand, IT is observed the analytical yie 


a 


point is reached one- -and- -a half cycles before the is terminated. 
The definition of failure is expected to affect dissipated energy and i 
cumulative rotation at failure, but in most cases it has little effect on computed = 
Fig: 9 shows the load-deflection curves for specimen 10 (S10) in the a 
set of experiments (18). The difference in the two specimens is only in their ° 
size and steel reinforcement; specimen 10 is a full-scale model. Pinching of — 
hysteresis loops is also present for this specimen. On the other hand, a comparison — 
of Figs. 9(a) and 9(b) reveals that the computed loads at peak displacements 
of each cycle are generally lower than the actual loads measured in the test. 
This has } caused the analytical 1 normalized dissipated energy (197) to be Toughly — 


exhibits a sudden drop in its load-carrying capacity beyond the seventh cycle 
of loading. The cea ersse: point is defined to be at zero wrnpercenerl 


1 


7 S10 Tested by Scribner and Wight foe 


Table 1 and “experimental. energies for the ten 
; specimens in the tests by Scribner and Wight which are included in this study. 
_ Specimens 4-8 are half-scale models, whereas specimens 9-12 are full-scale 
models; all represent an exterior joint in a tall building. Dissipated energies 
listed in Table | all represent values up to the theoretical failure points. Dissipated 
energies compare well for the first six specimens, but much. lower dissipated — 
energies are .re computed from the model for the last four specimens. The difference, 
as in the case of specimen 10, is due to underestimation of the yield load 
: by the analytical model. Unfortunately, dissipated energy is rarely mensused 
quasi-static cyclic load tests. For other specimens in the sample where energies 
- were measured in experiments, analytical and experimental results compared 
Based on the study of analytical and experimental hysteretic 

_ the sample of RC specimens, the following conclusions are reached: 


1. The single component model is sufficiently accurate in por 

behavior of RC members (see Figs. 7, 8, 9). There is especially good agreement 

between the model and capectmente when the: flexural mode of deformation 


| 
| 
Load-Deflection Curves for Specimens 
a 
i 
| 


is dominant and shear and slippage "deformations are minimized. This proves 
that the Takeda model can Benue reproduce deformations of RC members 
under cyclic loads. RC members. Ductility HY 
4a 2. Rotation ductility, which is a popular measure of member damage in RC 
_ frames, is generally computed to be higher than actual experimental ductilities. = 
- This is due to the fact that the elastic stiffness of a RC member is reduced 
by shear and slippage effects; therefore, the yield rotation is underestimated. 
It is also found that the overall inelastic behavior of RC members is insensitive 7 
_ to the initial elastic stiffness, i.e., , hysteretic behavior and stiffness degradation a 
beyond the yield level are little affected by the elastic stiffness of the member. 
— SY Curvature ductility in a single component model can be computed only 
by assuming a bilinear moment- -curvature relationship (Eq. 4). On the other 
hand, using stiffness- characteristics is obviously in conflict 


in in RC members. Thy pats snd 
TABLE 1.—Comparison of Experimental and Analytical Normalized en ~— 


te 

6.01 


490 


damage. Since there was agreement hetwoen and 

experimental loads at peak displacements of each cycle, flexural damage ratios 

‘ = for the sample are expected to be accurate. The only other source 
in FDR (Eq. 6) comes from flexural stiffness (K,), but flexural 


_ 5. Dissipated energy (£,,) and normalized cumulative rotation (NCR) reflect 


_ cumulative fatigue-type damage in RC members. For experiments where energies 
_ were listed, generally good | agreement was found between analytical and experi- 


Table 2 lists the five damage indicators, computed at failure from the be ancient i 
‘model, for all of the members in the sample. Unfortunately, all of the computed - 
damage indic show large is in part due to intrinsic differences 


| 
| 
| 


in member resistances. If one of th the indicators at failure, 
rq such as p,, were ‘stable, it would be very easy to predict the failure of RC | 


TA E 2.—Damage Indicators | for Specimens Tested in Laboratory sted 


Specimen _ 
= R3 


H9 


T3 = Ma, et al. 
-©H7-H9 = Hanso~, Conner (7). 


= VBC3 = Bertero, etal(3). 
*FI-F4 = Fenwick, Irvine (5). OM 


4 ‘P43 = Popov, et al. (16). Wied 
*§3-S12 = Scribner, Wight (17). 


roe mo failure of RC ‘members has to be predicted via a 


probabilistic model. In order to construct a damage model for RC members, — 


> 
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SEISMIC DA DAMAGE 
one would like to use as few and as informative Pore iat as possible. a 


. As mentioned before, the set of damage indicators presented in this study represent ‘7 
two different types of damage in RC members. Ductility and damage ratios 4 


E indicate a suddenly inflicted damage due to large deformations, whereas energy 
cumulative rotation represent cumulative fatigue-type damage. A high degree 
_ correlation (0.95) is found between rotation ductility (u,) and FDR for the 
sample. From the foregoing review, FDR seems to be superior to 1, because 
it also takes into account strength degradation in a member. Also, a high degree © 
of correlation (0.98) is found between normalized dissipated energy (E£,,) and — 
normalized cumulative rotation (NCR), , which i is to a large extent a direct result 
of the hysteretic behavior in the analytical model. Quite arbitrarily, itisdecided [7% 

=: Using E,, and FDR as damage state parameters, one may set up a probabilistic 
model of failure i in RC members (2). Damage may be thought of as ei nd . 

7 process which causes deterioration of a RC member until conditions are reached 

_ for the failure to occur. The analytical model gives information on E, and 
FDR as a function of ome. One is interested in the probability of failure a 

_a member at any time, given the knowledge of load history on the FDR-E, : 

_ plane. The probability may be computed if one has a knowledge of the so-calle d 4 
“hazard” or “risk” function on the E, -FDR plane. Fora one-dimensional control 
parameter T (e.g., E,, or FDR), with f(t) being the probability density function 
of T to failure and F(t) = f(x) drt, the hazard (risk) function, A(t), is defined 


in which ‘the prime denotes differentiation. ‘Conversely a 


Clearly, A(t)dt is the probability of failure before T + dt, given pe at 
- time t, and F(t) is the probability that failure occurs at or before ¢. The present 
case is more complicated, because the notion of the hazard function has to — 
be extended to two dimensions. Such an approach has been used in Ref. 2. 
Results of the laboratory experiments are then used to estimate parameters 
> of the probabilistic model, It should be noted that any probabilistic model of 
failure in RC members will incorporate both the uncertainty in 1 member resistance ~ 


— and structural properties, and also modeling errors in calculation of damage > 


sail model, with account ‘shear 
deformations and slippage of longitudinal reinforcement, has been used to analyze ~ 

a set of quasi-static cyclic load tests. Rotation ductility, curvature ductility, : 
flexural damage ratio, dissipated energy, and cumulative rotation are computed — 
from the model for each one of the Iti is that the 
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eycli lic behavior of RC members. The 1 match and 


. results is especially good when flexural deformations are dominant. Based on — 
i accuracy of results and physical grounds, flexural damage ratio (FDR) and 
_ dissipated energy (E,,) are chosen as damage state parameters | for a probabilistic — 7 
model of member resistance. The FDR and E,, represent damage due to large - 
deformations and cumulative fatigue-type damage respectively. On this basis, 
it is believed that damage can be predicted from analytical results with some 


_ This paper is based on part of the first writer’ s Sc.D. thesis done under 
7 ‘the supervision of the other two. The support of the National Science Foundation a 
under Grant ENV 77-14174 is gratefully 
: 1. Atalay, M. B., and Penzien, J., ‘‘The Seismic Behavior of Critical Regions ¢ of Reinforced 
Concrete Components as Influenced by Moment, Shear, and Axial Force,”’ Earthquake 
_ Engineering Research Council, University of California, Berkeley, Calif., , Dec., 1975. o- 
Banon, H., ‘*Prediction of Seismic Damage in Reinforced Concrete Frames,” Publica- 
4 ? tion R80-16, Department of Civil Engineering, M. I. T., Cambridge, Mass., May, 1980. ’ 
Bertero, V. V., Popov, E. P., and Wang, T. T., ‘‘Hysteretic Behavior of Reinforced 
Concrete Flexural Members with Special Web Reinforcement, Earthquake Engineering 
7 Research Council, University of California, Berkeley, Calif., Aug., 1974. ‘ela : 
4. Clough, R. W., Benuska, K. L., and Wilson, E. L., “Inelastic Earthquake Response 
we of Tall Buildings, ”* Proceedings of the Third World Conference on Earthquake Engineer- 
s Fenwick, R. C., and Irvine, H. M., “Reinforced Concrete Beam-Column Joints for 
Seismic Loading,”’ Bulletin of the New Zealand National Society for Earthquake 
__ Engineering, 1977, Pt. I, 10, 3, 121-128; Pt. I, 10, 4, 174-185. | ee 
a. Giberson, M. F., ‘Two Nonlinear Beams with Definitions of Ductility,”’ Journal 
of the Structural Division, ASCE, Vol. 95, No. ‘$82, Proc. ‘Paper 63 saeies Feb., , 1969, 
. Hanson, N. W., and Conner, H. W., “‘Tests of Reinforced Concrete Beam-Column 
Joints under Simulated Seismic Loading,” Research we Development Bulletin RDD12, — 
Portland Cement Association, Portland, Oreg., 1972. 
8. Kent, D. C., and Park, R., ‘Flexural Members with Confined Concrete,” Journal — 
a the Structural Division, "ASCE, Vol. 97, No ST7, Proc. Paper 8243, atone 
9. Kiistii, C., and Bouwkamp, J. G., “Behavior of Reinforced Concrete Deep Beam- _ 
Column Subassemblages under Cyclic Loads,’’ Earthquake Engineering Research — > 
Council, University of California, Berkeley, Calif., May, 1975. = | 
10. Litton, R. W., “tA Contribution to the Analysis of Concrete Structures under Cyclic 
Loading,” >i presented to the University of California, at Berkeley, Calif., 1975, a 
~ . partial fulfillment of the requirements for the degree of Doctor of Philosophy. 
. Lybas, J. M., and Sozen, M. A., “Effect of Beam Strength and Stiffness on Dynamic = 
7 Behavior of Reinforced Concrete Coupled Walls,’’ Department of Civil Engineering, — 
of Illinois at Urbana-Champaign, IIl., july, 
2 Ma, S. H., Bertero, V. V., and Popov, E. P., “‘Experimental and Analytical Studies 
on the Hysteretic Behavior of Reinforced Concrete Rectangular and T-Beams,”’ 
‘ _ Earthquake Engineering Research Council, University of California, Berkeley, Calif., at 


Mark, K. M. S., “Nonlinear Dynamic Response of Reinforced Concrete Frames, 


thesis to the Massachusetts Institute of Technology, at Cambridge, Mass. 


= 
= 


, and Sozen, M. i 
ar During Earthquakes,” * Department of Civil Engineering, University of Illinois . 
_ Urbana-Champaign, Ill., Nov.,1972,00 
1s, Otani, S., ‘Nonlinear Dynamic Analysis of 2-D Reinforced Concrete Building Struc- 
tures,’ Proceedings of the Third Canadian Conference on Earthquake 
16. Popov, E. P., Bertero, V. V., and Krawinkler, H., “Cyclic Behavior of Three R.C 
Flexural Members with High Shear,” Earthquake Engineering Research Council, 
University of California, Berkeley, Calif., Oct., 1972. 
17. Scribner, C. F., and Wight, J. K., “Delaying Shear Strength Decay in Reinforced 
_ Concrete Flexural Members Under Large Load Reversals,’’ Department of Civil — 
Engineering, the University of Michigan, Ann Arbor, Mich., May, 1978. 
: 18. Takayanagi, T., and Schnobrich, W. C., ‘Nonlinear Analysis of Coupled Wall 
Earthquake Engineering and Structural Dynamics, Jan., 1979. 
19. Takeda, T., Sozen, M. A., and Nielson, N. N., “Reinforced Concrete Response 
i to Simulated Earthquakes,’’ Journal of the Structural Division, ASCE, Vol. 96, No. _ 
ST12, Proc. Paper 7759, Dec., 1970, pp. 2557-2573. 
20. _Viwathanatepa, S., Popov, E. P., and Bertero, V. V., ‘“‘Seismic Behavior of Reinforced 
Concrete Interior Beam-Column Subassemblages,”’ Earthquake Engineering Research 
Council, University of California, Berkeley, Calif., June, 1979. 
21. Whitaman, R. V., Reed, J. W., and Hong, S. T., “‘Earthquake Damage Probability 
a Matrices, * Proceedings of the Sth Wor World C d Conference on atic Engineering, P Paper 


Appenpix .—NorTATION 


are used in this paper; 


E. = normalized energy; 


cracked member stiffness; 
probability density function at 1; 
cumulative density — at t; 
flexural 
reduced secant stiffness; 
member length; 
maximum moment; 
moment at time, 7; 
-tatio of second slope to initial slope; 
parameter of shear or slippage t hysteresis model; 4 


= rotation ductility; 
maximum curvature; 
= plastic curvature; 
= rotation increment time TtOT +: 
yield rotation; and 


plastic r rotation. 
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Air- -supported enclosures have been in use only since World War II but t they ig 


are now a common sight in most cities, where they cover, e.g., tennis courts, 
offices, exhibitions, and greenhouses. However, the building codes for these 

_ membrane structures are, in most western countries, not yet completed partly 
due to the lack of documented case histories, and partly due to the continuing 

research into the structural behavior of these forms. 
“4 One structural aspect that has been given only limited attention is the — 

7 - of collapse through accumulation of rain, ice or snow. Failure in this mode 


has not been widely recognized (3)* but it has been a ‘contributing factor in 


_ several failures known to the writers and is a possibility | of which igners 


z Previous work in this field has been limited ve cylindrical inflatables 
(), this being the most common geometry and the one that lends itself to a 
the simplest analysis. Experimental models have indicated the basic correctness | 
of the predictions of the theory for cylindrical structures (2). The next most 
~ common inflated form is the part-sphere and is the simplest shape having double 
curvature. ‘Interesting: investigations into the axisymmetric deformation and 
wrinkling of such a shape have been presented (4) but the present paper is, 
to the writers’ knowledge, the first to present theoretical and experimental results 
on the ponding stability of these domes. 
_ The following section outlines the assumptions and theoretical development 
7 of the governing equations for both the uniform membrane and for a system — 
of one-way acting cables. It also shows how the ‘influence of an irregularity — 
in the form of an axisymmetric bubble may be considered. The section entitled 
_ “Numerical Solution”’ gives a brief description of the solution algorithm and 
computer program together with some relevant numerical results. In the section - 
"Sr. Engr., DAF Indal Ltd., 3570 Hawkstone Road, Mississauga, Ontario, LSC 2V8. 
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This paper is part of the Journal of the Structural Division, Proceedings of the American 


_ Society of Civil Engineers, @ASCE, Vol. 107, No. _ ST9, September, 1981. ISSN 0044 


= 


| 
| 


on experimental models, the experimental tests are we Guede and 1 results 


Its 
these tests are compared with those from the theoretical model 


the initial geometry is a perfect sphere of radius R, the material is. inextensi 

and weightless, and that the deformations are axisymmetric. Fig. lis a cross. 

~ section through the central axis showing the initial and deformed outlines. In — 

the deformed state the central portions of the cap are deflected downwards © 
and result in circumferential wrinkling due to the assumed inextensibility. This 
_ wrinkled zone, in which the membrane acts in the meridional direction only, — 

_ and in which the radial tension per radian substended must remain constant, _ 
extends as far as the “transition level.” The location of this level depends a 


«FIG. —Cross Section through Sphericai Cap with Central Pond 


_ Noting that at the edge of the pond the membrane is horizontal, the overall - 
vertical equilibrium of the cap within the transition zone is described by dt 


in which r = radius; d = diameter of the pond; p = internal over-pressure; 


7 T = radial tension per radian substended; and @ = angle of membrane with 
the vertical. For continuity of slope at the transition | level, cos 2.2 
Eq may be used to solve for the transition level radius, 


At any other point within the transition zone, defined by : a dettes be, 
1 may be used to write write the angle of as 


24208 


‘Writing an element of membrane length as as ds = dr /sin o, Eq. 3 n may be used 
to describe the total length of membrane in the transition zone as 


st9 
development _is based on the assumptions that ss & 
) 
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7 tt addition there i 
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differential equation to the within the 


of derivatives of y as 


1+ {— 


which may be numerically integrated to (half) of 


pond to be be obtained from 
to assumed inextensibility the total meridional unchanged 


S an equilibrium | condition to be stishe at the center 


a however, such irregularities are often encountered when fabrication is 
nd this development the irregularity is regarded as a part spherical cap of © 
radius, nw. which is less than the radius, R, of the remaining cap and has 
a span, 2/’, at its junction with the remaining cap (see Fig. 10). When a central - 
i vertical load is first applied, the deformation is that defined by a cap of radius 
R’ but as the deformation and pond diameter increase to engulf the smaller 7 
cap the problem reverts to the original one with the single exception of the — 
(inextensible) length between the transition level and the centre now being slightly 


longer. Eq. 8 ‘must now be adjusted to read a 


density of ponding medium, and x and = coordinates onginating 
a This equation may be written in terms J 
: 
+ 8, = Rsin" 
| the cap where aC 
cual 
th 
— Irr metric irregularities can be admitted 
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"One-Way Action. the membrane is s made up ofa number of 
a 

cables between which the membrane spans, then a description of the ponding — 
deformation of a structure follows very similar } to that described 


a section - FROM PHOTOGRAPHS 

«FIG. 2 —Outline « of 1. 43 | m Radius Dome 


3. Results of Central Weight Pond Diameter 


two-way action of the membrane, and the unknown, may be 


the known (half) span, /, of the spherical cap. Eq. 4 therefore becomes 


while Eqs. 


uired to determine the ‘maximum value of the central weight; 


ire st9 
A Iterrnatively the total 
| 
pe 
q 
4 
req 


the | is slow and that the internal will remain ‘constant. 
To use Eq. 9 to calculate W the value of the meridional tension, 7, must 
be known. For a chosen value of pond diameter, d, T may be obtained by 
an iterative scheme (such as Newton-Raphson) which attempts to satisfy the 
condition of constant meridional length ih represented b by Eq. 8 or modified as — 


in in Eq. 10 to account for irregularities. tos rarer. 
i The membrane length within the pond, S, , is obtained by a numerical integration — > 
2 scheme after the values of dy/dx have been determined by solution of the 


differential equation, Eq. 6, using a Runga-Kutta fourth order algorithm. The 


"integration (using a Romberg integration algorithm) of Eq. 4 and in the case 
of the one-way oe si membrane te limit of i f integration 7 extended to include 


the entire cap. 
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CRITICAL CENTRAL WEIGHT, W, Newtons 


10000 20000 


_ The typical a of the numerical results obtained is ve in Fig. 3, which 
= demonstrates that a limiting. value of the central | weight does exist and that 

‘ if the pond diameter is increased further then the central weight must be decreased. — 
$ Fig. 4 shows the effect of the internal pressure and radius (with fluid density 
that of water) on the critical central weights, and indicates the greater susceptibility " . 
“a theory) of structures with small radii of curvature. ‘Fig. 5 plots the influence > 


of fluid density and radius on the critical « central weight when the mt 


overpressure is kept constant at 200 Pa (0.81 in. of water). ‘It is not surprising 
to learn that fluids with higher densities pose the more serious hazard. Bis 
_ The influence of irregularities in the initially spherical dome is demonstrated 
in Fig. 6. The extent of the axisymmetric irregularities is not given in the form 
_ Of radius and span of a central bubble but instead is described as a percentage — 
increase (over the perfect sphere) in the total membrane length within the ponding 


_ and transition zones. This would correspond to a series of axisymmetric ripples ae 
in the initial membrane. Fig. 6 shows that an increase in the meridional aaa 


| 
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‘ 200 q 
Of Pressure Radius FIG. 5.—Effect of Fluid Density on Criti- | 
on Critical Central Weight calCentraiWeight 
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= structure may therefore be described as sensitive in 
=z Fig. 7 gives results for the dome having one-way meridional action only : and ; 
_ shows the influence of pressure and span on the critical central weight. It is 
important to note that such a structure is not initially a sphere but is flattened 
at the apex and, strictly, should not be included under the title of this myer. : 


ot ad fluid density (water) = 9600 Pa Ip 
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A 6.—Effect of Irregularities on Critical Central Weight 
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system sal cables onto a spherical cap and to realize that the meridional tensions - 


are now functions functions of of the span as well as initial radius. 

— 

Fabrication. —Tests were carried out on two air-supported spherical caps, 
sz made from axisymmetric segments of a polyester fibre reinforced PVC 
--membrane (approx 6.8 N/m’, or 0.14 lb/ft”). The first was already available, 
_ having been initially used for other paapeats and had a a cadies of 1. 43 n m (4.69. 
7) with a total span of 2. 34 m 


| 
ma 


ol, 


7 »* The first model was made f rom m 36 segments whereas the second dome c one 7 
= ¥ only. The smaller number was considered to approximate a sphere accurately | 
enough and to offer less opportunity for error in fabrication. A template for 
the segments for the second model was made from hardboard and a 1/2-in. 
/(12-mm) overlap v was allowed for bonding. This | bonding was achieved using 
prior treatment of surfaces with a solvent containing equal volumes of isopropyl 
_ alcohol and tolenol, followed by A4924 adhesive and clamping for a minimum > 
of 20h. A specially made oa was used to carry out this clamping [see _ 


inches of 
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’ FIG. 7.—Effect of Pressure and Span on Critical Central Weight—One- Way | 


% Both models included an air- -tight floor and were secured around their perimeter 7 
te a 3/4 inch plywood base. The central load was applied through a 1/4 in. 
diameter steel plunger (weighing 4.12 N, 0.93 lbs), which was contained by 

a guide located in the center of a steel member spanning above the dome. 
The plunger rod supported a platform at its upper end to which additional 

_ weights could be added and a pointer was ; attached so that the central deflection 


could be read off a scale connected to na cross beam. Some of these details 


_ Both short and long-term ‘tensile loading tests were carried out on 1 material a 
samples a as well as on test thew 
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FIG. 9.—Typical Test Results of Pressure Versus Central ian 
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.—In eo simulate one-way meridional action, a set of eight — 
"steel cables, joined at the apex by a 80-mm (3.1-in.) diam ring, was placed © 
over the first (1. 43, m radius) dome. Each cable was secured to the plywood i, 


- of bonded strips to ensure that the membrane and cable system deflected together. » 
of the cable system are shown in Fig. 8(c). 
3 Testing Procedure.—The domes were inflated by means of an electrically 
driven portable compressor and the overpressure was measured by a manometer 

_ Although the most common ponding medium available to full-size structures 
is water, it was not possible to contain such a fluid without considerable spillage, 7 
_ due in part to the irregular furrows that the deformation caused in the membrane. S 
; Therefore a a fine quality, Ottawa sand was used, with a density of 14,300 N/m° | 
Ql ib/ft°). This allowed exact measurement of the ponding medium to be 
- made and levelling of the medium could be carried out at lesiure. nd Wh setae 
To simulate the pond diameter versus central weight relationships shown in 
Fig. 3, it would have been necessary to maintain a constant nent Since 


44, 


10 —lrregularity rity as  Axisymmetrc P Part- of Radius 

this could not be done accurately, the denned weight was instead kept constant 

while the pressure was reduced and the ponding fluid maintained at the appropriate 

level. During this operation the central deflection was noted and a graph of 


a pressure versus deflection was obtained as depicted in Fig. 9. The critical pressure 7 
for each value of the central weight was obtained the 


main difficulties during the testing were: to imperfec- 
tions in the shape of the domes the poad formed was not circular; (2) imperfections 
in the membrane shape resulted in difficulty in judging the correct level of — 
_ ponding medium; and (3) small air leaks in the first (1.43 m radius) dome limited 
“the pressures under which that dome could be tested. Nevertheless all results 
obtained were reproducible. At values of the central weight lower than the 
self-weight of the plunger (4.12 N), small dead weights were placed at the 
apex of the dome. 
& The initial shape of each | dome was s examined photographically to obtain an : 
. 10 shows a tracing > 
of the outline of the first (1.43 m radius) dome and comparison with a circular > 
outline. It is difficult, however, to quantify the extent a — pe at 
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to theoretical predictions for and for with a 2% 
in meridional length in the displaced zones. Fig. (11) also shows theoretical ; 
and experimental results for the one-way action on a sphere of initial radius 
1.43 m. The eight radial cables placed on this dome were tightened by means 
of turn-| -buckles until further tightening made no difference to the experimental ; 


results. (These cables are visible in Fig. 8(b) and (c).) However, as may be 
noted from Fig. 11, the presence of these me made little difference to the 


Asystem of meridional cables was also on the second, 0.8 m, dome. 
This time the strips which loosely secured the cables to the membrane were 


omitted and the result was that the the weight o of the ponding medium ——— 


FIG. 11 Results from FIG . 12.—Experimental Results from 


the membrane, separating from the the cables the and 


_ The experimental results differ considerably from the theoretically predicted © 
values for perfect spheres and suggest that failure may occur at values of central. 
load far less than expected. The only satisfactory way of explaining this apparent — 
discrepancy is to suppose that the sphere is not perfect since, indeed, the | 
predictions are extremely sensitive to small increases in the meridional length 7 
available to form the pond. This small increase in length may be interpreted _ 
as an equivalent decrease in radius, which may be shown by calculation | to 
‘No adjustment has been made for the self-weight of the membrane which 
was measured as 8 N/m’. However, near the apex where the membrane — 
overlapped considerably, the self-weight may have averaged about 14 N/m’. 


results should, therefore, especially at lower pressures and 
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. 4 smaller pond : sizes, be corrected accordingly. This a the effect of utente 
somewhat better agreement with predicted values.§ 
_ The membranes were assumed inextensible and material tests carried a 
supported this assumption within the range of membrane loads involved in these 
experiments. However, it is expected that a theory which includes the elastic 
deformation of the membrane will require higher central loads to cause failure 
_ The experimental results obtained using a cable system on the 1.43 m dome — 
give critical weights far lower than expected. The reasons for this are firstly 
a that the wide spacing of cables (only eight were used) did not prevent considerable — 
4 circumferential stresses from remaining; ‘secondly, the flattening of the surface 
in the region of the apex, caused by the one-way action of the cables, also 
_ produced considerable wrinkling of the membrane and these wrinkles collected ns ; 
the ponding sand and hence increased the deflection due to ponding. Indeed _ 
this wrinkling also occurred to some extent on the domes without cables and “y 
‘no doubt « contributed i in the: same way tothe increased displacement an and premature * 
collapse in those experiments. The ; experience with the cable system on the - 
q 0.8 m dome, in which the membrane was free to separate inwards away from a 


the cables, suggest that this is the most dangerous situation of all despite the 
theoretical prediction of added stability due to increased meridional tensions. — 
Indeed the theoretical results are based largely on the meridional tension 
that can be developed and therefore predict that low profile pneumatic structures _ 
will have greater stability against ponding than will more curved shapes. However — 
the theory also points to the importance of the concentrated central weight | 
in initiating ponding and it is the flatter surfaces that are less able to shed 
- external loads in the form of lumps of ice and snow. Unfortunately this type 
a: concentrated loading, unlike an internal load applied t toa cable intersection, 7 
e. g. is difficult to estimate a and it is only observation of full scale ¢ structures 
7 ir. poor weather and under low pressure that will indicate the magnitude of 


the discrepancy between theoretical and experimental models prevent any firm 
design tules from being formulated, the following recommendations are mate: u 

A It is possible ‘that the danger c of collapse by p ponding i is even greater r than 

theory predicts for surfaces cf double curvature. gar nae 
a + These structures are extremely sensitive to surface irregularities near bd 


3. Further model tests should be m in exact control or measurement 


the surface geometry is included. hey — 
= Observations should be made of full scale inflatables (of either single or 
double curvature) to determine the likely accumulation of snow and ice which 

et . The exact interaction between cables and membrane under ponding condi- 
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a The following symbols a are used i in 


span of part-sphere; 
internal over-pressure within inflatable; 
radial distance from central axis. of symmetry; 
radial distance to transition level; 
= meridional distance along deformed ‘membrane; 


_T = meridional tension per radian subtended; 


= diameter of ponding medium; 


central point weight; 

= horizontal coordinate from edge of pond; 
vertical coordinate from edge of pond; 

angle between meridian and vertical; and vax 
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OPTIMIZATION I PROCEDURE FOR PULSE- -SIMULATED 


Sami F. Masri‘ and Frederick B. Safford? Bol leubividn dass 
on od). basen of lo mie ot 


“ | ‘The capability for simulating the response of structures, in-place and as-built, 
- to transient dynamic loadings (such as earthquakes and blast loads), is useful 
for testing structural adequacy, for improving mathematical models, and for | 
_ investigating the response of equipment in a structure (9). In addition to various 
types of vibration generators that are appropriate for certain classes of structural 
systems, large testing facilities (which have limited availability) and ground- 
explosion approaches (which are economically and practically limited) can be 


used for dynamic tests on equipment and structural systems (3, 


pulses with random amplitude to represent | the effects of on dynamic. 
systems (8). ‘Scruton and Harding used a crude explosive charge to excite a 
chimney in order to determine its damping characteristics (20). 
oe simple mechanical pulse-generating device of fairly recent development 
is capable of producing short duration forces of large magnitudes over a wide 
' frequency range that can be controlled to satisfy multimode system response 
(16). Such force pulse generators have been successfully used to simulate the 
| in- -place > motions of up to 500 Hz in equipment weighing up up p to 200, 000 Ib pand 
to measure system impedance functions (17,22). aera? ~ 
_ The procedure for developing a suitable pulse train to simulate the effects — 
of arbitrary dynamic environments on a structure is shown in Fig. 1. A recent 
= demonstrated the ability of force pulse excitation to simulate structural 
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5 acting on structural systems has been investigated both analytically and al 
i 


vibrations produced by continuous forcing functions: as those occurring 


_ during earthquakes (11). During the course of this study, an optimization program f 
was developed to obtain the pulse arameters that yielded the best possible _ 
The computationally efficient optimization scheme of Ref. 11 has been shown a 
~ to yield good results when applied to a variety of multidegree-c -of-freedom (MDOF) — 
systems. However, this computational efficiency was achieved at the expense — 
; s the pulse characteristics: all of the pulse trains at different locations within 
‘the structure were assumed to be scalar multiples of a nominal pulse train. 
This constraint restricted the individual pulses within trains at different locations 
from : selecting the optimum pulse initiation time as well as pulse duration for — 
‘While the limitations previously mentioned may not be very restrictive for _ 
x low-order systems, their effects become more pronounced for high-order dynamic 
: _ systems with large spatial distribution (e.g., tall multistory buildings). In addition, 
as the size of the structure to be tested increases, the energy requirements — 
become more demanding. The lessening of the constraints on the pulse trains 
SYSTEM 


will permit more efficient utilination of the energy content of ‘ins excitation, 
9 thereby leading to a more optimum selection of pulse train parameters. me i 
_ The present study is concerned with the development of an efficient optimization a 
desitthins for selecting the characteristics (initiation time, amplitude, duration) / 
| the pulse trains acting on a distributed system whose structural dynamic 
response to some arbitrary « continuous criterion excitation is is to be simulated 
by the pulse trains. The structure does not necessarily have to be linear, and AS 
the pulse trains may be spatially distributed inthe structure. 
_ The following section presents the optimization algorithm. The validity of = 
this procedure is demonstrated in the section headed ‘‘Applications,”’ by applying | 
it to two systems: (1) A hypothetical 25 DOF system that is representative 
of typical modern tall builc.ngs; and (2) an actual | three- -story f frame structure — 
that has been subjected to. exhaustive analytical and atic studies at 
tte University of California, Berkeley (UCB). ; 


; oe FIG. 1.—Optimum Response Simulation of Multidegree Systems by Pulse Excitation { 
+ 
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_ Statement of Problem. —Consider a linear MDOF system (continuous or with 

lumped parameters) subject to several pulse trains applied _ at locations j = edie 

+, , NFRLOC. Assume that the response of the system, X,(t), at locations, 2 
i , NRESP is desired (see Fig. 2). Let the impulsive ‘Tesponse at 
location. i i due to an impulse at location, be denoted by AY (ty. 


id 


PULSE LOCATION INDEX 
4,04 = COUNTER OF PULSES AT 


Ss Due to ‘the linearity | of the system, the total arene at location i can be. 
found the of different pulse trains from different 


(t)= Excitation at location jand time 
Assume now that desired response at some points of eon te in the structure 
specified and and denoted ted by X(t), i = = 1, 2, The “‘goodness of 


| 
1 
} 
| 
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in which the weighting factors that can be be adjusted to 
or de-emphasize the significance of the fit at different points in a structure, 
since a good fit at some points may be much more important for simulating ~ 

certain phenomena (e.g., , damage and malfunctions). 

: ‘The specified, or criterion, response X(t) are those recorded in the structure _ 

: (as during an earthquake) or obtained from applying known excitation forces 
to the model of the system. Since the objective is to find a pulse excitation — 
that produces a response X,(t) as close as possible to X(t), each — 
of F(t)is restricted tothe form 


of F(t) is res 


in which A = of kth pulse applied at location /; = initiation 


. time of kth pulse applied at location j; w{’ = width of kth pulse applied : at 
location j; and u = = unit step f function; u(a)= ifasO 


Once the linear system fonctions have been determined and the « criteria 
on response established, the problem reduces to finding the set of numbers {t’’, 
AY, wi}, k = 1, 2, ..., N for each of the NFRLOC pulse application vod 
as In order to get some idea about the magnitude of the solution task, consider — 

* typical case where three pulse trains (NFRLOC = 3) are applied to a ‘structure a 
_ whose criteria motion is to be matched over a time span of about N= 20 a, 
/ funcamental periods, 7,, of the structure. Suppose it is decided to use about — 

- one pulse per system time » T, . Thus, at each pulse location, 20 individual — 
bs pulses need to be determined. Since each pulse is characterized by a triplet 


of numbers (1, AY’, wi}, the total number of parameters to be determined e 


(NFRLOO » x > x Triplet = 3 20 x 3 = 180 parameters . 


this optimization procedure i is 


Although a a large number of optimization procedures are available (7), none 


of them can realistically tackle a problem of such magnitude. The next section ; 
fc the procedure for formulating the identifiction problem in a way tha that 
is feasible for implementation with efficient optimization techniques. a 
P Development of Optimization Algorithm.—The basic idea behind the present 
optimization procedure is shown in Fig. 3. In Fig. (a) of this figure, the solution — 
of pulse train parameters for a system with three pulse application locations 
and 20 pulses each in a train is posed as a single optimization problem involving y 
_ the simultaneous optimization of 180 distinct parameters that cl characterize c 
pulse trains. Fig. 3(b) shows how the size of the problem can be ‘reduced to ot 
a manageable level by sequentially ooneny 20 optimization problems, each with 
‘By using the procedure depicted in Fig. 3, an optimization problem of size 
(NFRLOC) x (N) x (Triplet) c can be handled as N sequential problems, each 
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involving (1 (Triplet x NFRLOC) p parameters. This drastic reduction in the order 
_ Of the problem (from an impractical size to a large size) allows the application 

of modern to select the “best” train parameters. 


or 


o4 
T (t/Ty) ao 
ENT SEGMENT to | 


(b) Sequential optimization: 20 operations with 


FIG. 3. —Segmentation of Optimization Ported for System with 3 Pulse Locations — 


location j can be expressed as now that the input atjis 
into two time periods: from tor= TSTART, from 
= TSTART tot = = TEND | as shown in Fig. 4. stalemaniadal aS 
ne excitation at location , j during the | period, 


4 
Simul taneous optimization: 1 operation with 
.... 
: 


TSTART, has been determined and will not be ‘altered. ae ont number of these _ 
7 pulses that have already been determined is NPOLD”’. The problem now is 


the period, t = TSTART to t = TEND S TMAX, in which TSTART = starting : 
time to the present optimization time segment; (0 = TSTART = TMAX); TEND a 
= ending time of the present optimization time segment; (TSTART = = TEND 

= TMAX); and TMAX = upper cutoff limit of optimization time. 


the above definitions, the total response can be expressed as 


to optimally select a total of NPNEW”’ pulses to be applied at location j = 


= XOLD‘(t) + XNEW%(t) 


which | = A? (1 
which XNEW/(t) = 


< 


Stas explained in Fig. 5; T,= 


4 masts hal - 
FROM PREVIOUS SEGRENTS 


t< 


FIG. n of of Pulse Train FIG. 5.—Evaluation of Respon 
ul 


io 

at Location XNEW”’ Due to Single Puls ing 
APOLD\” = amplitude o of pulse num number, k, ‘applied at location j during period, — 
0 = t = TSTART; DPOLDY’ = duration of pulse number k, applied at location — 
Jj during period, 0 = t = TSTART; TPOLD¥’ = initiation time of pulse number, 

_k, applied; at location j during period, 0 = t S| TSTART; APNEW = amplitude — 
of pulse number, n, applied at location j during period, TSTART <t = TEND; 
DPNEW'” = duration of pulse number, n, applied at location j during period, 
TSTART < ¢ S TEND; and TPNEW\’ = initiation time of pulse number n, 
applied at location j during period, TSTART < = TEND. | 

Once X‘/(t), the contribution of the jth pulse train to the response at location — 


| _ i, is computed from Eqs. 7 and 8, the total response at i can be found successively 
b 


+ 


7 


@ 
ng period TSTART J 
| 
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DSPHAT a Nn = 


ad 


in which At, = digitization rate; NTMAX = maximum number of digitized 
in the criteria record; and M MAXL = maximum digitized ‘points 


) START 
"SEGMEN 


«FIG. 6. 6.— Digitized and and Impulse Fun: Function ton ay 


FIG. —Evaluation of Cost Function 


_ rate of At,, which is small compared to the shortest system characteristic time — 
NT) 1) Ar,], NT=1,2,...,. NTMAX ........ (10) 
— 
10-12 can be used ion integrals appearing in an 
4 
toab 
an. 
q 
| 


ny of Cost Function. —Consider the portion of ‘the response shown 
Fig. 7. Let n = index of Gghined data points, n 


_ = total number of data points in present optimization segment. = 
The cost function, J, is set equal to the root mean square deviation error 
a between the criterion response, X,, and the pulse response, X,, summed and 7 
weighted over all response locations and for all data points starting — NTSTRT | 


___ In order to eliminate the possibility ‘that the sequencing of pulses does violate 
_ the physical constraints, the components of the parameter wer. F must be ; 


constrained to ‘satisfy th the following conditions: 
‘TPNEW NP) = = TSTART ‘(for all all 
TPNEW ( (J, NP). + DPNEW TEND (for all NP) 


_[TPNEW (J, NP 1) + DPNEW (J, NP 
[TPNEWU, NP) - 

An APNEW NP) < pag 


pulses at a given location; and DTPMIN = the minimum Lalitha spacing | 
of consecutive pulses applied at any location within a sn optimization 
_ Random Search Algorithm.—Random search optimization algorithms have been 
«a explored and utilized for more than 20 yr. Brooks described a random > 7 
search optimization method in 1958 (4). Analog computer implementations of 
random optimization schemes were published as early as 1959 (6,13). Since _ 
_ ome variety of implementations and their convergence properties have been 


explored. Rastrigin has recently published a bibliography of 760 


in which @ = (a,, a, ... ay)” is a vector of unknown parameters, proceeds 


as follows: (1) An initial parameter value a° is chosen and J(a°) is evaluated; 


(2) trial points a’ « ,, in which , = the given permissible region in the 


_M-dimensional parameter space, are selected from an appropriate probability | 
density function defined over (3) a successful point is one for which 
a Rather than using “pure random search,’’ most algorithms are based on a 
“random creep”’ procedure in which exploratory steps are confined to a 
_hypersphere centered about the latest successful point a’. However, convergence 


by may be extremely slow, since no allowance is made for 


g basic algorithm for 1 minimization of a non- negative ‘cost cates J(o), 
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variations in the nature of the cost function surface as progresses 


The algorithm (12) employed i in this study i is another approach to the determina- 
tion of the locally optimal : step size. Rather | than using a fixed-length step, 
- steps used are random in both length and direction. The adaptation is based . 
7 on the selection of the optimal variance of the step-size distribution as the r 
search progresses. Large variances o” are desirable in the early, exploratory 
a portions of a search. However, in the vicinity of an optimum, a smaller value 
of a” will decrease the probability of overshoot. rns 
UA second and equally i important aspect of the algorithm is that the periodic — a 
step size variance, e, also locates new areas of bacene aad) and — the 


SELECT INITIAL VALUE OF 
4 PARAMETERS AND VARIANCE = 


LOCALLY OPTIMUM VARIANCE 


= N EVALUATIONS | 
THE GLOBAL RANDOM SEARCH 


~ tends toward the global optimum. The overall flow chart ‘of the adaptive random 


“search algorithm used ii in this is given in Fig. 8. 


P LICATIONS 
& _ The system shown in Fig. 9 is a 25- -story office building. designed in accordance 
with applicable building code provisions for recommended lateral force require- 
“ments (2). Modal analysis of this building, treated as a linear elastic strecture, 
_ yields the mode shapes and natural frequencies shown in Fig. 10. 
To: show the simulation procedure | outlined in the preceding, it was decided 
=. use four pulse-excitation locations at modal antinodes (Floors 8, ‘13, 7 
and 23) and to attempt to match the response of two locations (Floors 13 and — 
23). Due to the linearity of the system, its transient response to pulse trains 
could be determined by using the convolution integral approach. The necessary 


‘displacement functions were determined are shown 


[| 7 
| 
AG. 8.—Overall Flow Chart of Adaptive Random Search Algorithm 
: 
a 


(a) Building model 


9.—Plan and Elevation cample FIG. | ilding Characteristics 
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i in which h(t) denotes t ‘the response of , due to 


El Centro 1940 motion ‘was wsed as 


and shock spectra shown in Fig. 12 as solid lines. The earthquake criteria response 
“% was simulated by four suitable pulse trains using the optimization algorithm 4 


4 outlined above. The simulated response is superposed (dotted lines) on top © 
criteria response in » and the time histories of the four required 


2 


0 
13th floor displacement 


(a) 23rd 


time- 


FIG. of and Simulated of 25 System at 


23rd and 13th Floors to El Centro — 
- It is clear from the comparison shown in Fig. 12 thet a good oe is - obtained 
tomes the criterion and simulated response. Note that this simulation of the 
_ motion over a period of =20 sec ¢ required ~ =13 pulses i in each of the four ‘ne 
Referring to Fig. 13 and noting that pulse trains F(t) and F 2(t) correspond — 


baw 
‘e the simulation forces acting on response Locations 1 and 2 (story Number if 


= and 13), it is seen that more simulation force is needed at m,, than at 


In the present example, the error weighting functions, e,, appearing in 


‘Eq. 4 were all equal; consequently the optimization tends 


pF an equally good fit for the response of both locations. 
ame the nature of Eq. 4, only the displacement response is being matched, 


unit 
meters 

eee 


TABLE of Pulse Train Requirements to Excite 25- 25- Story Office Building 


“Number of Average, Maximum, 
pulses seconds in seconds force né 

0.08 
0.35 
_ Summary and y \ | 06 | 0.10 
average 


13. —Optimum Pulse Trains for 25 DOF Structure 


f one is interested in matching other spepense measures, the € expression | for 
the cost function can be modified to become 
— 


A 


“PULSE. SIMULATED RESPONS 
El Centro 1940 Earthquake (Pulsers Located at 23rd, 18th, 1 3th, = say: nee 


zz 


Maximum, ‘Minimum, in Average, in ‘Total, 

Pound force force pounds per pounds per second 

577,000 | 100, 1,788,900 | 
1,000,000 100,000 180,490 23463000 


1,000,000 204,930 2,664,050 


2,049,500 

170,000 

(756, 160 N- (9.4.x 10° N-sec) 


symm. at. 
2 TYPICAL FLOOR PLAN 
: G BRACKET | ore. all | ig 
| 


FOOTING || 
WLOx49 


FRONT ELEVATION SIDE ELEVATION 


ie 
in which 5, = wei hting factors associated with the dis lacement, 
in che; €25 3, enting pia 


velocity, and acceleration ‘response. This approach normally implies a larger 

> 
number of pulses of shorter duration. Thus, the 13 pulses used in this application 
_ for displacement matching could change to 25 pulses with a significant reduction 
im average pulse force required. 
- Pulse trains to induce criteria motions are not unique, as indicated previously. 
Constraints for number of pulses, number of pulse units, maximum and aliintn 


force, pulse duration, initiation time, and location om be stipulated a priori 


a 


a | 
4 "1,000,000 =| 15,000 
4 
q 
ar 
a 
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to accommodate pulse generator performance capabilities. An algorithm to 


‘ optimize for the foregoing variables (constraints) i is currently under "development. 
No effort was made to optimize the pulse train requirements as given in Fig. * 
13, which is summarized in Table |. Force requirements given in the table 
indicate 6 pulse generators will be required for the 23rd floor and 10 units 
for the 18th, 13th, and 8th floors for each direction, for a total of 72 units. 
Performance studies indicate that pulse units _ of 100,000 Ib (444,800 N) force 
capacity provide ; adequate wee--y A for mobility in field tests and particularly 


TABLE 2. —Comparison of Three | Natural 


Results of Present Study 
with UCB Frame and SAP6 = Results of Ref. 21 (Table 4 4, p. p. 99) x rire 


Frequency, Period, Calculated Measured 
inhertz | in seconds frequency in hertz” frequency in hertz 
0.4387 (2.295 
0.1258 8.345 
0.0633 


FIG. 15.—Mode | Shapes of UCB Frame Determined by SAP6 


UCB —A three-story steel frame structure consisting of | two 

single- bay, moment-resistant frames has been extensively tested by means of 
the earthquake simulation facility at the University of California, Berkeley (UCB). 


Detailed and analytical studies ‘Fegarding the subject frame 


In the present study, a eaten purpose finite. element computer program 
-SAP6 (18) was used to analyze a linear model (henceforth referred to as the 


| 
Mode 
— 4 
“4 1 
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of the mathematical model of the UCB frame were re developed from information — 

~ given in Ref. 21. The frequencies and mode shapes of the UCB frame and 

_ the corresponding measured and calculated frequencies given in Ref. aaa al 

Pulse trains were to be applied at ‘thet three floor locations. The needed impulsive q 

response functions were analytically determined and are shown in Fig. 16. The o 

criteria response to the El Centro 1940 earthquake was , likewise analytically | 


determined by using SAP6, and the results are shown as solid lines in response 


= pel 


Ob 
FIG. 17.—Comparison of Criterion and FIG. -18.—Optimum Pulse Trains: for 
Simulated Response of UCB Frame UCBFrame OFT 
time history of Fig. 17. The adaptive random research optimization procedure 
was again used to determine the required pulse trains. The resulting simulated — 
motion and the three required pulse t trains for each floor location shown | 


' A summary of the pulse generator performance requirements is shown in 


a -4 7.5x10> 7.5 x 10 = 
[ 
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_ TABLE 3.—Summary of Pulse Train Requirements to Excite Three-Story UCB Test 


Average, m, | Average, 
in seconds | ins in seconds | poun r = 


SSS 
Pa 


a BLE 4. —semmary os of Pulse Train Requirem > Excite Thee e-Story U UCB tory UCB Test 


me s in nds | in seconds 


(8.4 x 10° N) 
4,076 
(18 x 10° N) 
6,387 


(28. 4x 10° N) 
- and are mounted in opposite directions, for a total of six pulse generators. — 
‘Units located on the first floor are required | to produce thrusts up to 25,000 


| -Ibf (111,200 N) with pulse durations as short as 9 msec. oo | 


Alternate studies were conducted using pulse generators on only one floor 
"(either the first, second, or third floor) and also on any two floors at a time. 
The latter study, of course, leads to six combinations. Restrictions were placed | 
on the pulse 1 units not to exceed 10, 000 lbf thrust. The results of the study 


in Table 4. The differences resulting from. the ‘tables a are e substantial. The third 
floor requires a maximum pulse of 4,322 lbf (19,200 N) and a minimum pulse _ 
duration of 90 msec for units placed on the third floor. When the optniention 


algorithm is complete for pulser locations and impulse minimization, | further 
__ The programmable gas “pulse generator shown in Fig. 19 is currently under a 
construction. Initial use provides for cold gas but the system is adaptable for 
both steam and chemically generated hot gas. Thrust amplitudes are controlled 
in the on-state by positioning the metering plug for flow control. Off-state for — 
- the pulse occurs by signaling the hydraulic actuator to move the metering plug 
to seal off gas flow at the nozzle. This device will be “used for earthquake 7 
and antiearthquake investigations on the UCB structure shown in Fig. 14. In 
_ these tests, response motions of the structure will be compared when: (1) The 


Floor | 
3 | 0.046 03 0.100 | 0.039 | 2,430 
Average and 75 | 0.048 0.100 0.009 4,507 
— 
| 
seconds 
g 7 
q 
| 


-PULSE-SIMULATED RESPONSE 


poun pound- pounds per second pounds per second 


‘Li 


09 


Total, in 


pounds per second 


a9 x 10° om No (B41 N- x 10° N-sec) pal 
(45 x 10° N) 0348 Ne sec) x 10° N-sec) 
984.N-sec) | 6 x 10° N-scc) _ 


PROGRAM 
PARALLEL VALVE SYSTEM. 
“ico 


SERVO SLAVE VALVE 
SERVO PILOT VALVE oe 


_ FIG. 19.—Programmable Gas Pulse Generator Cold Gas System _ 


structure is vibrated by the shake table as base motion; and (2) the structure 
is vibrated with pulsers attached to the frame and \ with the shaker table blocked. - 
Additionally, the structure will be vibrated by table for several 


one 
5. euch we he q 
Structure to El Centro all Three Floors 
bee 25,470 puto din 3,825 = 
@25x10?N-sec) 
q equirements if Units are 
Maximum, in T Minimum, in in 
7 
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‘train characteristics ‘to > be used for simulating the re response for general = ; 
of structural systems to | arbitrary dynamic environments is presented. The 
optimization procedure uses an adaptive random search algorithm tat in- 


corporates a periodic exploratory search for the optimal step- size variance, which ; 


significantly improves the convergence characteristics, 


The method is applied to two example structures: (1) A 25- -story, 
Bey and (2) a moment- resisting frame that has been extensively investigated 
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APPROXIMATE DESIGN BY GEOMETRIC SERIES 


~ By ASCE snd Nicholas Willems,” M. ASCE 


Often” a designer has to | minimize the Structural weight ‘of a given frame 

; configuration for a given loading. The most natural and rigorous way of dealing 
with such problems is to make use of mathematical programming methods. — 

_ It is also well known that strict application of available constrained minimization . 

_ techniques invariably has failed to produce fully satisfactory results and that 
this has led to design Procedures b based on approximation concepts (5). By 

: constructing a high quality approximation of the minimizing process, a practical 
and meaningful design is generally obtained. In practice, the results of a structural 
_ weight optimization normally serve only as a first approximation of the final 
design because of various practical design and construction considerations. 
‘Therefore, it is important for a practicing engineer to have access toa sufficiently ~ 
accurate optimal design method which uses an algorithm and | can be ‘executed 

fairly simply on a desk calculator that c can be. programmed. It is s further desirable 
that the input data be basic, i.e., contain only geometrical and load data. As 

input, neither the starting values of member cross sections nor their mutual — 

rf ratios should be required. Such an algorithm will be much easier, quicker, and : 
cheaper to use than more exact methods which make use - of mathematical i 

4, and 8). The weight of each frame member depends upon 


has to represent an automatic successive reduction of characteristic moments, 
without violation of equilibrium, resulting in the successive creation of plastic 

oi 


hinges and a corresponding collapse mechanism. 

Itw was 1s previously shown (6) that the plastic moment can n be ‘used 

¥ for automatic optimization by using an iteration procedure. Later (7), it was 

also shown that correction and carry- over moments can be made to form terms 

of infinite convergent geometrical series with known exact sums. If this is the 

ease, iteration loops can be eliminated and substituted by closed-form formulas 
Representing the sums of an infinite number of iterations. This drastically reduces [| 

the required size of core memory and the operational time. This way it | _ 
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5 a maximum bending moment acting on that member. The algorithm, therefore, 7 


| 


~ possible to sabes: an optimization problem of an unsymmetrically loaded four-story, 
four-bay rectangular plane frame using a minicomputer with only a 4K direct © 
- read- write memory. Use is also made of partitioning of multistory frames (2). ~ 
7 Starting with the top story and moving downward, each story is optimized — 
— This procedure is equivalent to the assumption that in each column © 
_-. every story a plastic hinge will develop just above the beams of the — 
below it. Therefore, this assumption is also the first condition of reaching a 
global optimum. Often, this assumption proves to be correct, and the optimization | 
procedure, in that case, , could re result in a global optimum. The second condition — 
in reaching a global optimum re requires that the true . collapse mechanism of the © 
= frame be the one obtained by combination of beam and sway mechanisms. 
‘This requires two successive hinges in each beam, two hinges in the last column, 
and one hinge at the bottom of all: the other columns of that story. If any 


approximate value, but with small error. Comparative analysis has showin such 


an error to be less than 4%, The | explanation why this error is “so small is 


_ For the sake of simplicity, the method is restricted to plane rectangular frames — 
composed of prismatic members subjected to static loads. Vertical loads are 
always applied at the beam midpoints, and horizontal loads, which are assumed 
to act from left to right, are applied at the floor levels. 
_ Use is also made of the simplified plastic theory with its usual approximations — an 
limitations, i.e., work- hardening of steel, no occurrence of frame 
instability, and small deformations at collapse. 
_ A further limitation is that the first admissible moment configuration must 
- be such that bending moments at the ‘beam centers due to gravity or combined 
loading a are larger than the ‘moments ‘at both beam: ends. For normal icine 


conditions under or existing building codes and frame configurations, this is always 


linearized objective function, the weight o of the. 
_by the sum of the products of the length and the required full plastic moment : 

= capacity for each frame member. Therefore, reducing plastic moments obviously 

__ will reduce the overall weight. At collapse, a true plastic moment diagram must 

. : satisfy equilibrium and have a sufficient number of plastic hinges for the formation — 
ofa collapse mechanism. If all three characteristic beam moments are successively 
corrected so as to equalize two out of three beam moments, i.e., create two 
adjacent beam hinges, and if the largest of the three beam moments is constantly 
being reduced by this equalization, the total weight of the frame will also vedi 
decreasing. If this procedure is repeated an infinite number of times, the possibility 

: exists of reaching a global optimum | provided the the two conditions, mentioned - 

os before, are satisfied. If this is not the case, however, a local one will be reached - 

q but of sufficient accuracy. The frame weight will be reduced by iteration cycle, 
and a collapse mechanism will also be created. 
_ For each frame story, the first maximum beam moment reduction and the 

“creation of two beam hinges is performed in the first bay, | J = 1, when . 


‘suitable correction moment, is first applied ‘to the ‘Tight-hand 


i 
| 
| 
» 
| 
q 
| 
| 
i | 


“end- -moment and then, for the sake of equilibrium, ca carried over with a a carry-over oll 


>. ~ factor (CO) of —1/2 to the midbeam moment. Once two moments are equalized, 


= will remain equal because, as will be seen, they will be changed in each 
~ future cycle by exactly the same amounts. To satisfy equilibrium of the common 
a beam joint between the first and second bay, this first correction, Da, 1), has 
~ to be substracted from the left beam end-moment of the second bay. Now, 
_ to satisfy beam equilibrium of the second bay, a part of this correction, —D(1, 1), 
4 will be carried over to midspan, with CO + 1/2, and the remainder to the 
Tight- hand beam end-moment (with a CO = 1). If two-thirds of —D(1, 1) are 
- carried to the midspan and only one-third to the right-hand beam end, both 
_ moments at the midspan and the right end will be reduced by the same amount 
of —1/3D(1,1). This correction from the second bay will effect the third bay 
- in the same manner and will be carried over to the midspan and right- hand © 
: beam end, after another reduction by 1/3. This value, so reduced, will now 
_ pass on to the next bay, and so on. The first correction, D(1, 1), and its reduced 
values will thus flow along all beams from left to right, through the last column 
from top to bottom as one and the same correction, back to the first column, > 
3 and finally to the left end beam moment of the first bay as a correction of 
— (-1)” x (DAL, 1))/3"~ ‘o. in which n = number of bays. This is the end of the | 
first cycle of equalization of moments in the first bay. If this same process 
is now repeated, an infinite number of times, all the corrections at any one > 
and the same section will form an infinite decreasing geometric series due to 
the equalization of two ‘moments in the first bay. The sum of all the terms: 
of this series can easily be obtained. Next, two beam moments, at ‘midspan 
_ and right-hand beam end, in the second bay are equalized and corrections flown 
around the frame an infinite number of times. Their sums for each section 
are obtained. This same process is now repeated for all n bays of a particular 
_ story, resulting in nm sums of n infinite decreasing series for each characteristic _ 
section. As will be shown later, these sums at at any one particular section luckily » 
form elements of a matrix that can be generated easily. Therefore, there is 
no need to generate infinite series and find their sums. All that is needed is 
i. generate two such matrices—one for midspan beam moments and another 
for right-hand beam end-moments. Each row j of the matrix elements represents 


the correction sums for the corresponding bay j. 


The whole method can be presented by the general flow chart (Fig. 1). This 
be 
method will be first applied to a single story ‘multibay frame and then to a __ 
_ multistory frame. As an illustration of the method a completely worked ‘* 
example is also given. are 
_ Moments are taken as positive if they act clockwise at member ends and | 
a. t beam midspan if they are sagging moments. Bending mo moment diagrams are 
drawn on the tension side of each member. 
Single Story.—Let us first consider a multibay single- story frame (Fig. 2). ad 
The bays are numbered j = 1, 2 ... nm. Each bay has seven potential critical 
“sections where plastic hinges might develop. moments at at these locations 
are called M(j, 1), M(j,2) ... M(j,7). | SU 


ering ‘the ‘simple beam and sway ‘(panel collapse) mechanisms an ‘initial set of 
‘bending moments is obtained. The necessary correction moments to establish | 


joint equilibrium are carried back from corrected joints to all midbeam moments — 


| 
| 
| 
| 
| 
4 
| 
\ 
q | 


Read-in frame and | 


Create first admissible cont 
a 


Generate matrix elements, | | 
i.e., the sums of 


Write final moments 4 


Find all beam 


. 1.—General Chart 


except in the first bay, where both end ‘moments are are corrected ~ the ‘same 


emoumt. Therefore, the corrected first admissible moments are it 


sw moment, 


ibs 
La 
3 


ity 1) x eid ré 
4 
in which j = 2,3. nat ‘midspan. & 


which j = 1, 2... at right-hand beam end. 

3) = (4) 


in which j = 2, 3 ... n at left-hand beam end. Instead of assigning equal portions 


of the total sway moment to each column, it is also possible to assign — Wl) 


ore —To start the optimization procedure, the beam moments in 


the first bay (Uj = = +d at locations S and 5 (Fig. 2) a are equalized by 3 applying 


a correction moment: 


This correction, D(1,1), is added to the right-hand end beam moment, = 5), 
and, for beam equilibrium, carried back to the midbeam moment, M (1,4), with — 

d a —1/2 carry-over factor, i.e., two plastic hinges in the first bay are formed — 
and the midspan moment is reduced. Actually, this means that half of this 
correction D(1,1) is substracted from M(1,4) and the whole amount added to 
- M(1,5). The disturbance of equilibrium of the joint at the right- rhand end of | 
the first beam is therefore J D(1,1). To ) restore ¢ equilibrium of that joint, -D(l, 1) 
is added to the beam moment in the next bay, M(2, 3). Next, to restore beam fa 4 

- equilibrium in the second bay, —2/3 D(1,1) is carried over (with a carry- -over 

7 ~ factor CO = +1/2 to M(2,4) and —1/3 D(1,1) (with CO = +1) to M(2,5). b. 
_ This means that the midbeam and right- -hand end-beam moments of the second - 
_ bay are reduced by the same amount of 1/3 D(, dD. To restore joint equilibrium © 
now at tl the. right-hand end of the beam of the second bay, +1/3 DA, ‘is 

added to M(3,3) of the next bay. To establish equilibrium of the beam, here 


4 
i 
| * |~ Mil, 4) 4) MU,4) pdr | 


.* the third bay, this ounainn. 1/3 1/3 DA, 1) is again divided into two-thirds 7 


and one-third, as was done previously in the second bay. . Half of it, ie., [1/(- -3)? r. 
D(I,1) is added to M(3,4) and the same amount to M(3, 5). 
i, This process is again repeated for the next bay, and so on, , until the last 
bs bay is reached with the correction [D(1, 1)/(-— 3)"""] added to M(n, 4) and M(n, 5). a 
When the equilibrium of the last joint at the right-hand beam end of the last, — 
nth, bay is restored by adding to M(n,6) (—1)”- [D(1, 1)/(3)""'] to the top 
of last column, the sway equilibrium of the whole frame has been violated 
and needs to be restored. The sway equilibrium condition requires that the 
sum of all corrections added to any of the column moments be zero. To achieve 
= and also to achieve a “‘minimum column design,’’ the same correction 
. applied at section (n,6) is added to (n,7), and with opposite sign to sections ; 
(1, 1) and (1, 2) of the first column. These corrections affect, therefore, all beams | 7 
~ from left to right, then the last column, and finally the first one. In this way, 


the original minimum sway moments ‘remain intact in “all columns except the 


first and the last ones. This represents ** minimum stanchion design.”” 


_ This completes the first iteration cycle by arriving at the top section of the 
(fies column. The second cycle starts with reestablishing equilibrium of — a 
joint between first column and beam by to M(1,3) the « correction with 


a the first bay, one-third of that correction, i. 1)": [D(l, 
(—1/3)" 1), is added to both moments M(1, 4) and M(1, 5), already equalized 
in the first cycle of iteration. Therefore, they remain equal, i.e., act as hinges. — 

The successive restoring of the joint equilibrium at the right-hand beam sides" 
as well as the sway equilibrium in the last and first columns is again. carried 

out by a procedure identical to the one used in the first cycle. Each next 

cycle is performed exactly in the | same manner. This is done an infinite number 
of times. In this way, due to the equalization of two moments in the first 
bay, all the corrections at any frame section form terms of an identical decreasing 
geometric series, which has a known sum for an infinite number of terms. 

For or example, for a two-bay frame (n= = = 2) 

3)= cil) +— DU, 

and M (2,5) with a new correction factor of the same form as one, 3.7 The The correction 

factor for the second id bay is calcul is calculated 2.2 

1,2) = — [M(2, 4) 

MQ." 4) “wen, 5) are the first admissible moments, ¢ given in Eqs. 2 and 3 

due to the fact that both moments were previously changed by the same amount 
during the infinite correction cycling process with D(1, 1). The same operation © 
is Tepeated as } before with D(A, 1). This then for all until the 


“| 
. 

| 

— 
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q last is also applied by an infinite number of times 


around the frame as before. ‘There are now n infinite geometric series at each 
_ characteristic frame section. . Knowing the sums of all n correction series at 
each critical section, it is possible to write expressions for the final moments ie 
at every critical section in closed form. In case of beam moments, MUj,4) — 
and M(j,5), these sums are grouped as elements. (j,k) of a (nm X n) matrix 
[A] and [B], respectively. For any bay j, instead of sums the SS 
whole row j of matrix elements is used. Any element in that particular row, 
ja represents t the sum of the infinite ‘number of cycling corrections started 


Therefore, ‘like j, k = 
write final moments apply for an odd and an even total number of bays, mos 
Cc 


GAL gin: 
olumn Moments mA, 1) = MA, 2) = M(1,1) 
Column Moments MG, 6) = MU, N= M(j, 6) = 
in which j = 1, 2, 


‘Column man Moments (1 6) = M(n, 7) = M(n, 6) + AM(n, 6) 


‘Beam Moments M(j,3) = - (see bom 


Beam Moments at d-span MG, 4) = M(j,4) +4 M(j,4) 


in which j goes fr om ton 15: 
Beam M Moments 4) = M(j,4) Jj + ky 4) 


in which ae coefficients a,, are elements of a square m X m1 matrix. ‘This 


coefficient matrix A has the property 


Therefore, it is sufficient to form the first row and first column of this matrix. 


The first row elements are 


ve 


in which j = 2, mentioned cartier, ave have the” 


= 


F 
with egualizauion of two beam moments in the bav corresponding to column 
10) 
| 
12) 
(13) 


DA, +300. 5) — 3? D(1,6) + 3° 7)] 


coefficients are of another square matrix with the 

same property expressed in Eq. 15. The difference between these two matrices 

| in the elements of the principal diagonal. All the Tot on the i 


= =3" 


in n which j jandk = 1,2, ... n. All ‘the other e elements of ‘this matrix are —— 
the same as in the matrix for M(/,4), described previously. — . 
= even, the final moments 1)... 7) 


= 


Columa Moments 6) = 
= (n- I)as sive 
Column M(n, 6) = = M(n, = -M(n + 6) = 
it is easily that the only difference between Eqs. -23 and the 
corresponding Eqs. 8 and 10 is that each term under le summation be 4 has 
an opposite sign, and the denominator is changed from 3" - 
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j and the columns the subscript k. For example, for a seven-bay frame (n_ 
7) the coefficient matrix for M(j,4)is ee hepois i 
| 
q which rows j and columns k go from 1-7. For each bay j, the corresponding | é 
row j has to ements between 
— 
7 
4 


4 *an 


Beam Moments MA, = 2) in Eq. | 18 
= MU, 3) = 


The coefficients a, Gy» , are elements of a matrix similar to 14] for n = odd except 
for opposite signs. Therefore, the first row elements ares” 

in which k = 1,2... n. - The first column elements are ou sd nisl 


y-! 
=(- x 


The b,, are of a matrix similar to he previous matrix Bi. in 


on was an odd number. Eq. 20 does not apply for the elements are the same 
ast those ‘ove in = 26 and 27 for matrix ee 


ibe Story i+] _ 


ea 


dite to the top story with h respect to first ads admissible t moments, re 
Any Other Story (i = 2, 3 ... .m) Yo 

First Admissible Moments. Again, in the formation of the 

_ Sway Moments.—The sway moments of any ith story are pare from the . 

:: sway moments of the story just above it, story i — 1, (Fig. 3) plus an additional 


wane fo: apy 
animes 
= M(j, 5) + AM(j, 5) = MU, 5) + — (28) 
4 
| | 
Move 
» — Whatever wa aid far for ingle_cto frame in 


= 2) = MU, 6) = MU.7) = - 


Beam corrections to the beam moments due to the restoration 


of joint equilibrium now are larger than those in the top story because the | 

column moments of the story above il and of the ith story itself need to a 
? be compensated for in the second step ) of the procedure to obtain first admissible _ 

M(1,3) = 

in which j = 2, 3 4 Fare: > 

in which j = 2, 3 ... —_— 


MU, 5) = C(i- 


in which j = 1,2. 
_ Optimization. —Before the optimization of story i can start, the joint equilibrium . 
at the top of the first and the last column of that story must be reestablished. 
The moments M(i,1) and M(n,7) of the (i - 1) story have undergone during 
ee — 1)th story optimization the corrections AM(1, 1) (Eq. 8 or 21) and AM(n,7) 
(Eq. 10 or 23). These corrections are now added to M(1,3) and M(n,5) with © 
opposite sign and carried over to the beam midsection, i.e., to M(1,4) and 
wi After these two corrections the optimization procedure is the same as for 
a single story as explained earlier. Occasionally, it will happen that in the last 
=n) the moment M(n,4) < M(n,5). In that case, the 
therefore this correction is zero. When the last two moment 
at the end of the optimization process are equalized and “equilibrium: at the 
& of the last column is reestablished, to satisfy sway equilibrium the negative 
of the same amount is added to M(n,7). » 
 -Egs. 8-28 are valid for the final moments of the ith story. Only C(1) has 


correction coming from the ith story. Therefore, we have 
) 
q 
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matrices for any previous story are the same for any 
other story provided the number of bays, n, for these two stories is the same. _ . 
if n(i) # n(i — 1), new matrices ‘Must be be generated using the same formulas. ‘4 


my 


4 Loading and Fi Final M (10 


—Example. First Admissible Moments 


up to a four-story, four-bay -unsymmetrically loaded frame and up to 15 bays 7 Z 
Example. —To illustrate the application of the e method : and check its accuracy 
- anexample was taken from Professor John Munro, of Imperial College in London, hy 


England (10). In Fig. 4, the two-story, three-bay oma, its its” 


i. programmable desk calculator with 4K R-W memory. This program can handle ey 
. 
. 
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final bending moments at are re shown. ‘The ‘weight”’ obtained by Munro 
was 12.040. The writers achieved a weight of 12.054, an error ‘of 0.12%. Professor 
Munro first used a so-called ‘mesh unsafe’’ program followed by ‘‘nodal unsafe’’ a 
program. In each case, the kinematic variables were evaluated directly and 
the static variables were obtained as th the simplex multipliers at optimality (Figs. 


1. The method ciilenal in ai paper is an approximate one; however, a 
_ comparison with several problems solved by using exact mathematical program- 
ming showed that the error never exceeded 4%. For an initial design approximation | 


The closed form formula develop ed on the basis of infinite, decreasi 


this is considered to be sufficiently accurate. = 


geometrical series cut reduced computer memory capacity and execution time 

to such an extent that a minicomputer with a 4K memory is sufficient to optimize 

a four-story, four-bay rectangular multistory frame. 
3 Although the described cycling process might look complicated, the use 

_ of equations to determine the final moments is simple. For this reason the = 
development of a program for a any available desk calculator i is easy. 
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the starting sizes of frame elements might use this program is 
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‘Seren i story number; = 
column number in matrices, [A] and [B];_ 
= number of bays in any story i 
= column wind moments in ‘story i; 
i,j) basic correction moment; 
rey at = span length of a bay 
first admissible column moments, in which 
M(j,3), M(j,4),M(j, 5) = first admissible beam moments, in which 


M(j,3), M(j,4),M(j,5) = optimized beam moments; and 
Tor AM = total moment change due to optimization. _ 


re 
however 
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Qu 


The 
| the a 


16501 Ee SEPTEMBER 1981 


EVALUATION oF INELASTIO. 

oF UATION OF F INELASTIC 
SEISMIC DESIGN SPECTRA" 

By Stephen A. Mahin," M. ASCE and Vitelmo Vv. Bertero,? F. ASCE 


‘Because of economic considerations, structures are usually designed so that 
some of the. energy input during. severe earthquake ground motions is dissipated _ 
by means of inelastic deformations. _ These deformations, however, must be — 
- controlled in order to prevent collapse due to exhaustion of the structure’s 
energy dissipation capacity or due to excessive lateral displacements. Moreover, 
_ any economic savings made possible by allowing inelastic response should be 
reconciled with potential costs « of damage to the building’s ‘structural : al system, 
- nonstructural elements, and contents. Proper selection of seismic design forces, 
accounting for the expected intensity of ground shaking and the acceptable 
severity of inelastic response, is a complex task, and current building code 
recommendations (38) offer only qualitative guidelines. A number of other ee 
methods, which explicitly account for these and other parameters, have been 
‘Tecently Proposed. However, because of of the variety of such methods ds available 
and the many ‘simplifying assumptions 1s generally involved in their development, ES 
it is not clear how much confidence a designer can place inthem. _ rieatiniain pet 
The purpose of this paper is to review general methods available for specifying — 
seismic design forces for structures that can be allowed to suffer limited amounts 
of inelastic deformations during severe earthquake ground motions and to examine 
bs ‘some of the problems that exist in their specification. Following a a brief review 
of available procedures, these problems are illustrated by evaluating | the reliability : 
of two representative methods. This evaluation focuses, in turn, on three 
interrelated problem areas: (1) The ability of the methods to limit peak inelastic 
deformations to specified values is examined for a wide variety of ground motion 
records; (2) the sensitivity of computed inelastic response to uncertainties in 
structural damping an and hysteretic: characteristics is investigated; and the 
i. Fenn of the maximum displacement ductility factor, i.e., the ratio of maximum a 
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to yield as the primary design criterion is assessed in view 
*. information obtained on other key response parameters, such as maximum and 
permanent displacements, maximum cyclic inelastic deformations, , numbers of 
_ yield excursions and reversals, and energy dissipation demands. To simplify 
‘ this presentation, emphasis will be Eanes on single degree-of-freedom (SDOF) © 
OF Metnoos For Osrainine Inevastic Desicn Resronse SPECTRA 
Preliminary design forces can be derived from inelastic po response ae 
_(IDRS) which usually plot, as a function of period, the yielding seismic resistance 
coefficient or pseudo-acceleration required for a specified level | Of 
_ response. A number of methods exist for constructing such IDRS. — 
‘Time History Analysis Methods. —Nonlinear time history response 
: provide | the most reliable means for constructing an IDRS for a given ground P 
_ motion. For SDOF systems, the equations of motion may be conveniently — 
- nondimensionalized in terms of the fundamental period, 7; the viscous damping — 
ratio, €&; and a parameter, n, equal to the system’s yield resistance, F,, divided _ 
by the product of its mass, M, and the peak ground acceleration, "A (9, 20, 


‘= 27, and 33). Results of parametric studies based on such formulations a” nm, 


be used to estimate the strength required to limit the maximum displacement _ 
i: ductility demand, or another suitable response parameter, to a specified value 
* for the earthquake record considered (7-9, 18, 19, 22, 28, and 39-41). Statistical 

or probabilistic analysis of these estimates for an ensemble of recorded or artificial 

(13) accelerograms with characteristics appropriate for a site can be used to 

_ develop IDRS corresponding to the specified maximum inelastic re response (8, 

20, 24, 27, and 33). While IDRS constructed in this manner are powerful tools 

_ for preliminary design, the efforts needed to obtain the appropriate ensemble 


for most design projects, 

_ Random Vibration Analysis Methods. —Random vibration theory can be used 
= estimate | inelastic response in probabilistic terms. While an exact solution _ 

for nonlinear systems has not yet been derived, a number of approximate : solutions od 

_ Methods Based on Linear | Elastic Response Spectra.—Because of the whites 
practical difficulties involved with the preceding methods, it is more common — 
- in practice to use approximate methods that modify a smoothed linear elastic 

_ design response spectrum (LEDRS) developed for a site. One approach to this 

problem is to replace the inelastic system by an equivalent linear elastic system 

and use the LEDRS directly (16, 17, 33, and 34). The effective stiffness and ae 

_ viscous damping coefficients for the linearized system can be derived on the 


basis of a steady-state corresponding toa specified 


motions are re still required to determine ‘the ‘effective stiffness and damping. of 
In other methods, the LEDRS is modified to obtain an IDRS using factors is 
that depend on the acceptable level of inelastic response. . The simplest method 


fer doing this is to Teduce- the ned forces obtained i the — by A 


| 

I 
| 


a factor that is veibeihieas of structural | period. The tentative e design n recommen- 

_ dations of the Applied Technology Council (ATC) are representative of this 
: - approach (36) and are used as the basis of some of the following examples. — 
The ATC suggests that, at the point of first significant yielding, a structure 
‘should be able to resist a total lateral seismic force equal to ‘its total a 


“which i is given by Phi vali 

~s 

RT 

need not 2. SA for rock or | firm soil in which 
2 Tepresents | the effective peak acceleration; A , Tepresents the effective peak 

Be. aig -related acceleration, R is an empirical seismic response modification — 

coef ficient intended to account for the damping as well as the inelastic def ormation 

capacity of the structure; and Ss = a seismic coefficient accounting for the 


| qu 


FOR SPECIAL MOMENT- RESISTING 
REINFORCED CONCRETE FRAMES 


R=70 Cy=60 2 


of a Ci times the displacement when first significant yielding occurs, 
: i.e., the nominal displacement ductility factor equals C,. In some of the ng oe 
examples, these variables were specified to correspond to reinforced concrete, 
_ special (ductile) moment-resisting frames located on stiff soil in a region of = 
Severe seismic exposure (A, = 0.4,A, = 0.4, Cc, = 6, R= and S$ = 1). 
In other methods for obtaining IDRS from LEDRS, factors have 
been analytically derived by evaluating the response of elasto-perfectly plastic 5 
_ (EPP) SDOF systems, and are usually expressed as a function of period as — 
well as of the maximum specified displacement ductility, 1 (2, 20, 25-29, 33, _ 
40, and 41). Since these types of reduction factors are based on results obtained © 
for idealized structural models subjected to limited numbers of ground motion — 
= records, designers have usually been cautioned against applying them to structures - 
that exhibit Substantially different hysteretic or damping characteristics or 
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SEPTEMBER 1981 
Based on analytical oat Newmark and Hall (28) concluded that: (1) i 7 
short period structures, any significant reduction in design forces required fo 
elastic response would result in unacceptably large ductilities; (2) for moderate 
” 4 period structuses, the energy absorbed by an inelastic structure at its maximum 


4 7 ; and (3) for relatively long period — 

structures, the maximum displacements of elastic and inelastic systems were : 
nearly equal, so a strength modification factor of » va would be appropriate. 
This well-known method (28) is representative of this general approach for = 4 : 
_ obtaining IDRS from LEDRS, and is used to establish the strength of structures &§ 
considered in some of the subsequent examples. In most of the examples 
considered, a maximum displacement ductility factor of 4 is specified, and the 


LEDRS is developed, (28) assuming a ratio of peak ground velocity, V, to 
acceleration, A, of 48 in./sec/g (1,220 mm/s/g) and 5% viscous damping. — 
The IDRS corresponding to these values is shown in Fig. LOPS eevee 
~ Application to Multiple Degree-of-Freedom Sy stems.—Most IDRS are not strictly 
applicable to © conventional multistory buildings since it is usually very difficult 
7 to idealize such structures as EPP SDOF systems (21, 28, and 30). A number 
of studies of inelastic shear building models have been performed (2, 20, “a “4 
z and 42); however, general agreement on methods for deriving inelastic design 


story shear for ‘such has not yet been achieved. A number (of 
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Although convenient for preliminary design, | such approximate methods have 
been found to have a number of defects (3, 7, and 15). Additional information — . 


_ regarding problems in applying IDRS to multiple degree-of-freedom systems 
is presented in Bertero, et al. (7). nib 
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structural response, is ; desirable to assess the reliability with which 

TABLE 1 —Standard- -Type Ground Motion Records Considere Ps ‘ 


Earthquake Component acceleration” 


Centro) (1,305.3) 


Centro) 

1949, Western 


limit inelastic deformations to specifi 

ductility demands are presented for ideal EPP SDOF systems, “ranging in period : 
0.1 sec-2.0 sec, designed using IDRS suggested by Newmark and Hall 

(28) and by the ATC (36). Comparison of ductility demands obtained for a 

_ variety of ground motions with the specified ductility values gives an ara 


of the reliability of these representative types _of IDRS. these 


= SEISMIC DESIGN SPECTRA 
_ investigators have either proposed or derived approximate seismic design forces § 
_ for multistory building structures based on elastic root-mean-squared modal _ 
superposition techniques using the initial dynamic characteristics of the structure Jf 
_and DR ahtained far an nO avetem (2 A 16 97 32 25 an 2 7 
| 
il 
(610.6) 
38.90 


“the gi ground m motion records to have the same peak 
Newmark-Hall Type IDRS.—-Computed ductility demands corresponding to 
the basic Newmark-Hall IDRS with specified ductilities of 1, 4, and 8 are shown 
; in Fig. 2 for the SOOW component of the 1940, El Centro record. ‘Ductility — 
demands vary ‘considerably with period and can differ significantly from the - 
specified value. For this record, ductility demands are usually conservatively 4 
less than the specified values; the maximum discrepancy increases with increasing 
values of specified ductility so that , for a specified ductility of 8, the actual 
“ductility « demands are as low as 43% of this value. 
TABLE 2.—Average ‘Ductilities (and Coefficients* of Variation given between) for 
= Systems with Periods between 0.1 and 2.0 Seconds, nas to Standard 


= 48 
inches per inches al 
‘ second per second per 
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__ displacement 


Cyclic 

displacement 
Permanent 
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“Coefficient of variation given asa percentage, 


_ Results for Standard Ground Motion Records.—Ductility demands were ee 
computed for 10 standard accelerograms representative of severe ground motions 
recorded on firm ground at moderate epicentral distances (Table 1). The results 
are plotted in Fig. 3 and summarized in Table 2. For a specified ductility of 
4 the average maximum displacement ductilities computed for these records | a 
Fig 3(a)] are generally smaller (29%) than the specified value. The computed 7 
displacement ductilities are very sensitive to the ground motion record as well I 
as to —_ period; values as a, as 8 misnin as low as s 0. as 0.5 were obtain were obtained in some 


} 
= oF. 
| EPP degrading} EPP "model — 
6) | 6) | | 
(103) | (83) (104) 37. | 4751 
co | | cm | | 


_ cases. This large variation is indicated by the coefficients of variation plotted ' 
in Fig. 3(a), which average 54% | over the entire period range (Table 2). For — 
a specified ductility of 8, the same general trends are observed, as” can be % 
seen in Fig. 3(c) and Table 2. In this case, the ductility demands average 22% 
less than the specified value, indicating that the basic Newmark-Hall type IDRS 
7) _ somewhat less conservative for higher ductility values, as has been noted _ 
(20, and 29). Moreover, ductility demands in excess of 27 were computed in 
some cases. Such high ductilities may be difficult to develop in practice. > Anse 
1 (St) oor (AGG) mad 
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to noah PERIOD, “PERIOD, SEC. * PERIOD, SEC 
(a) EPP MODELS, 4 DEGRADING MODELS | (c) EPP MODELS, pe 
7 AIG. 3. —Response of SDOF Systems Design Using Newmark-Hall IDRS with € = 
0. 05 and V/A = 48 in./sec/g (1,219 mm/s/g) Subjected to Standard Records: (a) 


4 be One possible rea reason for the relatively large difference between the ¢ specified - 
and average ductility demands in the intermediate period range is that the LEDRS_" 
used is based on a V/A ratio of 48 in./sec/g (1,220 mm/s/g) as commonly 
_Tecommended for firm soil Sites (28). The average V/A ratio for the standard — 


this Seriod range. "When the correct average V/ A ratiois used to derive 
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the IDRS the ductility , shown in in Fig. 3(0) Table 2, average 
Results for Near-Fault Ground Motion Records.—The importance of —_ 
“estimation of the LEDRS i is s also illustrated by the response of "Systems to > records 


grams s and theoretical studies 0, dM, and 37) indicate that ‘such ground motions = 
-may occur near the fault rupture, as well as at some other locations. While | 
Z such excitations are not critical for most elastic systems, unusually large ductility 
demands can develop in yielding systems when the duration of the pulse is - 
long relative to the structure’s natural period (7- 1-9). . For example, the derived | 
Pacoima Dam (DPD) base rock record G2) shown | in Fig. a) contains two 


= DUCTILITY, = ‘ 


TIME,SEC 


ACCELERATION AND VELOCITY COMPARISON OF REQUIRED 

FIG. 4.—Results for Newmark-! Hall IDRS Using he Pacoima Dam Base R 
Record (S16E): (a) Acceleration and Velocity | Time ‘Histories; (b) Comparison 


bound estimates (1, 12, and 28). For a a specified ductility of 4, the basic 


_Newmark-Hall IDRS is unconservative, especially for periods g greater than 0.4 
sec, [Fig. 4(b)]. For a specified ductility of 8, ductility demands of nearly 
30 were required in some period ranges. If the EPP SDOF systems considered 
were designed using a modified Newmark-Hall IDRS based on a V/A ratio 
@ 105 in. /sec/g (2, 670 mm/s/g), which corresponds to the actual accelerogram — 


when the specified “ductility is 4. However, when the specified Guim 
ductility is 8, use of increased design forces corresponding to a realistic V/A 
_ ratio does not achieve conservative results in the short period range. ; 
ATC ‘Type IDRS.—The tentative ATC design recommendations utilize a 
constant reduction factor over the entire period range, whereas the reduction _ 
factors for the Newmark-Hall type method vary (Fig. 1). For co comparable LEDRS _¥ 


ove 


i pu 
ock 
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in the short period range. The effect of this difference is examined below. 
_ Results for Standard Ground Motion Records.—The ground motion records 
listed in Tabie 1 were used to compute the inelastic response of EPP SDOF 
_ Systems designed according to the ATC IDRS shown in Fig. 1. In these analyses, 7 
viscous damping was assumed to be 5% of critical. As indicated in Fig. 5, 
the average displacement ductility values. are slightly less than the specifies — 
value of 6 for periods in excess of 1.5 sec. However, ductility demands increase 
significantly with decreasing periods, especially for periods less than 0.4 4 
The average ductility” over the entire period range considered is 9.54 (Table 
13.3 od: lo oft 
batten 
[SPECIFIED ~~. B 


or 
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EPP MODELS STIFFNESS DEGRADING MODELS 
5.—Response of SDOF Designed Using ATC IDRS with R =7a 
= 6 Subjected to Standard Records, € = 0.05: (a) EPP Models; (b) men 


\f 
2), and at a period of 0.1 sec an a ches ductility in excess of 60 ie obtained. 
5 While the coefficients of variation are only slightly larger than those at 
4 


7 and 


_ with the Newmark-Hall IDRS, the large average ductilities in the short and 
moderate period ranges result in very large variations in ductilities. For example 

for a period of 0.1 sec, computed ductilities ranged from 2.8 to more than 
130. While it may be possible to detail short period structures to develop such | ¥ 
large ductilities, and a number of factors not considered in these analyses, 
such as the use of minimum rather than probable material strengths in rar 7 
« conservativism in mong and detailing due to practical considerations 


| 
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imitations on anit sizes and strengths, contribution of nonstructural 
parsed soil-structure interaction, etc., could reduce the actual ductility os 
demands on short period structures, it is evident that a constant reduction factor, 

such as used in the ATC method, does not give the designer consistent control 3 
over the expected structural damage. 

INFLUENCE OF Viscous DampinG ON Ducturry oat 

Small amounts of viscous damping can substantially reduce the 1 response of 

elastic systems (7,28). The effect of damping on nonlinear systems is apparently 1 

accounted for in the ATC method in the reduction values specified for various | a 

types of structures, but no guidance i: is provided on the levels of damping as ssumed. — 

In the Newmark-Hall method, the effect of damping is directly accounted for ‘ 
in the construction of the LEDRS from which the IDRS i is derived. _ Thus, 
this method assumes that viscous damping affects inelastic systems in exactly a 

_ the same way as it affects elastic response. However, this is not the case, 
as indicated by the maximum displacement ductility demands shown in Fig. 
6 for EPP systems designed using Newmark- Hall IDRS and constructed using ~ 


Siig , derived for the basic ground : spectrum 1m shape and 2%, 5%, and 10% 


damping. Although a - ductility of 4 was specified in constructing each IDRS, | 
the design forces become increasingly unconservative as the viscous damping 
ratio increases. The tendency of such methods to overestimate the effectiveness _ 
of large damping coefficients has been previously recognized. For example, _ 
is Veletsos, et al., (41) suggests that the basic method not be used for pe 
ratios in excess of 5%, and modified reduction factors, including the effect 


+. damping, have been recently derived based on analytical results for different — 


ensembles of records (20,29). 


OF HysTeReTic Cuaracteristics of Ductiuty Demanos 
_ Actual structural systems seldom exhibit ideal EPP hysteretic behavior (21 
: 30, and 43). Variations from this ideal can have a significant effect on th 
-structure’s inelastic For example, 7 8 compass 
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FIG. 7.—Effect of Deformation H: Hardening and Softening on Ductility Demands—DPD 
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% 
FIG. 8.—Effect of Mechanical Mode! on Ductility Demand of Structures —— 


demands for ideal EPP models with those for bilinear mechanical models s having ; 


_ post-yield deformation hardening or softening equal to a fraction, p, of the ‘ 
structural ‘stiffness. The plots in ‘Fig. 7 are for structures with — 
4 Hall IDRS, in both cases ; subjected to the DPD record. The results wath 
that the effect of moderate deformation hardening, p = 5%, is only marginal 
- and generally conservative. On the other hand, moderate deformation softening, 
p = —5%, representative of geometric nonlinearities, tends to significantly ; 
increase deformations in all but the most flexible systems. For example, according © 
to Fig. 7, for y = 0.5, deformation softening structures with periods less than 
0.8 sec would collapse for this earthquake. In the case of deformation less than 
_ systems designed according to the basic Newmark-Hall IDRS, the DPD record 2 
results in significant increases in ductility demand for all periods and — 
_ would result for structures with periods between 0. 5 sec and 1.6 sec (Fig. _ 


‘pec carefully pr Pec when specifying a acceptable values 1es of ductility for ust nl 
with IDRS based on EPP systems. 
_ Structures can exhibit significant stiffness deterioration when cycled in the 
_ inelastic range. As shown in Fig. 8, the ductility demands for a stiffness degrading 
_— subjected to a particular ground motion can differ significantly from 
those obtained for EPP systems in some period ranges. However, it appears 
_ that on average, for the standard earthquakes listed in Table 1, the differences © 
are generally small. This can be seen by comparing the average ductility demands a 
for stiffness degrading systems with corresponding values for EPP systems, = a 
for r example, results for systems designed according to the Newmark-Hall DRS, 4 ’ 
_p = 4, are shown in Fig. 3 and in Table 2. Similar comparisons can be spor’ 
for the ATC IDRS in Fig. 5 and in Table 2. The displacement ductilities and 
: associated coefficients of variation obtained for stiffness degrading systems 
are both slightly smaller, on the average, over the entire period range than 
they are for EPP ‘Similar been made by other 
(29, 31). a) 


= 
CONSIDERATION OF OF Response 


While the maximum displacement ductility factor is a common and useful — 
index of the severity of inelastic response, other response parameters may be 
more meaningful in some cases. To assess this, information on various response 
_ parameters has been obtained for SDOF systems designed using the Newmark- 
Hall and ATC IDRS, e.g., see Fig. 9. Average values and corresponding — 
“coefficients of variation obtained for some of these response parameters using 
. the ground motions listed in Table 1 are shown in Figs. 3 and 5 and in Table 2. : 
Cyclic Displacement Ductility Factors.—When significant inelastic displacement 
reversals occur, the maximum displacement ductility factor does not realistically a 
_ measure the maximum inelastic deformation suffered by a system. Real structures. 
have limited capacities to sustain such cyclic deformations (6). To quantify 
magnitude of the largest inelastic excursion, a } cyclic displacement ductility 
- factor, Bes can be defined (21) as S indicated i in Fig. 10. Values of H< will — 


> 
pre 
> 
lm = s The average cyclic ductility developed by EPP systems designed using the J 


DESIGN 


Newmark-Hall IDRS is 14% | greater than the sidiands maximum displacement 

_ ductility for a specified ductility of 4 and 30% greater when the specified ductility 
; is 8 (Table 2). For the ATC IDRS, the differences are greater, especially for 

r Displacement and Drift Limitations.—Because of a desire to mitigate the adverse — 
ef ffects of geometric nonlinearities or to control damage to deformation sensitive — 

% onstructural components, it may be necessary to limit maximum displacements _ 
or interstory drifts. The curves in Fig. 9(b) clearly indicate that the maximum | 
displacements, here normalized by the peak ground acceleration, increase = f 
with increasing period, even though the maximum displacement 
_ ductilities are nearly constant or even decrease for periods in excess of 1 sec. 

Consequentiy, specified ductilities may have to vary with period to account 
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9 —Response El Centro Record (S00W)—Newmark- Hall IDRS, 
§= 5%: (a) (c) Yield Events 
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for r displacement limitations. For e: example, for a maximum allowable displacement, 
Urimie » the corresponding limiting displacement ductility, 


Pim: = limit ) 


average story drift index, 8, equal to the limiting displacement divided by the 


height of the structure, a Assuming T = aH®, in which a and B = empirical 


Inq practice, limiting displacements are often specified in terms of an — 


constants (36), Eq. 2 can be rewritten in terms of 0: 


EEC, 
— = 
) 
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in which y the maximum displacement of a multiple -of- 
system to that of an equivalent SDOF system. For a single-story structure, — 
equals unity, whereas for multistory systems with uniformly distributed mass 
_ Vibrating with a displaced shape varying linearly with height, y equals 2/3. 
= deformation limited ductilities are plotted in Fig. 11 for the case where 
= 0.025, B = 0.75, y = 0.67, and © = 0.01. The » values correspond to 
a IDRS in Fig. 1, and A equals 40% of gravity. For these values, relatively | 
, large ductilities can be tolerated for periods less than 0.3 sec. For systems — 
: "designed in accordance with the basic Newmark-Hall IDRS with a specified — 
ductility of 4, none of the maximum displacement ductility demands computed 
_ for the standard records in Table 1 exceed this drift limitation. However, the — 1 
"average ductility demands for systems designed using the ATC IDRS exceed 
the drift limitation for periods less than 0.2 sec and demands for some records — 
exceed me limit for up to 1.1 sec. drift limitations, rather 
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FIG. 1— —Comparison of Ductility Demands for EPP SDOF Systems with Ductility 
Corresponding Drift (a) Newmark- IDRS, = 4, = 5%; (b) atc IDRS, 
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The values of ductility than can be in design would be even 
_ than indicated in Fig. 11 for higher values of A or lower values. of 8. Sie 
Permanent Displacements.—Following an earthquake, inelastic systems may 
retain a permanent displacement. As indicated in Figs . 3 and 5 and in Table — 
_ 2, these permanent displacements (,) are generally a significant portion of 
_ the maximum displacement. For EPP systems, permanent displacements average = 
more than 40% of the maximum displacements. Moreover, the coefficients of 
ae associated with these permanent displacements average about unity— - 
easly twice the value associated with maximum displacements. Permanent 
_ deformations are somewhat smaller for stiffness degrading systems. In some — 
: cases, , large permanent displacements can make structural ‘repair infeasible and 
can accumulate during severe aftershocks leading to structural instability. 
. Number of Yield Excursions.—As shown graphically in Figs. 3 and 5, all systems 


De Leeann in these figures yield a large number of times. The average number 
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of ‘yield events tends increase with "decreasing pe: period, in very 
_ Short period range. Yielding occurs about the same number of times in each 
direction, even though the maximum displacement is usually considerably larger 
in one direction than in the other. This results in a significant number of yield a 
for EPP systems when large ductilities are " specified with the Newmark- Hall 
_ procedure, Fig. 3(c), or when large ductilities develop in the short period region — _ 
_ Energy Dissipation Demands.—Because of the numerous cycles of reversed 7 
yielding, the total hysteretic energy dissipated by a system can be significantly — 
larger than that inferred from its maximum displacement ductility demand or 
se from its maximum cyclic displacement ductility demand. A i 


_ comparative index of hysteretic energy dissipation is the equivalent energy 
_ dissipation ductility factor (21), » ,, numerically equal to the displacement ductility — 
of a monotonically loaded system that dissipates the same energy, and has 
the same yield strength and initial stiffness, as the actual system (Fig. 10). 
Analytical results are shown in Figs. 3 and 5 and in Table 2. For the systems - 
designed using the Newmark-Hall IDRS, Fig. 3 shows that relatively large 
dissipation demands are made on systems with periods between 0.2 sec and 
1.0 sec. Systems in this period range develop relatively large numbers of yield 
pe and reversals. The large number of yield reversals generally associated _ 
with the Newmark-Hall IDRS when the specified ductility is 8 [Fig. 3(c)] and 
5 with the ATC IDRS, especially in the short period range [Fig. 5(a)], also results — 
im large energy dissipation demands. Stiffness degrading systems dissipate about — 
_ the same amount of hysteretic energy as EPP systems, even though they fo 
7 not reach their full yield capacity nearly as often, as shown in Figs. 3 and | 
5. Energy is dissipated hysteretically by EPP > systems only when the full yield 
capacity is reached, whereas it is dissipated on nearly all cycles following the 
first yield for the stiffness degrading models considered. In view of the large 
energy dissipation demands that sometimes occur, caution should be exercised - 
in selecting ductility or reduction factors for use in the design of structures | 


General methods for ying seismic design forces for structures 
_ that can tolerate limited amounts of inelastic deformations have been reviewed. ; 
- The reliability of approximate methods suggested by Newmark and Hall (28) — 
and the ATC (36) for obtaining IDRS by modifying LEDRS has been evaluated — 
in a series of analyses 0 of idealized SDOF systems subjected to a variety of J 
_ The results indicate that inelastic response is particularly ‘sensitive to the 
actual excitation input as well as to the dynamic and mechanical characteristic 4 
of the structure. The scatter in peak response values obtained is so large that : 
the | | approximate methods considered « cannot be viewed as reliably limiting 
maximum ductility | demands _ to specified values, even for ideal EPP SDOF 
_ systems. Nonetheless, such methods can provide very valuable guidelines for 
4 design, if their limitations and relationship to the overall design process is fully , 
recognized. In particular, the large” scatter in results its reemphasizes the 
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for proper selection of structural sy systems and materials, and for detailing to 
- insure ductility capacities far in excess of the values specified in constructing 
Displacement ductility demands computed for EPP SDOF systems designed — 
_ using the Newmark-Hall IDRS with 5% damping were on the average less than: wid 
the specified values. The Newmark-Hall method became less conservative as 
the values of the specified ductility and viscous damping were increased. 
Differences between the values of A and V assumed in constructing the LEDRS | 
and those corresponding | os the actual ground 1 motion experienced led to very 
Computed displacement ductility demands for EPP SDOF sy stems designed 
using ATC recommendations also varied considerably with ground motion and 


values was not uniform over the period range considered. In particular, average 
- ductility demands substantially exceeded 1 specified values in the short period 


Results indicate that IDRS obtained from LEDRS tend to overestimate the 
effects of viscous damping on inelastic response. Stiffness deterioration and 
moderate amounts of deformation hardening generally appear to have small — 
effects on ductility demands of SDOF systems. Moderate amounts of deformation 
softening, however, were found to” increase ductility demands substantially in 
For certain types of structures and coumpencien, 2 response parameters like = 
maximum and permanent displacements, hysteretic energy dissipation demand, v 
and number and magnitude of inelastic reversals may have to be considered 7 
in establishing seismic design forces, in addition to the maximum displacement 
ductility. Results indicate that drift limitations can control the design of structures, * 
especially short period structures designed according to the ATC method. 
Moreover, permanent displacements averaged more than 40% of the maximum 
for EPP SDOF systems designed by either method. 
Additional studies are needed in several areas. In particular, it may be possible 
_ to reduce the scatter of response values by normalizing ground motions in terms 
of parameters other than peak ground acceleration, e. g., peak or incremental : 
ground — peak ground displacement, spectral intensity, Ayre’s intensity, 
ete. (2, 16, 17, 26, and 24). More reliable methods are needed to account for — 


Be ee in response associated with uncertainties in determining the initial 

A period of a structure as well as its damping and hysteretic characteristics. Methods © 
should also be developed to account for and ‘mitigate the potentially adverse 
- effects of deformation softening. The effect of severe, long duration earthquakes - 

4 intense aftershock sequences, or both, on design requirements should be © 


investigated. The reliability of available method for establishing critical seismic. 


design forces for multiple degree- -of-freedom systems should be evaluated. 
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: | structural period. The ability of the method to limit ductility demands to specified 2 
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effective peak ground acceleration; 4 

effective peak velocity-related scceleration; 

coefficient for estimating displacement, using Ref. m3 

seismic design coefficient, using Ref. 5; si shell 


“height of structure; 
system mass; 


- inelastic deformation capacity of structure, using Ref. 5; nee 
seismic coeflicient related to soil profile characteristics at site, using 


= fundamental period of system; 


parameter relating displacement of multi- degree-of-freedom system 
that of equivalent single- degree-of-freedom system; 


cyclic displacement ductility factor; 
equivalent hysteretic energy ductility 
= permanent displacement ductility factor; 


nee = maximum displacement ductility factor corresponding to attainment 
of limiting displacement; and 


= 
the nent ior 

_— af _ eat rate of post-yield deformation hardening or softening as fraction of 

oa 


= 
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CELLULAR AND VOIDED SLAB BRIDGES) 
By Baidar Bakht,' Leslie G. Jaeger,’ and Mo S. .Cheungy 
_ The terms ‘“‘cellular bridges’’ and ‘‘voided slabs’’ are used in this paper for 
& bridge types shown i in Fig. 1(a) and 1(b) respectively. In spite of the apparent © 
differences in the two types of | Structures, their structural behavior can be 
mathematically modeled by the same methods. nt 
_ The main feature which distinguishes these sti structures res from other common — 
_ bridge types is the deformable nature of their cross sections. The need for 
a method of analysis which will take into account the deformable nature ae 
7 cross section has been stressed by many (5,10,11,14,17,20). 
yay _ There is a general agreement that neglect of cell distortion can lead to serious — 
underestimation of live load longitudinal moments and shearsincellular structures. __ 
a In the case of voided slabs, however, opinions differ with regard to the influence q 
of cell distortion. Some maintain that transverse cell distortion in voided slabs | 


is negligible, while others contend that the analysis must take into account 
_ As described by Bakht, et al. (paper to be published in the Canadian — 
of Civil Engineering) there are several readily available rigorous methods of 
ps which can take into account the effects of cell distortion. However, 
many bridge designers in North America, being accustomed to the simplified | 


- methods of analysis of AASHTO (American Association of State Highway and 


: method has been edie for the analysis of cellular and voided slab bridges. 
This paper presents the basis of this method and its development details. — 
te method provides a measure of the increase in longitudinal moments and 
s due to the effect of cell distortion. It can be used in conjunction -— 
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(16), would prefer to analyze these structures by simplified means without 
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“cell distortion. 
The response to load ofa plate that is 
on two opposite edges is governed by eight parameters, namely D,, D,,. 
«=, ,, D,, L, and b. The development of a simplified method for such an n orthotropic 
is feasible if the parameters which characterize its behavior can 
be reduced to a manageably small number. Many simplified methods have been 
4 developed for determining the transverse distribution of longitudinal moments 
in bridges that can be idealized | as conventional orthotropic plates, by reducing» 


b. CROSS. oF A VOIDED SLAB BRIDGE 


FIG. 1.—Typical Cellular and Voided Slab Bridges 


the eight structural to following two dimensionless characterizing 


0.25 
implication of the two parameters is in Ref. 
13. This shows the transverse distribution of midspan longitudinal moments, 
“a M ,, due to a single concentrated load applied to two ‘‘different looking”’ bridges 


- whieh, , however, have the same values of a and @. The two structures were 


It is shown that the pattern of M, distribution is the same in the two bridges ; 
in spite of the differences in the aspect ratios of their planforms. —__ anita + 
_ Simplified methods based on a and @ are applicable to bridges that can be x 
realistically idealized as orthotropic plates. If there are several transverse - 
‘diaphragms in a cellular bridge, and if these are sufficient to prevent the cells 7 


‘from in the transverse direction, the structure can be idealized 
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increasing the transverse flexibility of the structure (in other words mm 
the ability of the transverse medium to distribute loads in the transverse direction). — 
Consequently, higher concentrations of longitudinal moments and shears _ 
those predicted by an analysis which neglects cell distortion. 
2 Consider a transverse slice of a cellular structure, simply supported at its 
extremities and subjected to a concentrated load as shown in Fig. 2. It con } 
_* readily seen that the increase in deflections caused by the flexure of a es 
and webs about their own axes, is related to the transverse shear rather than 
transverse moments. Therefore, “the representation of cell distortion by an 
ft equivalent reduction of the flexural rigidity of the transverse medium, as proposed _ 
some research , Aster and is likely to le to lead 


erroneous ous results. bang, bn 


Pat 


_ CELL 


CELL DISTORTION 
PERMITTED 


RESTRAINED 


bea be accounted for by an equivalent shear area in the transverse direction. This _ 
approach has since been adopted and verified by others such as Sawko (18), 
_ Hook and Richmond (12), Robertson, el al. (17) etc. Based on the premise _ 
a — structure can be idealized as a shear- weak orthotropic plate, 


_ Robertson, et al. (17) have | givena harmonic a analysis of shear-weak orthotropic 5 
plates, and have derived two characterizing parameters whose values depend © : 


B- the number of harmonics taken; this feature makes their proposed charac- 
unsuitable to be the basis of a simplified method. ee ¥ 


Ke 
NES 
— 2.—Effects of Cell Distortio 


... 
aie plates that are simply supported on two opposite edges. Seatac 4 
of these plates is characterized by three nondimensional parameters, a, 0, and 
_ 8. Where a and 6 are defined by Eqs. | and 2, and 8 which reflects the influence | 


of finite shear area in the transverse direction, is given by: 


in which, ‘the transverse shear rigidity = 1 the product of shear “modulus 


and the equivalent transverse shear area per unitlengthe 
_ _ Massonnet and Gandolfi (14) further tried to simplify the problem by reducing — 
the three characterizing parameters to two. This simplification was based on 
_ unrealistic assumptions, one of which was that the deflections, w, of a a shear 


weak plate can be divided into two parts w, and w,, and that the ratio, w, [Way 7 


— LONGITUDINAL MOMENT 
COEFFICIENT 


DISTRIBUTION COEFFICIENTS 


FIG. 3 —Distribution Coefficients | for Deflections an and Moments 
‘Temains constant at all points in the plate. The term, w,, is the deflection - 
of the plate with infinite shear area and w, is the change in deflections due — 
to the finite shear area. The fallacy of this assumption becomes obvious if 
the total deflections across the width of the plate are examined. The total area 
under the curve representing the transverse distribution of deflections must 
—- constant irrespective of the values of transverse shear areas and transverse — 


a fact not borne. out by the assumption of (w, /w3) being constant. en OK i 
_ Despite the shortcomings of the Massonnet and Gandolfi postulate with regard . 
to the simplification of the shear weak plate theory their derivation of the 
_ three characterizing parameters was based on sound assumptions. = 7 
a-6-5 Characterization. —Definition: If a If anumber of which are physically 


a 
| 
| 
| 
| 
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@aert the same pattern of distribution coefficients for longitudinal moments 
_ then a-6-8 is said to characterize the distribution of longitudinal moments. 
+ _ It is noted that the same three parameters may also characterize the distribution 
_ 4 other structural responses such as longitdinal shear and deflection. This does _ 
not imply that the pattern of deflection is the same as the pattern of longitudinal — 
a “moments and shears; each response has its own pattern and this pattern is > 
3 for that of a, 8, and 


‘TRANSVERSE STATION POSITIONS 
Since (ie 
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TRANSVERSE STATION POSITIONS 


- moments in a certain bridge. These coefficients were obtained by ‘the alrea already - 
_ well established grillage analogy method for cellular structures (11,18,21). It 
is recalled that the distribution coefficient for deflections is the ratio of the 
actual deflection to the average deflection. From Fig. 3 it is obvious that the 
patterns of distribution for longitudinal moments and deflections are different. 


t In what follows it is shown that each of these patterns is characterized by 


“ 
Fig. 4(a) shows distribution coefficients for longitudinal moments for single 


consonant loads and for vehicle loads, in three physically different bridges, 


were obtained by the tien the grillage analogy method. The 
characterization of longitudinal moments by a-6-8 is convincingly demonstrated 
in the Figure. For the same bridges the distribution coefficients for deflections — 
were found to be virtually indistinguishable from each other, thus demonstrating — 
that a-6-8 also characterizes the distribution of deflections, 


- In conventional plates, the of shears is 


— 


Since in cellular structure D, is very nearly equal to D, (13), the value of 
_ yis constant (= 1.0) for all cases. Thus, y need not be considered as a variable. 
The analysis of the three bridges considered earlier revealed that 5— 

characterized the ‘distribution of longitudinal shear with acceptable as accuracy. 7 
The comparisons of the distribution coefficients for longitudinal shear in the 
Bay bridges are shown in Fig. 4(b). It is noted that the characterization fails : 
for single concentrated loads, but is still valid for vehicle loads. — : 


FIG. . 5.—Cross Section of Cellular ‘Structure 
Bakht, etal. (paper to be published in the Canadien Journal of Civil 
after an intensive literature search have shown that the orthotropic plate . 
for cellular structures could as follows: 
2 
D, = 0. SEt, 


‘D,, = Gt, 


7 Poisson’ s ratio | and other | acigten is as defined in Fig. 5. . For isotropic | materials: 


| 
| 
06 Oe 
ley 


Pe, the expressions for plate parameters in Eq. 1: 


Since v = 0.15 for concrete, a becomes approcimetely equal to! 1.0 for 


- Thus the conclusion is reached that for concrete cellular structures | the oT 


a longitudinal shears for vehicle loads. bide uit A 

Devetopment Backaro ND TY 

Concept of Magnifier. Since there are two (namely @ and 

LD: 3), that affect the distribution of loads in cellular structures, the distribution _ 

s these loads could have been geome with the help of charts developed 


reported by Bakht, et al. (3), for the a-6 space of conventional orthotropic 
plates. Such charts, besides being large in number, could not have readily provided 
information on the extent to which cell distortion causes a cellular rend - 
to respond differently from its counterpart in which cell distortion is absent. - 
dt was decided to develop « coefficients in in the form of magnifiers which could - 


= to those goovided in the Ontario Highway Bridge Design Code (16) and : 


be used in conjunction with | those manual and computer-based methods that 
neglect cell distortion effects. The advantage of such magnifiers would be — 

they could be used in conjunction with existing simplified methods such as 
those given in the AASHTO specification (22) and the Ontario Highway Bridge 
Design Code (16), that they would directly the influence of cell 

ms _ For or cellular structures the conventional orthotropic plate theory corresponds 
to a characterization a = 1.0, @ = b/L, and oe = 0. For nonzero values of 


rr) Maximum intensity of M, if cell distortion ‘effects are JN : 


nn magnifiers corresponding to longitudinal shears are designated by Bd 


Maximum intensity of M, if cell distortion ¢ effects are absent 


Effect of Number of Vehicles.—A three-lane cellular bridge was analyzed by ad 
grillage analogy under single Contententes loads and under vehicle loads with a 


1030 
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and when ‘the: effects of cell, distortion 1. The bridge « corresponded | to 6 = 
- 0.35 and 8 = 2.15. The results, which are plotted in Fig. 6, demonstrated that 
| value of hon decreased with increasing eccentricity of 1oad from the bridge 
7 - center line and also with the increase in the number of lines of wheels. When 
ae longitudinal moments under a concentrated load increases, to maintain overall _ 


staties, it is accompanied by a reduction of longitudinal moments in other - |g 
pioveit locations. It can be readily visualized that under the action of several — 
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DISTRIBUTION COEFFICIENTS FOR LONGITUDINAL MOMENTS 


cowl 
"concentrated loads the effects near any one concentrated load are relatively 
& severe than if there was only one load on the bridge. In the extreme case 
of an infinite number of loads on a transverse section, i.e., a uniformly distributed 
~ Joad, the values of A,, and A, would both be equal to 1.0. The dependence 


me A,, and A, on the number of vehicles necessitated that values of the coefficients 


be developed f for various loaded lane conditions. 


‘Ranges of a, 8, and 0. —As described earlier, a is ie equal to 10 


) | 
| 
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firs widths smaller than the span. Thus the maximum value of 8 was assumed ‘ 
_ to be 0.5. For the lower limit of 6 it is noted that the span will not normally 


, exceed 10 times the width, corresponding to @ = 0.05. A lower limit of 8 = __ 
0. 025 was therefore taken | in | order to | be sure of accommodating even m truly — 


a could not be less ‘than corresponding, value. for | a solid slab ‘for 


q 


* Substituting the expressions for v, and D, 
0. .15, the expression is obtained for solid slabs: 


a The emallest values of (b/L) and (t/L) were aa to be (1/20) and (1/40) 
respectively, thus giving the smallest value of 8 as 0. 005. However, trials with 
ar configurations of cellular structures showed that 5 would never | 


smaller than 0.1, therefore 0.1 was assumed to be the lower limit. Similarly 
a was never found to exceed 8.0. It was decided to assume 10.0 as the upper 


Effect of Small Changes in Vehicle Positions. the position of loads 
has a significant effect on values d,, and = the variations of vehicle positions — 
within a lane have negligible effects on the coefficients. Analysis of a large 
number of bridges, _ effectively covering the full 6- -§ space, by the grillage — 
analogy method, showed the effect of moving design vehicles within the design 

lanes to have alessthan one percenteffectond,. $= | 

Effects of Changes in Lane Width.—The method specified in the Ontario 

| “Highway Bridge Design Code (16), accounts for changes in lane widths through 
a correction factor. In the present study a a series of four-lane bridges having — 
lane widths of 11 ft (3.35 m) and 13 ft (3.96 m) were analyzed by grillage 

rs analogy to determine the effects of variations in lane width on A,, and A,. 
It was found that in all cases the increase in lane width from uM ft (3.35 =) - 
a to 13 ft (3.96 m) did not change d,, and r, by more than 4.0%. Thus it was 

i concluded that the multipliers need "be developed for only one lane width, the 
results being applicable to bridges with other lane widths. 
_ Effects of Axle Spacings.—A design vehicle represents a large \ variety of vehicles 
on the road. Therefore, any method of analysis that is applicable to the design 
vehicle should also be applicable to any one of the vehicles that the design 
vehicle is supposed to represent. I If this is not the case, than either the design 

vehicle or the analysis is at fault. Bakht, et al. are among many who have 

7 demonstrated that the > simplified “methods ¢ of analysis for structures that can 

be idealized as conventional orthotropic plate, such as those given the Ontario kod 
Highway Bridge Design Code (16), are vehicle independent. For cellular structures 


the effect of changes in vehicles was studied by analyzing a typical — 


ming thatv= 
as 
il 
7 
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structure under different vehicle loads. It was found that the pattern of distribution = 
longitudinal moments under at the same transverse 
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8.—Meshes used in Analysis by Grillage Analogy 
m), width ‘of 36 ft (10.97 m), and a a equal to 1.0. ‘The Tequired \ values of -. 
and 8 were achieved by adjusting the values of the grillage beam properties. — 
For a given case, A,, (or A,), was calculated as the ratio of the maximum —T 
longitudinal moment (or shear) in the bridge with cell distortion and the 
3 corresponding value in the bridge in which cell distortion was neglected. 
q values of A,, and A, obtained from the of three-lane 


€X,, and A, Values were developed first for three-lane bridges Dy analyzing 
_ 44 bridges by the grillage analogy method. As shown in Fig. 7, the bridges nd 
- _ effectively covered the full 6-8 space. All bridges had spans of 70 ft (21.34 


“space. terms, and “ed for two- and four-lane bridges were 
by analyzing fewer structures, as shown in Fig. 7. All structures were idealized 
by grillages having 8 x 10 meshes as shown in Fig. 8. From the calculated © 


wis 
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23 FIG. 9. —aeaetiiins for Longitudinal Moments and Longitudinal Shears 
values | of the multipliers, _charts \ were drawn in the form of con 
and A, on the 0-8 ne. Thee charts are — in Fig. 9. 


ane Proposed method is in the form of adjustments a ‘should be made 


from both simplified and refined as age 
7 = _ The adjustments can be made either directly to ‘the maximum “responses or 
to the distribution factors, D, that are calculated according to some simplified © 


| 


= 
» displayed a gradual change as 6 and 8 were varied; thus eliminating the need —s 
| 
; ‘> the effect of cell distortion in the transverse direction. The adjustments, which | 
| 


: 
Design intensity of longi- Design intensity of 
tudinal moment, including | | = =m m longitudinal moment 
cell ignoring cel cell distortion 


| including cell 


- Eqs. 12 and 13 are omniscient. Since the value of a multiplier is always greater ah 4 
than 1.0, both these equations reflect the higher local intensity of longitudinal } 
_ moment when transverse cell distortion effects are present. Similarly f for or longitu- 4 
Design of longi- Design intensity of 
tudinal shear, including | = =r, longitudinal shear 
cell distortion effect ignoring cell distortion 


D for longitudinal shear we D for longitudinal 
including cell distortion shear ignoring cell 


is equal to. 1.0 for cellular structures, 
t ™ quantity within parenthesis on the right hand of Eqs. 13-16 can be obtained | 


TABLE 1.— —Critical Number of Loaded Lanes for Cellular ‘Struc 


‘ 


6=0.0 


by renting the swectere as a solid slab. The Ms which 
are functions of @ and 5 can be read from Fig. 9 for the number of lanes ‘ 


the bridge and the number of loaded lanes we) 


_ In order to apply the proposed method it is necessary to know the number 2 
of loaded lanes that govern the design moment or shear. This information is 
implicit in the AASHTO simplified method (22). Information about the critical 
ome of loaded lanes can be obtained from Ref. 2, which demarcates a 
_a-6 space for conventional bridges according to the critical number of loaded 
lanes. These charts correspond to the multilane modification factors specified — 

the Ontario Highway Design Code (OHBDC) (16). However, since 


| 

q 


OHBDC modification factors are not much different from specified 
in the AASHTO specifications, the charts should also be applicable to the 
_AASHTO loading. According to these charts the critical loading condition for a 
—2- and 3-lane cellular bridges should be as shown i in Table 1. Since the OHBD te 
_ Code (16) method is vehicle independent (3), it is it is ae that it be u used . 2 
for both the OHBDC and AASHTO loadings. 
The magnifiers can be used to investigate ‘the 2 applicability of methods thet 
_ ignore cell distortion in the transverse directions, to structure, such as voided 
slabs, for which it is not clearly decided whether the cell distortion can or 
— be ignored. If, for a structure, the values of the mangifiers are found — 
to be nearly equal t to 1.0 for the critical loading case, th then clearly the ase . q 


need to take transverse cell distortion into account. » 


 Bakht, et al. (3) recommend the following equations for the 


These recommendations based ‘upon a of the relevant 


ature (e. a 9 


“ bee 2—Values ot! F, for S, of Voided Slabs (Reproduced from Ref. 8) 


0.205 


10.—Cross Section o of Voided Slab 


= 
| 


= F,Gt a7) 


F,= obtained from Table 2 and other notation is as defined in Fig. 10. | 
‘tS Having obtained the values of the various parameters, the values of @ and 
8 can be obtained using Eqs. 2 and 3 respectively, and then the method is 
similar to that used for cellular structures. 
The value of a for most voided slab bridges ranges between 0.85 and 0.95. 
Although this is somewhat different from the value (1.0) for cellular structures, 
TABLE 3.—A,,, Values for Voided Slab Bridges 


17 


8 


= 


wne 


it was peers. that the va values of the magnifiers would be selatively insensitive 


be can also be used Tor structures in which the value of a is somewhat different 
than 1.0. This does not mean that changes in a do not affect transverse distribution 
om loads. It simply implies that the increases in the longitudinal Tesponses due q 


case. Thus and developed for ca cases with a a equal to 1. = 
4 


To investigate the effect « of cell distortion in voided slabs the values of \ ~~ 
were obtained for slabs of various void to thickness ratios and aspect ratios 
_of the planform. To obtain the largest, i.e., the worst, value of 8 it was assumed 

ee (L/t) was equal to 20 and (L/b) equal to 2. For the more er value 


— 
mum | Usual | loaded _maxi- | Am for | maxi- | Am for 
t,/P, | 8 lanes | mum | usual’ | | 
1.05 | 1.03 | 1.06 | 1.05 
0.70 | «0.70 | 0.96 | 0.32 4110 1.09 
| 0.70 | 0.77 | 0.26 109 | 107 | 
q 

| 
&g 


SLAB BRIDGES 
of of (L/t) was taken to be the s same as above but (L/b) wa 


assumed to equal to six. The values of together with the 
t 


on on 

exceeds 1.04 in three- and four-lane bridges, ingioininn ‘that in these "bridges 
the neglect of cell distortion can cause the ‘governing longitudinal moment intensity 
to be underestimated by at most 4%. For single-lane loaded conditions the 


underestimation can be as high as 17%. hes 
a When a bridge cross section is subject to distortion, the pattern of ee 
distribution of longitudinal moments requires three nondimensional perameters 
for its characterization. ‘In the case of cellular and voided slab ab bridges one 
oe. these three parameters has an almost constant value. The same is true of | 
: «es shears; hence the pattern of longitudinal moments and shears is 
characterized by the othertwo parameters, 
_ This paper provides charts, plotted as functions of these we parameters, 
which give the contours of a “‘magnifier’’ \, being the ratio of the maximum 
_ intensity of moment or shear with transverse cell distortion to that without 
it. It is noted that this “‘magnifier”’ or intensity ‘multiplier, always has a value — 
larger than unity, which depends upon the number of lanes in the bridge and we 
- : the number of lanes that are loaded. In a given case it is straightforward, by _ 
_ referring to the appropriate chart, t to ascertain whether transverse cell distortion — 


"scope of the simplified fits of AASHTO (22) and OHBDC 16) to bridges’ 
having transverse cell distortion. 
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sy — in this paper: 
= longitudinal flexural rigidity per unit width; 


longitudinal flexural rigidity per unit width; j 
», = torsional rigidity per unit width; 
= 


transverse torsional rigidity unit length; 


= = rigidity per unit length; 
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Es de. 
= ponding 
= distance between flange centroids (Fig. 4); 
| 
acing of longitudinal webs or voids; 


diameter of circular void « or height of rectangular void; 
thickness of top flange and bottom of cellular structure; 
= parameter defined by Eq. 1; ——— 
ie parameter defined by Eq. 4; 
= defined by 


magnifier for longitudinal shear; and 
Poisson’s ratio, 
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Static BEHAVIOR OF BEAM-TO- 


W. F. Chen,' M. ASCE and K. V. ‘V. Patel, A.M.ASCE 


The cost of fabricated steel construction is controlled to a large extent i 


the connections. Simplicity of connection details may have a greater influence — 
on the total economy of the structure than would a reduction in weight of 

_ the main members. Before 1950, r riveting was the primary y method of « connecting 3 
structural steel, but, in recent years, connections using welding or high strength 

— bolting, or both, are the most practical. Connections, which require a combination 
of welding and bolting, are considered economical. Previous investigations on 
_ these connections conducted at Cambridge, Cornell, and Lehigh universities | 


“were summarized and discussed in Ref. 2. Recent ‘results. c¢ conducted at the | 
University of California at Berkeley = were reported in Ref. 13. The moment 
. connections tested at Berkeley were Subjected to cyclic loading, ———s 


_ Described herein are tests of 12 full-size welded or bolted symmetrically-loaded 
moment-resisting beam-to-column connections, or both, subjected to monotoni- 
cally i increasing static > loading to failure. These connections are of importance 
in the design and construction of steel multi-story frames. All specimens were 
_ designed incorporating all possible limiting cases in practical connection design. 
This includes (1) All connections proportioned to resist a combination of high if 
shear force and plastic moment of the beam section; (2) all specimens are made 
of the American Society of Testing Materials, ASTM AS72 Grade 35 steel; 
ae) high allowable : shear stresses of 30 ksi and 40 ksi are | re utilized for ASTM 
A325 and A490 bolts in bearing-type connections respectively; (4) some png a 
are designed with attachments having oversize holes or slotted holes, or both; 
and (5) some specimens are designed using friction-type joints and others — 
bearing- -type joints. This paper summarizes the highlights of test results of this | 
experimental program. In addition to the behavior and the manner of failure ; 
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oof ‘the beams and their connections to the columns, the implications of the | 
results on design procedures for safe, efficient, and economical beam- to-column 
connections receive particular attention, 
In 1971, atest program was initiated at Lehigh University to study the behavior 
and strength and to develop the design method of moment resisting beam- to- 
column flange connections (7,8). These are connections in which the beams 
are framed into the column flanges with the _ beams causing bending of the 
 cokmea about the major axis. The test program, consisting of 12 full-size 
_ beam-to-column connections, was under the guidance of the Welding Research 
Council (WRC) Task Group on Beam-To-Column Connections. The investigation — 
of the behavior of such connection is being done to aid the designer in his ‘Ss 
decisions on their use and to improve design rule. Detailed results of this _ 
e experimental program may be found in Refs. 3, 8, and 15. Some of the test 
results together with a preliminary theoretical analysis have been reported © 
~ elsewhere (6,11,16,17). A detailed elastic-plastic-fracture finite element analysis _ 


@ these specimens is currently underway at Purdue University. — isto? aq) #0 


The specimens were « designed a according | to the American Institute of Steel 
Construction, . AISC specification (10). The connections we were proportioned | to. 
esis the plastic moment of the beam section. Since the loading condition 

resembles gravity type loading (dead load plus live load), the load factor used 
Se 1.7. The stresses used in proportioning welds, shear plates, and top and 


bottom moment it plates were then equal to lI. 7 times those given in Sec. I. 


‘connections the design gn shear § stresses used v were paren to pe 7 times es 30 ksi ol 

40 ksi, respectively, instead of 22 ksi and 32 ksi suggested in the current — 

_ specification. The concept of this procedure will be discussed later, oa ot 

_ Member Size and Beam Span.— —The connection specimens were chosen to” 
have an appropriate combination | of a beam section and a column section that 
_ Tepresented the real interior beam-to-column connections in a multi- story frame. 
Three different sizes of | specimens consisting of Wi4 x 74, W24 x 61, and 

w27 x 94 beams connected to W10 x 60, W14 x 136, and WI14 x 176 columns, — 

respectively, were used in this test program. All of the beam sections are plastic 
design sections that satisfy the requirements of * ‘minimum thickness’’ (Width- 

‘Thickness Ratios), of the AISC Specification. Both w24 x 61 and W27 


— «94 shapes are the lightest in weight in each par particular group given ‘in the AISC | 


‘Another factor considered in selecting beam sizes is the way a wide- flange 

_ shape resists bending moment and shear force. It is wel known that the flanges — 

resist most of the bending moment, and the web | almost entirely carries the : 
shear force. The ratio of flange : area to web area . furnishes a an index to the 

- amount of moment carried by the web, which must be transferred to the flanges 7 
at the connections since the shear connections have negligible moment-resisting | 

capacity. The ratios of one flange area to web area, A,/A,,, for W14 x 74, 

W24 x 61, and W27 x 94 sections are 1.39, 0.44, and 0.60, seepectively. The € 
behavior of these sections be representative of a wide range of of wide- 


sections. 


| 
| 
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BEAM-TO- COLUMN CONNECTIONS 
ot was | deed to use column sections that did not need horizontal stiffeners. _ 


_ In each case the « columns selected were the least column size based upon the 
AISC tension or compression tng criterion. The s specimens were proportioned _ 


the factored shear capacity of eal pris bolts in beam web would be Sadhed 
a concurrently. Beam span then was simply the ratio of moment to shear values. 
a All test specimens were fabricated from wide-flange shapes made of ASTM — 


A572 Gr. 55 steel. High strength steel was used because there is a narrower 

: margin between yield and ultimate than | for lower strength steels. Further, the 

"plastic range or ductility for A572 steel is somewhat less compared to A36 
_ steel. Thus, if the connection behavior is a for A572 steel, the results - 
could be assumed to apply to lesser r grade steel. A detailed report of of material 

_ Fasteners and Holes.—ASTM A325 and A490 bolts were used to > assemble 
the joints. In bearing-type connections, the allowable shear stresses used in 
design for A325 and A490 bolts were 30 ksi and 40 ksi, respectively. The use 
of higher allowable shear stresses reflects the logical design criterion that would 
4 result if an adequate factor of safety were applied against the shear ah > 


of the fasteners. This design criterion based upon the results of a study 
7 of A7 and A440 steel lap and butt joints fastened with A325 bolts, and A440 F 


steel joints connected with A490 bolts (5). Tests have been subsequently carried 4 
out to substantiate the suggested design eqhesion, expecially the use of A490 rm 
bolts in A440 and A514 steel joints (9,18). 
Since both oversize holes and slotted holes are desirable to facilitate erection 
adjustments, and slotted holes may better facilitate the assumed distribution * 
of shear and moment at the connections, experimental justification is required 
f or beam-to-column connections assembled with high-strength bolts with enlarged — 
and with slotted holes. Previous research has indicated that oversize holes, 
‘sized according to bolt diameter, do not adversely affect the slip behavior of — 
_friction- type joints or cause undesirable bolt- tension losses (5). It was also 
observed that slotted holes did not affect the strength of bearing- type joints. — 
the test program, | 1/4 in. _ round holes were used in top and 
moment plates fastened with 1 in. diameter A490 bolts and "designed as a 
4 friction-type connection (Test C8 in Fig. 9). The use of 1/4 in. ‘oversize holes 
‘is the maximum size permitted in the current Specification (14). Slotted holes 
- were used in one-sided shear plates fastened with either A490-X bolts (C3 in 
_ Fig. 4) or A325-X bolts (C8 and C9). ‘The 1 remaining joints had round holes” 
1/16 i in. larger than the nominal diameter of the bolt. 
Both A325 and A490 bolts were installed by the turn-of-nut method. Washers — 
were not used for A325 bolts. In bearing-type connections, A490 bolts had 


“both the head nut. Nut rotation ‘from snug ‘tight condition was 2 ‘turn 
as required by the specification (14). All bolts were calibrated and installed — 
in the Fritz Laboratory, Lehigh University, | 
_ Welds.—The connection specimens were welded according to the American 
"Welding Society AWS Building Code (4). The welding process used for groove 
welds was the innershield procedure; the electrodes were E70TG (flux cored — 

a arc welding with no Bas shielding). The type of filler metal fee. beam 
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groove weld 
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groove weld 
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round holes aie. 
W 


W27 x 
Wi4 x 176 
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| W14 x 136 


M, = 8, 360 k-in. 
"moment plates 
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Fi in 1 1/4 in. 
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T COLUMN CONNECTIONS 
BLE 1 


Ti groove weld | beamweb 


= 374k (94.7% 


groove weld 


flange | groove welds in the flat position for moment ‘groove: welts 
_ in the horizontal position was NR-311; NR-202 was used for beam web groove 
7 welds in the vertical- -up position. The electrodes for fillet welds were E7028. — 
In determining the size of fillet weld, the allowable shear stress on the effective 
Nondestructive testing methods were employed to inspect the welds whe 
testing of the specimens. Groove welds were inspected by ultrasonic testing 
r g fillet welds by magnetic particle. Results of weld inspection were crabeted 
according to the AWS Code. Those rejected were repaired and atte 
inspected priortotesting, = | 
Descriprion OF Srecusens 
test n, designated as is summa n Table 1. Detailed 
- _ descriptions of each test are given as follows. fey 
The joint detail of Test C1 is shown in Fig. 1 were directly 
welded to the column flanges providing for plastic moment capacity. A one-sided 
shear plate fastened with three | in. diameter A490-X bolts was used to resist 
ae shear. The fillet weld connecting | the shear plate to the column flange 
was sized for vertical shear + only; the moment due to the eccentricity of the 
applied load was neglected. Horizontal stiffeners were designed according to | 
the AISC specification. roared the connection panel zone was under symmetrical — 
loads, a clearance of 1/2 in. was provided between horizontal stiffeners and 
_ column web. The size of fillet welds for horizontal stiffeners was determined — 
by computing the force taken by the stiffeners when plastic moment was attained 
_ Fig. 2 shows Test C10, which is a fully-welded connection and was used — ‘ 
‘as a control test. Beam flanges and beam web were connected to the column — 
a flanges by groove welds. An erection plate was tack welded to the column 4 
_ flange, and was used as the backing strip for the beam web groove weld. Test 
C2 is shown in Fig. 3. Its connection type is similar to Test Cl, the only “a 
difference being that horizontal stiffeners are not required. 


_ Test C3 is identical with Test C2 except that the one- -sided shear plate of 
Test C3 (Fig. 4) has slotted holes. The dimensions of these slots conform to 
provisions in the current specification (14). A continuous bar 5 / 16 inch in thickness 
2 and having a width equal to the length of the slot was attached ‘on the side ? 
of the slotted shear plate. (The addition of continuous bars aad single shear 
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BEAM-TO- COLUMN CONNECTIONS bin 


connections was ap by the Research | ‘Council on Riveted and Bolted 
Structural Joints at its annual meeting on May 12,1971), 
_ Fig. 5 shows the joint detail of Test C4. Moment capacity was provided 

"through direct groove welding of the beam flanges to the column flanges. Vertical 

Shear was resisted by at two-late welded stiffener seat designed according to 
Table 8 of the AISC manual. The. strength of this connection should be greater 


than that of Test C5 shown in Fig. 6. 
In the case of Test C5, the beam was connected to the column by groove 
welds only. It had neither an erection seat nor an erection clip. The purpose __ 
4 was to determine the actual capacity of the beam flange groove v welds. 
_ Fig. 7 shows bolted top and bottom moment plate connection Test C6. T 7 
___ plastic moment was carried by flange plates that were fastened with | in. diameter _ A 
—A490-X bolts. The design procedure follows the example of the AISC Manual. 
The bracket stiffeners were designed with the aid of Table 8. t# ratiasasao od} in 
Test C7 is - shown in Fig. 8. The vertical shear is resisted wy a one- sided 
shear plate connected to the beam web by three | in. . diameter A490-X bolts. } 
- Tests Cl, C6, and C7 were designed for the same queues of on of moment and 
_ Shear, and, therefore, their behavior should be comparable. tesT 
holes in moment plates. The use of 1 1/4 in. round holes for J in. _ diameter 
_A490-F bolts is permitted by the specification (14). There is no reduction i in 


ss Test C8 (Fig. 9) was designed as a friction-type connection — over-size : 


slip 1 resistance of the joint. The one-sided shear plates had slotted holes 
were designed as bearing- type connections. A continuous bar was also attached _ 
on the side of the shear plate (Fig.9), 
- Test C9 (Fig. 10) is similar to Test C8. For the purpose of comparison, the 
moment plates of Test C9 were designed as bearing-type connections Roving 
ound holes 1/16 in. in excess of the nominal bolt diameter. sp Ie » stent 
All 12 tests of the series of 12 specimens have been completed. ’ Table” 7 
4 shows all specimens along with beam and column sizes, amount of moment 
~ and shear they are subjected to, stiffening requirements, and beam - span. The . 
‘ type of connections studied were fully welded, flange welded with various means : 
_ of carrying the shear load, and fully bolted. Table 2 contains a brief description 
of each failure of the 12 specimens. The 12 specimens represent typical ay 
“4 _ connections for the lower stories (W27 x 94 beams connected to W14 x 176. od 
column), middle stories (W24 x 61 beams, Wi4 x column), and “upper 
Stories (W14 x 74 beams, W10 x 60 column) of a multi- “story frame. Each 
group includes a fully groove-welded connection to act as a control specimen — 


_ The discussion of these tests are presented in the following. The le 
divided into four parts we simil odd ea ai 
Fully welded connections (C10, C11, C12). 


a Flange welded, web bolted connections (Cl, C2 and C3). watt) a 


Fully-bolted connections (C6, cI, C8, , C9). ows 


| 
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“which must be ‘examined to insure an pr dip ponds 
— Welded Connections.—The fully welded specimens serve as control — 
3 specimens for evaluating the performance of several other connections of different 
_ joint design in the series. _ These fully welded specimens are considered the 
ideal connections” becouse they the connection of directly to 


The first specimen te tested and fully analyzed was (Fig. welded 
connection of the | largest | size considered in this study. The sections used in 
‘this connection were a W27 x 94 beam and a W14 x 176 column and were : 
groove welded. The load-deflection curve of Test C12 is shown in Fig. 12 and 
the connection at o end of test is shown in ‘Fig. iB. bb Fig. 12, P,, is the 


Crack at the fillet weld of the tension stiffener { 


Tearing of column web along web-to-flange juncture © a med 
Fracture occuring at the heat-affected zone of the tension flange groove weld 
Buckling of beam web and fracture of beam web initiated at the cope hole © ta 
Fracture of beam web initiated at the cope hole 


Failure of five bolts in one of the tension 
Buckling of column web and test terminated 

Six bolts on one of the tension flange connections and one of the wae bolts res 1? 

Crack at the fillet weld of the tension stiffener 
Crack at the heat-affected zone of the tension flange groove weld ts nageils 
Fracture at tension flange groove excessive we deformation 


maximum load of the connection under test, P, is the theoretical plastic limit 
load, and P, = =. -P,/1.7. The connection n behaved linearly elastic until the first 7 
yield lines began forming i in the compression web of the column. As the load 7 
increased, the yield spread in column web compression as well as tension regions. 
The testing was concluded due to a combination of excessive column web : 
formation and fracture of weld at tension flange and along beam web (Fig. 
14). From load-deflection curve, it is evident that this type of connection cin 
be used in plastic design as the plastic limit load, sufficient rotation capacity, © 
and adequate elastic stiffness are developed. The present AISC Specification 
(10) provides adequate rules in design of such welded connections. The extensive r 
yielding in tension and compression regions of column web as well as — 
of beams shows a balanced failure (Fig. 13). Detailed r results and ew of “_ 
_ The other two fully w 


: _ The overall behavior of these connections as characterized by their load-deflection 7 
s =€©=—-—s curves will be discussed along with the modes of failure. The load-deflection | 
Test 
(1) | | 
C3 
C4 
co 
C10 | 
C12 


(Fy 


} & 


- Strain Hardening 


| 
| a 


_ horizontal stiffeners were provided in the column wal for this specimen. The 
- sections used for Cll were W24 x 61 beam and WI4 x 136 column. The 
load-deflection behavior of C10 and Cll (of different sizes) is compared to im 
C12 on a non- -dimensionalized plot as shown in Fig. 16. Their overall behavior — : 
is similar, as well. These three connections are used as the control — 


in a their respective size groups to which all results are pects In Fig. ‘16, 


J 

oi 


a= 


= In summary, fully welded amine have the m merits for plastically designed — 
‘multi- -story frames if designed according to AISC Specifications. They maintain 
elastic stiffness up to working load, and starts loosing strength due to local — 
: ‘yielding. Lateral deformation of the column web in the compression region 
- causes unloading. The failure of the connection was due to a combination of 
. excessive column web deformation in the compression region and fraction of fy 


weld at tension flange. The inspection are sho own 
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13.—Load-Deflection Curve of Specimen C1200 
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FIG. 16. —Load- Deflection Welded Connections 


FIG. 15. Weld Fracture at 
at Tension Flange of C12 the 
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; on the beam web. This weld can 1 be “expensive when done as a field weld. b 
Cl 1) a bolted web plate to in 


WI0 Column 
A572 Gr. 55 


Hardening Consi 


4 FIG. 17. —Load- Deflection Curves of Specimens C C1 and 
Py. 
_ 
142176 
Ge 35 


and 


a 18. —Load-t Deflection Curves of Specimens C2, C3, and C12 fee 

place of the vertical weld, and behaved ve wt much like C10 (Fig. 2), its couteel, 
as shown in Fig. 17 where M, is full plastic moment. Both connections ~ i 

; horizontal stiffeners. A good correlation between predictions by simple plastic 

Ga 4 and results is obtained. This is is due to the the use of of 


| 
= 
— 
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horizontal stiffeners the rigidity of panel zone, 
assumption of the plastic analysis. 
The Specimen C2 (Fig. 3) and | G3 (Fig. 4) had i beam flange directly w wanes 
. to column flange and the beam ‘web was bolted to the shear plate. A 490 high 
4 strength bolts with shear capacity of 40 ksi in bearing were used. Test C2 7 
_ and C3 are of the same size, only difference being C2 shear plate had round 
a holes while C3 had slotted holes (Fig. 11). At the predicted plastic limit load, 
the bending moment was assumed to be carried by flanges due to strain- -hardening, : 
and t the shear force w was , assumed to be resisted by the web ¢ attachment. This “4 . 
_ assumption is examined considering the joint details used: (1) Beam web connected nf 
_ to column flange by groove weld; (2) beam web shear plate fastened by 
high-strength bolts in round holes; and (3) beam web shear plate fastened by 
-strength bolts in slotted holes. 
Figure 18 shows the load- deflection curves of C2, C3 and C12. Both C2 
C3 | showed adequate. ‘elastic stiffness under working load. The deviation 
of C2 and C3 from C12 was due to slip of joints that occurred above working — 
load. C2 failed due to tearing of the column web along web-to-flange junction 
while C3 failed when fracture occurred at the heat-affected zone of the groove 
_ weld at the tension flange. The bolts _— ae 0 the the complete 1 redistribu- 
of forces at maximum load. 

Tests demonstrated that the flanges were able strain harden: sufficiently 

te transmit the full plastic moment of the beam section. In Fig. 18, P,, is the _ 
plastic limit load counting flanges only; P, corresponds to full strength, and 

-P,, is the plastic limit load modified to include the effect of shear force. Quite — 

both connections ns C2 and C3 were able to strain harden sufficiently 

Tests also showed that connections slotted holes (C3) and round holes 
(C2) exhibited similar overall behavior. The use of slotted holes is desirable — 
to permit erection adjustments and may facilitate the assumed distribution of 
shear and moment at the connections. Detailed results and ‘ee of e 
part of study have been reported it in Ref. 6. i oe ie 

of beam « connections that are € subjected toa . shear force at the predicted onl 
limit load of not more than 60% of V,, where V, is the shear Saves that a 
full yielding of web. If the shear force is approximately equal to V,, the proposed 
- theory in Ref. 6 may be used. The slotted holes may be used in one-sided — 
"shear plates that are designed as bearing-type joints. Also, the column web_ 
: ea requirements developed for A36 steel are applicable for higher ‘yield 
levels (up to 55 ksi). The fillet welds connecting a shear plate to the column 
_ flange may be sized for vertical shear only and the moment due to the eccentricity _ 
of the applied load may be neglected. 
s Flange Welded with Various Means of Carrying Shear.—In this group, ~ 
. connection included are C4 (Fig. 5), C5 (Fig. 6) and its control test C12 (Fig. — 
11). In all three connections, the beam flange is welded to the column flange. — 
The primary difference among the three is the method by which shear is carried. Es 
The control specimen utilizes the beam web to carry shear, the C4 specimen 7 
carries shear by means of a beam seat having beam web unconnected, and 
the CS carries both moment and shear in the beam flanges. This three specimens 
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Test C4 failed with the beam buckling Fie. 
connection was then redesigned to include a beam web stiffener, labeled Test | 


+ 


* 


th 
He 
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C4R. In the redesign, the web was considered to act as a column with one 
end hinged and one end fixed. For connection the 
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BEAM-TO- COLUMN Ci CONNECTIONS 
¢ 3 toa sheds ll fracture, which initated at the cope hole of t penis tension —_ 


side. Testing was terminated when the tearing of the beam web near the k-line — : 
became excessive and the beam web buckled near the beam compression | flange 
_: load-deflection curves of C4, C4R, C5, and C12 are > compared (Fig. 
22). The fully welded connection (C12), and the stiffened seated beam connection 
> beam web stiffeners C4R are basically identical in their general chovier’ 
of the applied loads, and may be used interchangeably. The flange-welded only 
7 connection C5 attains 51% of its predicted plastic limit load based upon the 4 


. 21 ieiials C5 at End of Test 


COR (C4 


“FIG. 22.—Load- Deflection Curves = Limit Load; P 
Load, P, ae 7) (1 in. in. = 26.4 25. 4 mm; 1 kip = =m pure 

Detailed 


stiffened beam connection as C4 provided the possibility of the buckling of 
the beam web is checked and beam web stiffeners are added if necessary a 
C4R. Although the stiffened ‘Seated (C4) fails by excessive 


capacity, though it attains only about half of its predicted plastic limit load. 


| | 
— 
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Te 
Bolted Bolted Connections. .—Herein, the results of t tests of those specimens 
designed with high strength bolts for resisting moment and shear force are ; . 


a This group includes specimens for each of the following cases: (1) 


- Flange bearing bolted with a stiffened beam seat; (2) flange bearing bolted, 
web bearing bolted (both of the smallest combination considered; W14 x | 
beam section and W14 x 60 column section; (3) flange friction bolted, web 
bearing bolted; and (4) flange bearing bolted, web bearing bolted, (both of the 
-—Thewe comne size tested, W24 x 61 beam section and W14 x 136 column section). | 
These are labelled as C7, ues , and C9 respectively. 2 


23. — Overall View 


C6 (Fig. and C7 8) are designed to resist moment through 
~ bearing connection using eight A490 bolts. The moment plates are groove welded a 
q the column flange. Horizontal stiffeners are provided in the column web, 


by a stiffened beam seat fillet welded to the column slange for C6 while shear j 7 
is carried by three A490 bolts in the beam web for | C7. The control test. f 


6 and C7 is C10 (Fig. 2). ab or 
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BEAM- TO-COLUMN CONNECTIONS 


‘The C8 (Fig. 9) is "designed to resist moment using ‘fourteen 1 1 in, diameter 
nay a friction type bolts having | 1/4 in. oversize holes in the moment plate. 

- ‘The moment plates are groove-welded to the column flange. A bearing connection . 
using / A325 bolts in long slotted holes is used to resist shear. The C9 test (Fig. 

— 10) is similar to C8, but for the pu: purpose of comparison was designed as a el 
bearing type connection for resisting moment, meaning that fewer bolts (six o 
A490 bolts) were required in each beam flange. The control test for C8 and 
C9 is the welded test specimen C11 (Fig. 15). Sw! 

ow 


FIG. 25.—Prying Action of C7 Moment Plate 


— 


P 


4 CO-Fully Welded 
C7-Fully Bolted 


C6-Flonge Bolted w/ 


Fig. 23 shows overall view of C6 at failure. The failure of this test occurred — 
when the eight bolts connecting the tension flange to the moment plate sheared 
Off simultaneously. Fig. 24 shows overall view of C7 at failure. The failure 
: was due to the failure of five bolts in one of the tension flange connections. _ 
Two bolts on the outer line of a tension beam flange failed first. The bolts = 


- appeared to fail i in tension. A close look at the moment “say shows it bending — 


| 
cL 
| 
siding of beam web up to the end of the moment plates shows that plastic : 


nian, , if it it forms at all, will fo form at the end of pany moment plates and and not 
at the column face (Figs. 23 and 24). yor 
_ The load-deflection curves for tests C6 and a along with their control test _ 
C10 are shown in Fig. 26. Test C10 closely follow the predicted stiffness, -— 
then deviates from it in a smooth curve due to yielding. Initially, C6 and C7 
- follow the predicted stiffness, but at a load of approximately 150 kips, both 
tests abruptly display a smaller, linear stiffness. This is presumed to be due 
to the slipping of high strength bolts in the moment connections into 
_amdsome yieldingtakes place, 
‘The load deflection curves for C8, C9 and their control test Cll are shown — 
in Fig. 27. Test C9, the bearing-type of connection, follows the prediction until 
it suddenly follows a smaller slope due to slips. Test C8, the friction connection, 
- does not display the distinctive slope due to the slipping of bolts into bearing 
but rather follow smooth curve. The ‘initial stiffness of C8 is - slightly greater 
_ than predicted due to the longer moment plates needed for a friction connection. _ 
- Detailed results and analysis for this part of study have been pees | in Ref. 5 _ t 


© C9-Fully Bolted- 
a ce- Bolted - Friction 


FIG. 2. —Load- Deflection Curves—C8, co + De 


‘than the plastic- -limit ‘load of the beam. This is in contrast to prior tests on a 
moment connections with beam flanges welded to the column. The deflection a 
resulting from slip of bearing bolted moment connections may be an additional J 
factor to be considered in the analysis of the stability of frames. The bolted 
moment connections tested exhibit a maximum ‘strength that is at least equal - 
to that of the comparable welded connection and varies up to 30% higher than = 
the welded connection due to the strengthening effect of the moment plates. 
Bolted moment connections also exhibit sufficient rotation capacity when 
compared with welded connections. For the range of beam and column shapes 
tested, the results have shown the adequacy of the present stiffening requirements 
& the AISC Manual. When column web stiffening is required, the test data 


indicates that such stiffeners need not be welded to the column web for — 


| 
summary, the moment designed for bearing 
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BEAM. TO. COLUMN CONNECTIONS 


-Conciusions AND Summary 


series full- scale beam-t -to- column moment- resisting flange connections 


sity 2 are porns: here. It shows ‘the general ie of moment connections 
_ and peculiar change in behavior when the connection detailings are changed 
from fully welded to bolted and welded and to fully bolted. This testing program - 
ong with future theoretical and analytical work will provide a thorongh 
understanding of such connections and will lead to recommendations and 
guidelines for those involved in their design. gy 
_ The results of these connections will provide needed information for the 
analytical work for the critical parts involved in these connections. By using 
Rg these results, bench marks can be set and connections analyzed with numerical 
_ methods like finite element method using high speed computers. Thus these _ 
results will i help to develop analytical capability needed for design of these 
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DESIGN OF SLABS FO FOR CONCENTRATED Loaps 


By y Hans G Gesund,’ F. ASCE 
liiw teal besylans od of esbom odi al 

g * The Yield Line Theory i is a well established (4, 8, 9,10) method for limit analysis 4 


or design of reinforced concrete slabs. Solutions are available for isotropically 
and orthotropically reinforced slabs, supported with a variety of boundary — 
conditions, and subjected to either distributed or point loads (4,6,9,10). A few 
solutions are also given (8) for combinations of p point and distributed tow 
acting on isotropic slabs. 
_ Superposition in general is not valid in inelastic analysis, though it can be 
shown (4) that addition of the results of Yield Line analyses for individual — 7 
loading cases will Bive a lower bound on n the total bound load 
load can be determined in this manner depends on how nearly congruent the | 
4 individual oer mechanisms are. When concentrated loads are placed near J 
= : slab boundaries, the concentrated and the distributed load mechanisms may > 
b 2 be very different, and Superposition v will give poor results. An example will 5 
_ Reinforced concrete ‘ade: are normally constructed with orthotropic rein- 
; forcement and different positive and negative bending moment resistances, as : 
dictated by serviceability requirements and designer judgment. The reinforcement _ 
is also likely to be different near the boundaries and in the interior of the a 
+ and | under ‘some boundary conditions it may be difficult to provide om naqgedl 


7 in | any” of the yield hinges. The accuracy, however, with which the collapse 


: ‘the boundary. Concrete slabs also tend to have considerable self —" which 
wil reduce the magnitude of any concentrated load a slab can carry. cae | 
oo True point loads do not exist. All real ‘“‘concentrated’’ loads must actually 
be spread over some area. _ Depending on circumstances, the critical failure 
mode may then be either flexure, or “punching” s 4 shear. The shear problem > 
will not be addressed here. The punching shear provisions of the American 
Concrete Institute (ACI) Building Code (1) appear to provide adequate guidance © 
for safe design. The flexural problem will be addressed. There is also the possibility 
Prof. of Struct. Engrg., Univ. of Kentucky, Lexington, Ky. 40506. 
Note.—Discussion open until February 1, 1982. To extend the closing date one month, 
a written request must be filed with the Manager of Technical and Professional Publications, — 
ASCE. Manuscript was submitted for review for possible publication on July 25, 1980. — 
This paper is part of the Journal of the Structural Division, Proceedings of the American 
nt of Civil Engineers, Errata Vol. 107, No. ST9, September, 1981. ISSN 0044-— 
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= 
of interaction “between shear and flexural failures. The ACI Co Code f provisions 7 
deal with this - only indirectly, but it seems unlikely that flexural distress will 
_ enhance shear strength. It is therefore doubly important to aac an accurate 
assessment of the flexural strength. 


to be distributed over a small rectangular area of dimensions c and Xe, oriented 
and located with respect to the coordinate axes and slab boundaries as shown > 
in Fig. 1. The entire slab, including the area beneath P, will carry a uniformly 


dies 
J daidve & woul cogent 


= In the models to be analyzed here, a “concentrated” | load, P, will be assumed 


x 


bie FIG. 1 1.—Plan View of Slab and “Concentrated” basige od 

- distributed load of magnitude, w, per unit area, which may be self weight plus 
additional likely surcharge. It will, however, be assumed that << _ 

-Johansen’s (4) stepped yield criterion, , which has been validated by Lenschow - 

a and Sozen (5) and by Mills (7) among others, will be used here. The flexural — 

strength of the slab can then be expressed as follows with reference to the 

coordinate axes shown in Fig. 1: m = unit yield moment per unit length; k,m 

= negative yield moment per unit length provided by top tension reinforcing 

placed parallel to > the x-axis; m= positive yield m moment unit 


a 4 
| 
aly 
I / _ yield moment per unit length provided by top tension rebars placed parallel __ 
i 


CONCENTRATED LOADS 


= positive yield moment per unit length getedh by a 


to the y-axis; and ke m= 


4 


bottom rebars placed parallel to -axis. It is to 


aa to a total ‘available y yield moment mk, and a coefficient of 
ax to ads si yd > olgaston A 
neg: yield 


The vector moments per unit length acting ig along a Positive or a a - yield 
line making an angle with the x-axis can be written as _ ne in 


b 


whet: follows, “extensive use will ‘be made of energy ‘input and 
absorption in portions of elliptical yield fans forming in orthotropic slabs. It 


is convenient to present the relevant expressions here. Fig. 2 shows the geometry © 
nd notation of elliptical fan. In polar equation of 


sin Tite. + cos’ 
The he angle R, and the ine, ei 1 be found fr 


a Then if point oO drops a virtual distance, 8, while the e elliptical p periphery 1 remains 
_ in the original plane ane forming the fan, the virtual energy input of a uniformly 7 
distributed load, per unit area acting” on the fan between $, and d, will 


R,, sin® 
sin 
or E,=— war [arc tan (q tan, ) arc tan tan ,)] 
The virtual energy during the yield hinge rotations in the fan 


given | and Jaeger Translated into the present notation itis 


tan” ry + q’ tan? o,+1 / 


| 
| 
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mer 
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) lare tan (q  tand >... 


a Fis. 3 shows “i collapse ‘mechanism caused by a a single concentrated load, 
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Portion of Elliptic Fan 
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a uniformly distributed load, w per unit area, in which P. If the 
] sc distributed load were relatively large, a different mechanism, based 
on the uniform load and the boundary conditions, might govern and would 
have to be investigated. Such mechanisms are dealt with exhaustively in the 
—— ——. symmetry in the mechanism of Fig. 3, using Eqs. 6 and 
7 with , = 0 and o, = 1/2 while noting that here r = b and q = a/b, — 
and w = *)/P, a = a/c, and energy 


E 
2a, +). a (a 


+p 
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FIG. 3 —Symmetrical Interior Concentrated Load Mechanism 


virtual energies and noting that from Eqs. | and 2 +1 
al and +ki= + is to write ~ 


sie to we 


+p 


mew) B + ka + — af 


In this unbounded mechanism ainsi that | the corner fans will: 


be. circular arcs if the slab is isotropic, i.e., p= 1. Then from Johansen’s (4) 
‘affinit theorem, g = = and E 10 “ be rewritten as . 


— 


It is obvious that the load carrying capacity of the slab will be a minimum 

for maximum values of 8, i.e., that the collapse mechanism must become as 
large as “‘possible.’’ That means the mechanism will spread until stopped by | 

a geometrical obstacle or a reverse, i.€., positive since ‘the peripheral yield 
line is negative, moment field, or both. For a slab supported along its entire 
periphery by strong beams or walls, such as that shown in Fig. 4, the collapse _ 
mechanism will spread until it becomes tangent to one or more of the slab 
boundaries. Which of the boundaries will be reached will depend on the geometry, a 
the loading, and on the layout of the reinforcement. This is in contrast to_ 
the equation if Fig. 3 is | the collapse 


FIG. 4.—Collapse Mechanisr in Slab with ti Supported Boundaries 
a uniformly loaded slab_ supported on an interior rectangular 


en this ‘equation is minimized for Fr. B is found to be self limiting, v with — 
Fig. 4 represents a slab supported on rigid boundaries, which have been labelled — 
_ but that the positive moment reinforcement will be constant throughout the ed 
- though it may be orthotropic. This corresponds to frequent design practice. 


"oad acting onthe slab 
‘in 
— 

the mechanism always having a definite size whose surface area is smaller 
than the area of slab required to produce P,i.e.,4ab+c*><P/w. 
_A, B, C, and D for easy reference. To generalize the analysis, it will be assumed — 
‘that the negative moment reinforcement crossing each boundary can be different, 


| 
a 
| 
| 
Bi: 
L 
_ 


coefficients | be: k k, over over : 
‘boundary C; k,, over boundary B; and k, > over boundary D. Note that the 
_ subscripts correspond to the subscripts of the dimensions in the direction of _ 
the r reinforcement. It is convenient to define 
= (k,. k’)/(k,, + ki); Bk, =k, we, 
a =a, /e, 0. = b, je: ai lh, Je, 
= Li Li /c. The loading will be assumed as 


= 2 zk, l+n, 1B, a, 4 


+0.) + (a, +a; XB, + 
This equation must be minimized with respect to’ B,, and in order 
to obtain the lowest upper bound on P/(m>k,). At the same time, the ro medal 
constraints require thata, =A,,a,=A,,B, andB, <A}. The minimization 
- _ is best done aumecically, either r iteratively on a computer or by trial and error — 
on a programmable calculator. It will quickly be found that the mechanism — 
—— will’b be tangent to one or more support boundaries and, depending on the magnitude 
of w, may reach all the boundaries. ‘The large: w, the more the mechanism 
tends to spread, and a/f = V yp only when the boundaries force this relationship re. 
literature (4,8,9,10) presents a different collapse for point 
irs loads from the one rere in Fig. 4 when the slab is ‘‘simply supported” on 
_ one or more boundaries. However, use of Eq. 13 with k, or k, = 0 for the 
: appropriate boundary and the region adjacent to it, will give a lower value 
for P. (As a quick illustration, for an isotropic slab with equal positive and 
negative yield moment m’ per unit length, loaded only by a point load near . 
a simply supported edge, the conventional mechanism | gives P = 11.42 m ‘ + ; 
Use of the yen 6 formulation, with the appropriate negative yield moment 
= 10. 73 


are enlikely to provide much top reinforcement to a simply 
: 7 supported slab boundary. However, if such reinforcement is provided and the | 
a so-called ‘‘simple’’ support can resist some torsion, then these factors can be 4 
incorporated into k, or k, (whichever is appropriate) and this will 


_ The mechanism of Fig. 4 was based on the assumption that the slab boundary 
supports would be sufficiently strong to prevent yield hinges from forming | = 
them. It therefore does not provide any information on the flexural strength — 
7 
sé _ required in the supporting beams. If required beam strengths are to be investigated — o 
or if a free slab boundary is to be considered, two additional geometrically _ 
§ admissible collapse mechanisms must t be analyzed. They are ‘shown in ee 
Looking first at the mechanism of figure 5, let W be the weight of the stem 
_ the spandrel beam per unit length; M, be its negative yield moment to the 


left of P; M, be its negative yield moment to the right of P; and M’ be its 


Positive yield moment. ‘stem width of will be assumed to be 


A 4945 
: | m’.) in addition to giving a lower upper bound on the load, the mechanism _— § 
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slab yield moments, noting that k, should i be taken equal to ‘zero, >, though its 
value will not affect P. Then proceeding as before, the equation of virtual 

We 


mik, 


This equation, , also, must be minimized with re respect to the geometry, pes ey 
the — imposed by the boundaries. If w has a realistic value and the 


Positive ‘Yield Line 


Negative YieldLine 


slab proportions a are within a ‘practical range, it is very ‘likely that P will be | 
£ minimum when the mechanism has reached two or more of the ite 


‘Thus a, = A,, a, = Ray and B, = A}. Of course B, will equal If the 
beam strengthe are known, it will be neseeeney to check P from Eq. 14 against 7 
P obtained from Eqs. 13 and 15. If the beam strengths are not known but 
+ P and the slab strengths are, Eq. 14 can be solved for the beam strengths. 
If W, M’, M,, and M, are set equal to zero, the solution for a slab with 
a free edge is obtained. If the slab is continuous across the beam with a mirror 


= 
| 
| 
— 4 | 
a 
a 
wat 
2 
| 


“may may be solved. using half the actual beam strengths and weight. 
‘The mechanism of Fig. 6 is an alternate to that of Fig. 5. It is likely to 
govern if \{ is relatively small and the beam is weak or the slab edge is free. 
A large value of Hw, and M2 will also f aver its formation. In band case the equation — 
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of the ellipses as 


sin’ + cos’ 
4 = arc sin Vai sin’ 6, + 


and the spread of the eittedsbe in the y direction can be obtained as a functi ion - 


ad in the x direction from the geometry of the triangles as 
(V q sin’ + cos’ 6, sin @,- 


toad tends to cause it to spread i in all directions toward the actual slab boundaries. 
_ Thus one or more of @,, a,, and B, (and therefore also B’) will be nour 


by the slab dimension. However, if one dimension is very different from the 
others, as when the concentrated load is located near boundary A or C, or 
if the distributed so is very small, the mechanism may not spread to all the 
_ It is not mathematically feasible to obtain the q, and , directly from B’/a, 
Fig. 7 has therefore been prepared from Eqs. 16 and 17 as a calculation aid. 


Eq. 15 can then be ‘minimized for q,, 72, and either iteratively or by 
_ trial and error, or a combination of the two. Setting values for a, and B’ from i 
the slab geometry, possible combinations of q and can be obtained from 
Fig. 7 and tried in Eq. 15 until the lowest value of P is obtained for the given 
conditions. The Process i is not as } cumbersome and - time consuming as i it may 
and a very acceptable value can 1 usually be obtained in half a dozen or so _— 
tries on a programmable calculator or minicomputer, hat 

£ If the Strength of the beam along edge B is ists to be Small, or if its” 


need be investigated. If be eam is known to be very strong, 


ky 
wae. 
he | 
; ven by q, and q,, the ratios of 
Tene tae elliptic arcs are determined ty 
 &§ the angles and ,. The angles 6, and 6, can be obtained from the geometry > 
— 
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'7.—Solution of Eqs. 16 and 


Por mechanism of Fig. . can control In between, ‘all three mechanisms may 


If the reinforcement is ccaiuged symmetrically and the concentrated load © 


is movable, the worst location for it will be at and A; = 0 with 
weak beam, or at = and MM = 3 for § a strong beam. “Symmetry of course» 


greatly ‘simplifies the various ‘equations and reduces the number of iterations | ‘= 


or trials required for minimization, 
Consider a reinforced concrete slab, 20 ft x 25 ft (6.1 m x 7.6 m) in area 
and 8 in. (0.2 m) thick. The | ‘short sides and one of the long sides will be 
e continuous ¢ over supporting girders, while the other long side is Giacontinnces 


and cast integrally with a spandrel beam which has a stem width of 12 in ; 


_ (0.3 m) and stem depth below the slab of 18 in. (0.5 m). The torsional strength 
of the spandrel will be ignored in these calculations, but it could be taken — 
into account if adequate tonional ductility were assured and continuity rein- 

wz will be assumed that the slab : may have to. resist a a “concentrated” service 
load of 40 kips (18 tons) with a footprint of 12 in. x 12 in. (0.3 m x 0.3. 
m), which can act at any point in the slab. The structure will also have to 
carry its self weight, and the effect of additional distributed load will also be 
; investigated. Applying | ACI Code a) load factors, the load ratios to be used 
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= a If the strength of the spandrel beam were to be critical, the worst location | 


the concentrated load would be just above ‘midspan of the spandrel. 
If the strength of the spandrel were such as to force failure to occur in the aS 
interior of the slab, the critical location of the concentrated load would be 
on one center line and nearer the simply supported edge than the fixed edge — 
a in the other direction. Orienting the example slab in conformity with Figs. 5 
and 6 so that the spandrel beam is edge B of the : slab, the ¢ two: critical locations 
| the concentrated load will then be A, = A, = 12, A, = = and 
‘These critical locations indicate that the x direction catiborcusitied should 
_ probably be symmetrical. Assuming grade 60 reinforcement, reasonable rein- 
: forcement ratios and symmetry in the layout of the x direction reinforcement, wn ~ 
the qs, preliminary ultimate (~ yield) moments will be made available: _ 
M, = M, = 8,100 in.-kips (915 kN-m); M’ = 5,400 in.-kips (610 kN-m); mk, 
= 30 in.-kips/in.) (3.4 kN-m/m); = mk’, = 20 in. -kips/in. (2. 
kN- -m/m); mk,, = 25 in.-kips/in. (2.8 kN- sake and mk, , will be assumed 
to be = 0, though it may have a small value. — — 
a The following ratios can be obtained from the whenlie 


ap 


effect of sys "symmetry will be to make = M, 
7 = a, = a, >, = >, = >, and q. = q2 = q. . These "designations will be 
Table 1 presents the results of s sample ‘calculations for | the problem. The top 
part of the table contains all the trials used in obtaining solutions for the example 
with each of the three possible collapse mechanisms. All calculations were carried 
out on a TI-58 programmable calculator to ensure that the solutions would 
_ be within the capability of any engineering office. . Without the advantage an 
_ symmetry, however, Eq. 15 would require use of a TI-59 or equivalent. | eat 
_ Looking at the top part of Table | and the solution process in more detail, 
the first mechanism analyzed was that of Fig. 6 and Eq. 15. On the assumption 
that the mechanism would spread as far as possible, an initial value of B, = 
5 = 19 was chosen. The first trial for d was = 0 (actually 0.0001 was 
used to avoid calculation difficulties) for which q= = a/B". This led t to an initial - 


Fig. 7, ‘with B’ ra = 20/12 67, gq = 0.59. For this trial, Eq. ‘15 gave 
= 3.01. A further trial, with g = 0.55 and = 0.24 led to P/m<k, = 3.06. 
 \iinslieding the possibility that the mechanism would not spread all the way _ 
to the far boundary, B, = 18 was investigated and led to almost identical results, 
as shown in the table. Setting 8, = 14 started an upward trend in the solutions. 
_ The numbers indicate that the solution surface is almost flat in the vicinity 
of the correct answer, making it unnecessary to perform many trials. — ie gates’ aan 
‘The process was repeated with the mechanism of Fig. 5 and Eq. 14. Again, q 
the very flat solution surface and the tendency of the mechanism to spread : 
made thre tres sufficient. This me mechanism produced a lower value of of 


| 


Set of physical parameters is highlighted with an beside it. 
_ Next, the load was moved to the interior of the slab. The mechanisms of 

‘Figs. 5 and 6 will give greater load capacities as the concentrated load moves — 
away from the edge. Therefore, only the mechanism of Fig. 4 and Eq. 13 — 
- was analyzed for this « case. Parameter k,. was kept = 0, but the weight and ¥ 
_ omen resistance of the spandrel beam were no longer part of the problem, — 
since the beam was now assumed to be sufficiently strong to force an interior 
Slab failure. Five trials with the mechanism spreading to all the boundaries, — 
g two additional ones in which the mechanism was curtailed in size, were | 


sufficient to establish P/ m&k, =  & 16 as the low value, and the critical location _ 


= 12, A; = 7.5, and 1:5. 
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At this Point, th the preliminary design can checked. For t the slab itself, 
with P in the interior, m=k, = P, /2. 16 = 68/2.16 = = = 31. 5 in. -kips/in. 6.6 
af kN). The assumed mk, = 95 in.-kips/in. (10.7 kN-m/m), and the assumed 
reinforcement could be reduced by a factor of a about three. 
amuch heavierload could be carried. 
7 - Toobtain some “‘feel’’ for this type of problem, some of the physical aide 
- were varied. The minimum results for each mechanism are shown in the lower 


half of Table 1. First, w, i.e., w, was set = 0, a common assumption when — 
7 dealing with large concentrated loads . The result is a 28% increase in the 
concentrated load capacity. Doubling w on the other hand, reduced P by 22%. ct 
ieoningy, when w =¢ the two edge load mechanisms gave the same value 


re “ ry) evious one, and will therefore govern for the load at the edge in 7 
Ps 
of jimi ia | 


: 
for P but with two different size mechanisms, neither of which spread to the © 
- far boundary. When | w was doubled, so that w = 0.004, the mechanisms behaved — 
as with w = 0.002 except that B, = 18 gave ‘somewhat larger values of P/m<k, 
- Reduction in beam strength to slightly less than one half its original value 
- caused the edge load mechanisms to become critical. For w = 0.002 both gave 
the same concentrated load capacity, for » = 0 the mechanism of Fig. 6 governed 
and for w= 0.004 that of Fig. 5 gave the lowest strength. In most of these 
= the mechanisms did not extend to boundary D, but did reach boundaries 
, 5 Complete removal of the edge beam had similar effects. It is interesting to 
‘note, however, that the slab as originally designed could still safely carry . 7 
- "concentrated load approximately | 1 3/4 times the assumed P, plus its self If weight, 


To investigate the accuracy of the “superposition method compared to to the 


mparison of Eq. 13 and Superposition (Homogeneous, Isotropic 


Lal 


Position Eq. 13 
(6) 


Li the + L‘, with fixed boundaries and equal top and 
BR reinforcement, can easily be analyzed. Taking corner levers into account, 
the uniformly distributed load carrying capacity of the slab wie £2 0) will 


(4) = 11 When = | and assuming za, = 
(12 and w = 0, Eq. 13 gives P = 3. 31 m<=k,. Using superposition, the total : 
“required =(wL’)/11+ P/3. 31. Rewriting i in the notation used previously, 
expressi 


and using the aforementioned dimensionless ratios, this superposition 


‘Results use of this can be checked against puny 
obtained by use of Eq: 13 for various magnitudes o of w. | 
_ The upper half of Table 2 shows some results from both superposition and 
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‘CONCENTRATED LOADS. 


works very well ‘this case, in which the tend toward. 
_ The lower half of Table 2 shows results from both superposition and Eq. 
13 when the concentrated load is moved along the diagonal of the square plate — 
oh, = hj = 4 and A, = AS = 20. With w = 0, minimization of Eq. = 
gives P/ (m3k,) = 3. 625 and ‘therefore the superposition equation, Eq. 18a, 


The lower half of Table 2 shows that the are now no longer congruent, 
: and that for realistic values of w superposition will lead to a gross underestimate re 
: of the strength of the slab. It is interesting to note, also, the effect of the nd 
relative magnitude of the distributed load on the spread of the collapse mechanism 


as shown in Table 2 by the values of a, and B, for which P/(mzk) became 


DDITIONAL UBSERVATIONS 


‘Other calculations which have been carried out that the 

solution to Eq. 15 is almost totally independent of the values of k, / ky and »¢ al 
4 k.,/k Multiplying or dividing these parameters by factors of three produces — A 
7 less than 0.1% change in results. Changing k,_/k’, similarly produces approxi- _ 


; mately 10% changes in results. This is in agreement with previous findings — 
> 


the slab are much shorter than edges B and D, the mechanism will invariably 
spread to edge D, but may not reachedgesAandC. 22 
“a 2. Reducing the size of the footprint of the concentrated load will reduce | 
its permissible magnitude. However, this effect is also relatively small for practical | 
«a cases, since the footprint must be kept sufficiently large to prevent shear punching. ; 
— . more than one edge of the slab is to be simply supported on collapsible 
ea, the mechanisms of Figs. 5 and 6 must be applied to each edge in turn. © 
ih addition, an elementary overall collapse mechanism (4) with yield lines over 
and between column lines must also be investigated. hai sci 
_ 4. If the collapsible beam is an interior one, with the slab continuous over 4 
it, the mechanism of Fig. 5 will be critical for the concentrated load acting 
on top of the beam. For calculation purposes it will then be useful to = 
‘f with half the beam and half the concentrated load in conjunction with the slab — 


| for uniformly loaded slabs 3). It was also found that if edges A and C of 


B, = d' # 0. The larger B, the greater the likelihood that the 1 mechanism @ 

& will give lower carrying capacity than that of Fig. an 
7 6. Only rectangular footprints have been considered here for the concentrated — 
2 Previous experience with the flexural strength of uniformly loaded slabs _ 

_ supported on columns of various shapes (2) has shown that isotropic arrangement _ 
the reinforcement will lead to lower (on the order of 1%-3%) flexural 


| 
for orthotropically reinforced slabs, where appropriately oriented 


— ‘columns sometimes gave slightly lower strengths. The differences were always 
too small, however, to be of practical significance, and for normal design purposes" 


: it will be quite to substitute footprints area for 


Multiple Joads acting on as well 
as mop concentrated loads acting on slabs supported directly on columns = 


single load mechanism may be expected to form, and a reasonable approximation 
“4 to the slab strength may be obtained by smearing the smaller loads into ‘the 2 
- uniformly distributed load. In the flat plate problem, an intcraction on (between — 
7 load and column) mechanism has been identified and analyzed for some simple 
& _ cases. It leads to the lowest collapse loads for a fairly narrow range of geometries | 

and load ratios, with overall slab collapse mechanisms or pure flexural column 
punching t mechanisms controlling more frequently. 


* Yield Line Theory h has been used to derive expressions for the flexural capacity 
of orthotropically reinforced concrete slabs subjected simultaneously _ to a 
concentrated and a uniformly distributed load. The weight and flexural strength — 

_ of a boundary beam are included in the expressions. Example calculations have — 7 d 
shown the effects of the location at which the concentrated load is placed, 7 4 
of the magnitude of the distributed load, and of the strength of the boundary — 
beams. It was found that the collapse mechanisms have a tendency to . spread 7 
until stopped by one or more boundaries, and that the tendency becomes stronger . 

with increase in the distributed load. It was also found that the solutions are 
relatively insensitive to the negative to positive yield moment ratios of the 
slab and to the dimensions of the collapse mechanisms. An exploration of the 
possible use of superposition of solutions revealed that the results could be a 
quite accurate or ‘grossly inaccurate, depending on how nearly goemeutenlly 
4 
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ments OF the load Cause little difficulty: The loads closely Spaced 
/ _ and approximately equal in magnitude (as the load carrying wheels of a forklift) _ 
: they may conveniently be treated as a single load spread over the actual region : 
occunyv: if one loads is uch larger rs 
y 
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The following symbols are used in this paper: thee u 
label for one boundary are 
x direction dimensions of ‘collapse n 
label for boundary -slab—when slab has one 
ere ony in boundary supported on collapsible beam, Bis that bound- 
= y direction dimensions of collapse mechanism; Paced me 
label for one boundary of slab; 
x direction dimension of area covered by concentrated 


- D,,D,,D, = symbols used to represent other expressions in Eq. 15; 


, = virtual energy input by loads into collapse mechanism; — 
virtual energy absorbed due to rotation of yield lines lines 


collapse mec mechanism; at tient tapes. 


ij boundaries of yield fan; 
"negative yield moment per unit by top 
tension reinforcing bars placed parallel to x-axis; 
_— yield moment per unit length provided by bottom 
tension n rebars placed parallel t to x- ‘axis; 


tesion rebars placed parallel to y-axis; 
yn positive yield moment per unit length provided by bottom ia 
tension rebars placed parallelto y-axis; 
slab dimensions locating position of concentrated load 
= seb dimensions locating position of concentrated load - 


negative yield moments of spandrel beam; (im 
positive yield moment of spandrel beam; 
an m aa: unit it yield n moment per u unit length of slab; iat rabtcadoomy, 


“magnitude of concentrated load; 8 Septam | 


>| vee applied concentrated load; 
Tatios of x-axis to y-axis length of elliptical yield fans; 


i 
i} 
q 
7 
9 k 
| 
| 
Ry fadius of Curvature of ellipse making angle @ with x-axis, 


= per unit length of collapsible beam; 
magnitude of uniformly distributed load per ‘unit 2 area 
x,y coordinate axes; 


= = a/c, a,/c, a,/c, respectively; 
ww = virtual displacement of designated point in collapse mech- 
6,8, angles between radii and circumferences of elliptical yiel 
— . dee he = y direction dimension of the area covered by concentrated 
= Ly L,/¢, respectively; 
angles yield fan radii make with x 
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Recent experience that dusting a major earthquake stiff buildings 
having good ductility in ‘the postelastic range behave better than others (11). < 7 
_ For this reason, there is a strong renewed interest in reinforced concrete shear 

walls and braced steel frames for seismic design. Such structural systems are 

inherently stiffer than the widely used moment resisting frames and tend to 
reduce costly nonstructural damage during a strong earthquake. Moreover, often | 
it is economically impractical to use moment resisting frames ot Oe a narrow 4 
_ The key element in determining the behavior of a conventionally braced steel 5 
‘thai in resisting seismic forces is the brace itself (6). Based on experiments | “ 
with small specimens, the qualitative features of this behavior under severe — 
- cyclic loadings simulating earthquake action are now well understood (5,9,14). | 


4 However, quantitative information o on \ braces of sizes comparable to wren used 


‘t framing systems, such a as eccentric bracing (10), which « can be ‘used a as s alternative 


r sizes employed in practice are reported in this paper. The geometries of some 
= of the Specimens were so selected | that they also simulate Some frequently used ; 


; of data on n the inelastic behavior of struts s under severe reversing ‘loads. An 
approach for predicting analytically the deteriorating capacity of struts under 
extreme load reversals is suggested. Some practical implications resulting from 


this work are summarized at the end of the - 


_ The individual specimens were selected from standard structural steel shapes = 

B primarily on the basis of two criteria: (1) That the slenderness ratios of the 4 
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_ test specimens be appropriate to those used in practice; and (2) that the member 7 
"shapes and proportions be such as to represent typical brace and strut members o- 

in current use. The maximum size of the members was set by the capacity _ 

of the available double-acting hydraulic cylinder. 

_ Since the effective slenderness ratio of a compression member has been shown ; 

to be the single most important parameter in 1 determining its hysteretic behavior _ 
a, 8,9,16), care was taken “that the chosen Kars allow the Specimens to be 
ile with one another, as well as with members used in practice. A common $i 


slenderness ratio of 80 was used for specimens within each structural group 
to allow for comparison of results due to variation in cross-sectional shape. 
— In addition, slenderness ratio ratios of 40, very close to the range of plastic action, 
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FIG. 1.—Typical Details of eee (a) Pin-Ended; and (b) Pin-Fixed Ende 


= of 120, very near the elastic buckling range, were assigned | to several 


“3 To assess the effect of end restraint on brace behavior, two bounds on the 

possible end conditions were considered; both ends pinned, and one end pinned - 

and the other end fixed. It should be recognized that in actual structures, something | 

other than these idealized end restraints would be expected, depending on the 
frame configuration and loading, member slenderness, initial camber, and joint 

4 

- Included within the 24 selected s specimens were six different cross-sectional — 

_ shapes: wide flanges, thin and thick-walled square tubes, thin and thick-walled a 

soune, pipes, structural tees, and fabricated double angles, and double Gams. 
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3/8 i in. (9.5 mm) back | to back of angles and channels; ; for double angles the shorter _ 
legs areturned out. 


The built-up specimens were of special interest to see if the two elements would 

_ effectively act as a single member during extreme inelastic cyclic loading. Ss 

ms: Typical details of specimens are shown in Fig. 1. In all cases, the test struts _ = 
_ were welded to | 3/4 in.-2 1/4 in. (44 mm-S7 mm) thick end plates with full 
ms penetration welds. The specimens so prepared were attached to end fixtures 
with large high-strength bolts. The specimen length from pin to pin, or from 


the fixed end t to the pin, end clevises containing roller 


used in brace applications. The specimens were tested with different combinations 
of end conditions, K//r ratios, and cross-sectional shapes as previously men- ‘ 
tioned. A complete list of the test specimens with pertinent information is ai ae 


Test Set- Ur 


: commercially available sections of sizes and shapes similar to those typically 


2. Depending on the desired end ‘conditions, the 1 right- -hand ‘end shown in n Fig. 
2 was either pinned [Fig. l(a)], or rigidly attached to a foot frame or to a 
; steel beam prestressed to a concrete reaction block [Figs. 1(b) and 2]. The | 
other end of the strut was attached to a double-acting hydraulic jack that was 
‘pin- ~connected to another concrete reaction block. The jack-specimen 
a was made via a headpin assembly that was restrained horizontally by a ie 
arm and vertically by bolts in slotted holes. This arrangement allowed the primary 2 
buckling of the specimens to take place in a horizontal plane with negligi . 
small lateral movement of the headpin. a ee 


A load cell was inserted between the  spoemnen and hydraulic jack | to aie 


g used to measure the axial displacement between the ends of 1 a specimen, 
and linear potentiometers were used to determine the maximum horizontal = 
vertical displacements. The deflected shape of the specimens me the displacement 
of the head unit at selected points were recorded using photogrammetric 
‘techniques. Strain gages were placed at various locations along the — 
to monitor material yield histories (Fig. 1). Gl 

All specimens were subjected to quasistatically applied cycles. of reversed 


axial displacements. These cycles generally resulted in compressive loads 
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SEVERE CYCLIC LOADINGS 
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¥ FIG. 3. —Strut 2 22 in Test ‘Bay tara 
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followed by tensile ‘loads to cause yield in the 
os specimen. The displacement history imposed on one of the fixed-pinned specimens — 
q 


is shown in Fig. 3. While other specimens received similar displacement patterns, Ee. 
they were not identical. The 4 in. x 4: in. x 1/2 in. (100 mm x 100 mm 


x 13 mm) square tubes, for example, could not be caused to yield i in tension 
as the yield strength for this shape was 82 ksi (565 — requiring the development % 


° 


AXIAL DISPLACEMENT (IN) 
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. 4—Typical Sequence of Axial (1 in. = 4 mm) 


aa 8. The P-8c curves trace out the hysteretic loops. The areas enclosed al 
_ by such curves is a measure of hysteretic damping. Examples of ae 
curves for several members with different cross sections are given in Figs. yi 
. 6, and 7. A large number of additional examples may be found in Ref. 
2. The most striking feature of the these results is the indication that once 
a strut had buckled, during the subsequent cycles the same capacity in -ggrseerall 
— be reached. This deterioration of the critical buckling force P., is closely — 
related to the slenderness ratio of a member. Typically, the hysteretic loops” t 
for struts with large values of K//r showed a more rapid deterioration in - 
strengths than those with small values of [compare Figs. 


4 5(a) and 5(b)). . The Fatio of a member’ s in to 


OBB 
— 
| ecimen was exhausted, | 

dit, 


_ While the P-8 curves are of prime importance for the inelastic planar analysis — 
= _of frames, one should note the large lateral displacements which are associated 
le with a strut’s axial shortening, e.g., from Figs. 6(a) and 6(b) for Strat 4 it 

can | be seen that ‘corresponding to an axial displacement 5 of epgronimetely. 


= 


fe) 


ij 


FIG. 5.—Hysteretic Diagrams | of of Axial F Force Versus Axial Displacement for Struts Bs 


& in. fac) mm), the midspan lateral displacement A is nearly 9 in. (230 mm). 
Another important observation which « can be made from the P-A curves is the — 


. retains some midspan deflection. For the s same 1e example, note e that at zero ae 
wv - load following yielding in tension at the initiation of cycle 6 the residual deflection wt 
at midspan is 1/4 in. (6 — Thus, this specimen begins the 


i 
— 


compressive cycle as already a slightly curved member can 

s In order to explain the behavior of a strut subjected to inelastic cyclic loading, a7 
experimental data on the history dependence of the mechanical properties of —s 
a material are also required. This information is provided by the hysteretic 
curves for a material. A set of such curves for the W 6 x 20 specimens is” i> 


FIG. Fo Axial Versus. Axial and (b) Axial For Force Versus 


shown: in Fig. 8. Here, it is to note that upon the 
_ tangent modulus E, is considerably lower than the initial elastic modulus E. a 
While the number of cycles plays some role in this phenomenon, the — 
factor appears to be the amount of cumulative inelastic straining experienced 


— 


Ss the material before the beginning of a new cycle. The amount of the cumulative 


7 
= 


a cycle in question, e.g., for the second compress e 
equal to the strain from a to b plus the strain from 5 to c. In using this parameter, __ 


2 

a 


eretic Diagrams 
effect which ha: 
and includes cases when a strut is initially caused to yield in tension. 2 20 tt 
—— hysteretic loops shown in Fig. 8 are for a material which initially exhibited — 
a well defined yield point. This was typical for the steels in wide flanges, — 
_ double-angles, and double-channel specimens for which the yield strengths varied md 
Bis from 40 ksi-50 ksi (275 MPa-345 MPa). Most of the steels for tubes, however, _ 


# for the same level of stress the tangent moduli significantly decrease with an 
— : increase in the amount of cumulative strain, e.g., the corresponding tangent _ 
= § moduli in compressive cycle 6 for which Le ,= 0.054, are smaller than those | 
3 where Se, = 0.009. Bauschinger 


0.2% strain offset or by noting ‘the a attained ‘yield platens of a a strut 


during the tensile phase of a cyclic test. On these bases, the tube yield strengths 
varied from a low of 24 ksi (165 MPa) for Strut 21 to a high of 82 ksi (565 
_ MPa) for Struts | ‘B and et _A typical stress-strain diagram from a coupon test 


000 0008 


ain Diagram for W 6 x 
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STO. PIPE 31/2" 022 


STRAIN (IN./IN. OR mm /mm.) 


FIG. 9 —Monotonic Stress- Strain Diagram | for 31/2 in. (89° am) Pipe Coupon 


P for one of the tubes is shown in Fig. 9. As will be pointed out later, for struts 


made of such a materiel, in common with all steel struts with a previous inelastic = . 
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«FIG, 8 —Hysteretic 
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« ‘experimental initial strengths « determined « on ‘the. aforementioned 


’ = for the struts are also recorded in the table. Some representative cam 


pee | 
Square inch in ki i estimate 


0s 


calc 
= 


“ Average yield from coupon tests. 


7 “First yield in strut test. 
“Initial max A =0. 094 in. 


Initial max A = 0.16 
"Tangent Modulus theory. bea: 


Note: ksi = 6.89 MPa; | kip = 4.45 kN; and1lin.25.4mm. 


of struts at the end of tests are shown in Fig. 10. Fig. 11 shows two examples _ 
. of local buckling occurring in double-angle struts. This buckling develops in | 

the region of a plastic hinge at a stitch fastener in Strut .¥ and between en two 


— 
Strut 
number 
| “a 4 
T 
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such fasteners in Strut 9. In catia’ instance (not shown), during some severe 
cyclic excursions complete s at one of the stitches of Strut 20 had 


Initiat AnD Cycuic Bexavior oF Strurs 

Initial Buckling Loads.—In relating results” with 

= comparison of the carrying capacity of a strut to existing codes is particularly 
important. For the range of the slenderness ratios used in these experiments, 
the initial buckling capacity of a strut based on Eq. 1.5-1 of the AISC Specifications — 
(7) without the factor of safety y is appropriate. Using this approach, comparisons — 
with the experimental results are given in Table 2. Since in practice the yield = 
‘strength for the steels used in these experiments would usually be assumed © 
to be 36 ksi (250 MPa), a comparison based on this value of the stress is 

included i in n the table. The strut ut buckling capacities bases on m the ane formula > 


= In general, the experimentally « determined initial buckling leads for the struts 


are in reasonably good agreement with the predicted ones based either ono, _ 
36 ksi, or on the values of o, found from tests. However, improved correlations - 


: between theory and experiments can be found 1 by recognizing the nonlinearity _ 
the stress- strain -Telations for some ne steels, or by taking i into account 


is appropriate for resolving the first problem, and a well established procedure 
5) for determining the capacity of eccentrically loaded elasto-plastic columns ~ 
can be used to advantage for the second problem. As can be seen from the — 


OS 
=<. — 
: 
— 
4 


~ 


values given in the last column of Table 2, consistently better estimates of 

buckling loads can be found using these theoretically based refinements. 4 
of applying them is more considered in the section on ‘Analytical 

Predictions of Cyclic Buckling Loads. Wes 

_ Normalized Hysteretic Curves and Envelopes. —By dividing | the axial force wet 
in a hysteretic P-8 diagram for a strut by its “capacity P, in tensile yield, and 
; the axial displacement 8 by the displacement 8, at yield, one obtains a normalized ' 

; hysteretic diagram. Such diagrams exhibit the experimental results in a meaningful 
_ manner by eliminating the effects of variations in material property, cross-sec- { 
- tional area, and specimen length. Further, because < of an infinite variety of 
- cyclic patterns that may be applied to a strut, it is convenient to make use 

_of envelopes for a family of normalized hysteretic loops. In this analysis (ue 
; comparisons of cyclic behavior among struts is made using such envelopes. _ 

As an example, consider the envelope for the normalized hysteretic loops — 
for two identical Struts 3 and 4 shown in Fig. ». 12(a). As can be seen, the 
= of the envelope is essentially the same for the two struts which had . ; 

_ different loading histories—a fact of considerable importance in studying random - 
loadings occuring during severe earthquakes. This is true, however, only because 
_ the specimens experienced similar loading patterns. Both specimens were initially . 


com compressed, and each initially ae eeeaentinns its maximum buckling load. By contrast, ae 


| 
q 
— 
| 
| 


an strut, which initielly. was caused to ‘yield in tension, reached | only 
about 75% of this capacity [Fig. 12(b)]. This is directly attributable to the | 
j Bauschinger effect as a result of which the stress-strain diagram in compression 
is significantly rounded, reducing the elastic range. The maximum compressive 
load of 152 kips (676 kN) with the associated axial displacement for Strut 5 
- falls o outside the envelope of Fig. 12(a). This is due to the fact that the envelope — 


5 


FIG. 12. Normalized Envelope | for 3 and 4 4; and (b) 

_ ass there was established from the tests which began with the application — : 4 
_ of a compressive cycle. At the beginning of a second hysteretic loop [such | 
as at point c in Fig. 12(b)], the specimen’s cumulative plastic strain was equal 
to the distances ab plus be along the abcissa. By contrast, at point d, the 
a beginning ng of the first compressive cycle. for Species 5, the cumulative plastic 


MALIZED AXIAL LOAD P/Py 


q 
q 
a 


strain is given smaller. ‘distance aa ad. differences in the 

[ - history of loading are taken into account, the maximum first cycle compression S 
_ load for Strut 5 can be made to lie on the envelope for Struts 3 and 4. The 

“< _ deterioration of the buckling capacity of a a strut it due 1 to © previous plastic working 


An approximate 
"procedure for predicting the ‘buckling. capacity ty of struts subjected to random 
loadings is considered i in the section on Prediction of Cyclic 


ffects of Ratio, Boundary ond Cross-Sectional 
on 1 Hysteretic Behavior. —Thee effect fect of these ae on the hysteretic behavior 


(W6x20 STRUT 3) 


“ 


“<a 
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PINNED-PINNED 


FIMED-PINNED 26x 31/2 23/8 (STRUTS 


(@342 20276 stRuT 24) (STRUT IB) 


- Different Slenderness Ratios; (b) Struts with Different Boundary Conditions; and 
Struts with Different Cross-Sectional Shapes 

_of struts can be conveniently examined using the hysteretic envelopes displayed 

in Fig. 13. The dominant influence of the slenderness ratio on the hysteretic — F 
_ behavior of struts is clearly brought out in Fig. 13(a). _ For the slender er Strut 
11 the diagram is strongly biased, and the ratio between the capacity of the — 
strut in compression to that in tension is significantly smaller than it is for D 
_ the stockier struts. The initial buckling capacities of Struts 2 and 3 were nearly 

_ alike, however, upon repeated load reversals the maximum compressive loads 
for slender ut 3 rapidly deteriorated. In the limit, for a strut with 


| 
| 
* 
| 
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Fig. the normalized hysteretic — for a pin- ~ended specimen 


effective length concept i in n the inelastic range of material behavior for ‘cyclically — ‘ 
loaded struts. Photogrammetric studies reported in Ref. 2 indicate that during 
; inelastic buckling the curvature tends to concentrate in the regions of plastic 


hinges, but points of inflection do not migrate from their locations during the 


‘The effect of cross-sectional shape on the hysteretic behavior of struts can 
be compared with the aid of the normalized hysteretic envelopes such as shown 
in Fig. 13(c). This comparison together with other cases given in Ref. 2 indicates 
the following order in progressively poorer performance: tubes, wide flanges, __ 
tees, double channels, ane double angen. The characteristic properties of the 


to ‘fool buckling, lateral- torsional buckling, and, where applicable, ‘to 
_ With the exception of the struts date from thick-walled pipes, all specimens 
showed local buckling in advanced stages of cyclic load application. ‘Such local 
buckling occurred at extremely large lateral in the” Tegions of 
Lateral-torsional buckling was in some tee and double- -angle 
members. This occurs in singly symmetric sections if because of their geometric 
Proportions (3) they buckle in the direction perpendicular to their axes of ti 
symmetry. In such cases, flexural and lateral-torsional buckling take place 
simultaneously ‘causing ; a lower critical load than that which would develop — 
‘in pure flexural buckling. This tends to contribute to the somewhat poorer — 
_ pemteemance of tees and double angles in comparison with tubes and wide flanges. - 
_ The most significant effect responsible for the poorest performance of the 
- double. angle struts, and which is likely to apply to other built-up members, 
is local buckling of the individual members between stitches. Struts 8 8, 9, and = 
20 which had the propensity to buckle at right angle to their cross- -sectional 
axes of symmetry, and thereby also developing some torsion buckling, were 
particularly poor. The double angle Strut 8 had only one stitch at its midlength; ‘- 7 
= Strut 9 had two located at third points. In conformity with AISC | 


Tequirements, the slenderness ratio of individual between the fillers 


- the initial application of the compressive load, both of ‘Giomet struts behaved - 
well and their buckling capacities at first buckling load were good (see Table 
2). However, as severe cyclic load excursions were applied, the angles ode 
to buckle locally in the region of plastic hinges. An example of this behavior 
in advanced stages of loading for ‘Strut 8 is shown in Fig. 11(a), and for Strut 
9 in Fig. 11(b). From these figures, it is apparent that as cyclic load oo, 
proceeded, the flexural straining was concentrated in the middle of a specimen, 
and the webs tended to approach each other. This behavior was particularly - 
pronounced in Strut 9. Cyclic plastic working of a strut in a hinge region softens 

_ the material due to the Bauschinger effect contributing further to its deterioration. 
A complete failure of a stitch ‘Strut The Spacers” for this 


| 
j 
> 
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strut were located at venep points along the length and were - 1/2 in. (64 mm) ; 
_ wide attached with four 5/16 in. (8 mm) fillet welds 2 1/2 in. (64mm) long. 
ANatyTicat PREDICTIONS OF Cycuc Buckune 


Inthe analysis of diagonally braced frames for seismic excitations, the inelastic — 
_ buckling behavior of struts must be mathematically modeled. A number of 
Pa proposals for accomplishing this have been made (1,5,6,9,14), and the subject 
- continues to be an active field of current research. Here, the apne 
evidence on the behavicr of struts subjected to inelastic cyclic buckling is 
_ examined and analytical procedures for predicting approximate buckling capacities 
struts for the first few cycles are developed. 
_ As has been noted earlier,there are two main causes contributing to the often 
dramatic decrease in a column capacity for inelastic cyclic loadings. These are - 
4 the Bauschinger effect, exhibited by the steel subjected to inelastic load reversals, of 
and the effect due to the residual curvature of a specimen resulting from plastic 


during previous cycles. Each one of these effects can be 


4 


estimated by | using ‘the tangent ‘modulus i in the generalized Euler formula. Using 
this approach, the curve shown in Fig. 9 for Strut 24 was used to determine > . 
: its buckling load (see the last column in Table 2) with good results. The same 


extending the tangent or the 1 reduced (double) modulus approach to ) cyclic loading. 
To do so, however, requires that some approximations be introduced. 
_ Since the behavior at a plastic hinge is dominant in affecting the overall 
performance of a strut, it can be assumed (16) that the strains at the centroidal 
fiber of a hinge are decisive on the buckling behavior of a strut as a whole. 
> It is recognized that this is a drastic simplification of the problem. Significantly — 
_ different strain histories do occur elsewhere. Nevertheless, because of the key _ 
_ importance of plastic hinges on the buckling behavior of a strut, the proposed 
_ simplifying assumption seems reasonable. Next, each curve in a stress-strain 
_ diagram of a cyclic coupon test can be considered | ‘to represent a monotonic — 
_ test on a material with a previous inelastic history. This history dependence 
can be conveniently approximated and defined by the absolute cumulative plastic 
= - strain at the beginning of a loading cycle, e.g., consider the tensile loading 
curve for cycle 2 from a coupon test for the W 6 x 20 section shown = 
_ Fig. 8 reproduced in Fig. 14. At the beginning of this loading cycle, the absolute _ 
cumulative plastic strain Le, at zero stress is 0.0055 in. /in. (meter/ meter). 
This is the sum of strains along the abscissa in Fig. 8 from a to b to c and 


then to d. Similar curves can be isolated and identified with different amounts 
of cumulative plastic strain for the other loading branches of the coupon test. r 
_ After establishing a stress-strain curve identified with a particular absolute 
cumulative plastic strain, | as h has been done in Pig 4, the tangent moduli odened 


ot 
| 
by 
e | Kind Of stress-strain Curves are Clearly exhibited by steel in the post-yield range a | 
4 


in which E, and o., = the matching values of these oietinies. found a 


a stress-strain diagram such as Fig. 14. The calculated slenderness ratios K¢]r 
can then be plotted versus the critical buckling g stress ily One of the resulting 
curves obtained in this manner corresponding to the data given in Fig. 14 for . 
“a The family of curves in Fig. 15 identified with different amounts of cumulative © 
bs plastic strain | has been rd orion in the aforementioned manner. However, in 


=0 was fi found using the. AISC 


a 


0.002 


5 = formulas with no factor of ‘safety. The corresponding theoretical curve 


based on the simplifying assumption of ideal elastic-plastic behavior is known _ 
to be inaccurate in the relevant ran 


“ nge (13), and the two available experimental — 

points were considered to be insufficient to define the “required curve. With 

this data, for a selected slenderness ratio of a strut, its capacity for a sive 

_ cumulative plastic strain Le, # 0 divided by the capacity at re, = 0 gives 

reduction factor R » accounting for the Bauschinger effect. 


‘Since the Struts in this series of experiments tended to a residual 


appropriate t to saudi the reduced modulus theory rather than the ‘aforementioned — 
> 
tangent modulus daca for determining the cyclic buckling loads. The rationale — 


to the corresponding formation, one 
| 4 can make use of the ge ling slenderness 
ratios K4/r. The Euler formula in the appropriate form for this purpose reads 
| 


i process of the fibers on the convex side of a strut is likely to occur | 
earlier than it does in an initially straight member. 
4 By limiting the application of the reduced modulus approach to the wide 
flange sections used in these experiments, the procedure is very direct. Since 
all of these members buckled around their X-axes, on neglecting the contribution 
of the web, the reduced modulus E, can be taken as one for a rectangular 
section given by (13) 


Whee on ‘exchanging E, in Eq. procedure for establishing the 


tie buckling curves as has been done in Fig. 15 can be repeated, and ae 
. reduction factors based on the reduced modulus — can be determined. oa 
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giving the arrived at on the two 
bases to account for the Bauschinger effect as a function of the cumulative 
plastic strain are shown in Fig. 16. As to be expected, the tangent modulus 
approach indicates a larger reduction in the capacity of a strut than that predicted : 
_ Reduction Factor Due to Specimen Curvature.—The second major cause for — 
the decrease in column capacity during cyclic loading is that after an initial Co 
inelastic buckling cycle a specimen develops a residual curvature which generally — 5 
is not removed by the subsequent tensile yielding. Therefore, an inelastically — 
: cycled m member must be treated ir in th the analysis as having an initial curvature 


available 


tire 


“CUMULATIVE INELASTIC STRAIN, (IN/IN) 


FIG. 16.—Column Buckling Reduction Factors Due to Bauschinger Effect on Two 
‘Bases: Tangent Modulus and Reduced Modulus See 


@/S ECCENTRICITY RATIO 
cng es ad he 
FIG. 17. —Column Buckling Curves for Eccentric Loadings (15) 


fe RQQQ | 


in (13,1 15 for loaded elasto- plastic columns. 
_ the solutions obtained by Westergaard and Osgod based on von Karman’s concept > 
for inelastic buckling of eccentrically loaded columns are utilized. Some of 


the column buckling curves obtained by ‘them for eccentrically loaded struts — 
(15) are reproduced in Fig. 17. In this figure, e denotes an —— of the 


ratio of e to s defines the eccentricity — 
=) As an approximation to the problem being considered here, the experimentally _ 
Cone maximum effective lateral deflection A of a strut at the beginning | 


a 


fe) 
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| | tine 


er FIG. 18. —Column Buckling Reduction Factor Due to Eccentricity Effect — a 


3 ofa compression cycle can be taken as e. Adopting this approximation , the 
curves of Fig. 17 provide the necessary information for obtaining graphs for | 
the reduction factor R, as a function of the eccentricity ratio e/s = A/s. . 
Such graphs can be constructed in the following manner. Fora selected slenderness © 
_— Tatio such as 80, the critical values of the stresses on the e/s curves in Fig. J 

_ 17 are normalized using the capacity 0 of a straight column (e/s = 0). A continuous — 
curve connecting these points gives an R, plot for the selected column slender- 

ness ratio as a function of the eccentricity ratio. Two graphs of this kind are 
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Variable 


= 


ksi = 16.3 


x x : 70 | 16. 67 
1.00 | 1.00 | 100 | 1. 0.88 
| 


for E, an 


Note: | ksi = 6.89 MPa. 


er shown in Fig. | 18, Note how rapidly t —_ capacity | of a ctu decreases with 
an n increasing eccentricity ratio e/s = A/s. 
a Comparison of Analytical and Experimental Results. —By makighying on an init 

_ buckling load for a straight virgin column by the appropriate reduction factors 
'R, and R,, such as given in Figs. 16 and 18, a buckling load for the new 
condition is ‘feund. Some such calculations are summarized in Table 3, where 
the results for the first three consecutive cycles for Struts 2, 3, 4, 5, and 19 
are given. For the purposes of illustration, it was assumed that for Strut 2 _ 
the R, factors could be based on the graph of Fig. 16, although no cyclic 
= tests were made for the material of this strut. ed 


The calculated re stresses for the first cycle of the struts listed in 


a reduction factor R, was cutie. Further, since Strut 5 was initially perrae ~ 
ene force causing the member to yield, a reduction factor R , was employed. 

In these two cases as well as in all others, the experimental data were used 

to find A and Xe,. By applying in a similar manner the required reduction 

- factors to the second and third cycles, the corresponding estimates of the buckling _ 
_ stresses were found. As can be seen from the table, with the use of these _ 
factors the estimated buckling loads are in reasonably good agreement with 


q 
CycleNumber 
a 
0.48 
| 14.8. 
0.65 
0.80 
9.73 4 
| 


VERE CYCLI ADI J 
SEVERE CY CLIC LOADI DINGS 


at First Three Consecutive Cycles 
ned: 


: 


1 
aay aay 


0.97 84] 1.00 
9 
37. 32.6 
| 101 | 1.07 


- i vest emsing that the use of the reduction factor R, based on the reduced 
modulus concept in most cases leads to better venales than those based ce 
_ the tangent modulus approach. However, it is noteworthy that for Strut 5 better _ 
results for the first buckling load are obtained by using R, based on the tangent 
ws _ modulus procedure. This result can be anticipated, since Strut 5 initially was 
¥. caused to yield in tension, and it was ‘Straight prior to the application of a 
my, 
_ Considering the complexity of the p1 problem, the ; predictions for buckling loads — 
using the aforementioned approach may be said to be satisfactory and may 


prove useful in developing algorithms for determining the deterioration of the 


ONCLUSIONS AND Recommenoarions 4 ‘ane, Buckling 

} - Based on this investigation of inelastic cyclic buckling of struts of practical _ 
: sizes of various cross-sectional shapes, several conclusions may be reached = 
which have important design implications. yilvlaimg sdb 


. The conventional definition of an effective slenderness ratio K4/r deduced 


- on an elastic basis carries over into the inelastic range. The points of inflection 4 
on a deflected curve remain relatively fixed. 


_ 2. The effective slenderness ratio of a member ‘ appears to be the sles most : 


{191 sterensy 
| | | any 15) 
i 
q 


baste. 


"generate fuller loops than the more slender ones. The use of normalized hysteretic 7 
curves in comparisons is particularly advantageous because a number of variables 
_ are removed from consideration by this process. == oe | 
3. Hysteretic envelopes provide a convenient means for comparing specimens 
‘ ‘with different loading histories. They can be very useful in considering random 
loading effects on a brace during a severe earthquake. The use of normalized 
oh hysteretic envelopes is convenient for general formulations and studies. = 
_ 4, The hysteretic performance of a member is somewhat influenced by its | 
- eross-sectional shape. The major determining factor appear to be related to 
a member’ susceptibility toward lateral- torsional local of 


 -§. Stitching of built- -up critical c compression members for s service under severe : 
load reversals as currently specified in standard codes (7) is unconservative. © 

In the regions of plastic hinges, the individual parts of a member, due to softening . 
: of the material, have a greater tendency to buckle than envisioned by the codes. 

) Requiring the slenderness ratio //r for the individual parts of members between 
= ‘stitches to be less than | that of the ‘member as a whole and specifying minimum 
_ fastener strengths would help in the problem. Just this kind of a provision 
was contained in the 1959 AISC Specifications (12) and can also be found " 

in the current German ones (4). However, it would appear that for important 

a applications in seismic design, where severe cyclic loading of a compression 

member can be anticipated, in the regions s of potential plastic hinges built-up 
a members should be either avoided or very thoroughly stitched together, == 
_ 6. Significant reduction in buckling loads occurs during inelastic cyclic ieadion. 
The hysteretic loops displayed in this paper can serve as an aid for developing 
and verifying computer models of strut behavior. The use of the reduction 
factors 1 mentioned i in | the text may prove useful in such One 


; with the Bauschinger effect that occurs in 1 the plastic hinge regions ofa member. 

T he second factor accounts for the residual curvature that remains in a member 
following previous inelastic compressive cycles. 
7. Designers should be aware of the variability in the mechanical properties 
of commercially available steel used in building construction. The yield — 
ee rolled sections used in these experiments made of A36 steel varied over 

4s relatively narrow range from about 36 ksi (250 MPa) to 50 ksi (340 MPa). > 
_ Whereas and tubes made of Grade B steel and onal, 


MPa). 
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DISCUSSIONS 
| 
1 
us 4 4 we a A u a ar 
; 7 of manuscript including figures and tables); the editors will delete matter extraneous to the a _ 
_ subject under discussion. If a Discussion is over two pages long it will be returned for shortening. : 
Discussions will be reviewed by the editors and the Division’s or Council’s Publications 
> 


STABILITY OF FLEXIBLY Connecten PLATE Systems" 


fol following corrections should | be made to the original Paper: 
2077, paper title: Should read ‘‘Flexibly’’ instead of ‘‘Flexibility” 


2077, paragraph 1, line 4: Should tead “factor | is a function’ instead 


of “‘factor = a function” yet Sietrd, tien 
Page 2077, paragraph 2, line 6: Should read “welding and or 
pa 2077, paragraph 3, line 3 Should read “rigid” instead of “ridid”” gt ten 
“‘generally”” 

2084, paragraph 2, line 1: Should ‘read ‘ . For even only k,”’ 
Page 2085, 3, line 5: Should read “When”’ instead of ‘‘when’”’ 
2085, paragraph 4, line 1: Should read “local” instead of ‘“‘total’ 
4 Page 2087, Table 1 Note: Should read “ke are for a continuous (100%) 6.3 = 


2096, line 12: Should read ‘ instead of “a, = 8B, owe hag 


Page 2096, line 13: Should read “a, b,: 8,9." instead of “a, 6, = 
RESEARCH NEEDS IN STABALATY or Meta Structures” 


ian The writers wi wish to to congratulate Bjorhovde for the the nll State- of-the- art 
re report on the research needed in the field of stability 0 of metal structures. As 
the paper reflects, mainly, the viewpoint of the Structural Stability Research 
- Council (SSRC), the writers believe that their respective positions in Technical 
Committee TC8 of the European Convention for Constructional Steelwork (ECCS) 


a “October, 1980, by Milija Pavlovi¢, Donald S. Mansell, and Leonard K. stevens (Proc. F 
December, 1980, by Reidar Bjorhovde (Proc. Paper 15926). 
Asst. Prof., Univ. of Iniv. of Lidge, Belgium. 
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and the fact that the senior writer has for many years, been a corresponding — 
member of the SSRC, allow them to comment, somewhat, and eventually to 
the valuable information brought by their overseas 
isdone asrepresentatives of 
i Centrally Loaded Columns.—Indeed, the scientific study of the occurrence 
of residual stresses during the manufacturing process is needed in order to 
better know the connection between the structural properties and imperfections 
_ and the fabrication conditions. In this field, research on the effects of —— 
‘ on bucklingis under way in Liege, 
_ Simple methods of residual stress prediction should be welcome, 1 more 
especially, for welded built-up sections. A step forward has been made recently — 
in this field (42). The asymmetric shape of the cross section is not clearly 
known and a limited number of tests on asymmetric H shapes will be performed 
in Liege inthe nearfuture. = 
_ The problem of an out-of-straightness larger t than L/1,000 is encountered. 
sa that the residual stresses are the same as in regular shapes, a simple ~ 


method has been suggested (43) to evaluate the ultimate load of such columns. © 

wa is of the opinion that, as the new European buckling « curves are e more 
same way, one must be c cautious ‘whea defining ‘the effective. which, 

despite its criticism, remains a useful design tool. More importantly, the effective 
length of compressed chords of trusses for buckling in the plane of the truss, 
should be raised from 0.8 / to 0.9 /, as clearly shown by numerical simulation, 
taking into account elastoplastic behavior and large displacements. This has 
already been done in some national standards (SIA 161: Swiss “Standard on 
Steel Construction, 1979, and wa BS51- a Belgian Standard on Steel Construc- 
1977). 


. of a column belonging to a reticulated structure—the column being connected 
to other beams and columns at its joints—could be somewhat reduced compared © 
to the behavior of a centrally loaded column. The writers also believe that 
many structural columns should be designed as beam-columns, and that valuable — 
simplifications could be achieved by the definition of fewer column curves 
in these cases. Unfortunately, such thinking is not yet sufficiently supported 
by tests to be completely and accurately reflected in design rules. _ th 
Biaxially Loaded Columns.- -—Design rules for biaxially loaded that 
- are extrapolations of rules for uniaxially loaded columns have been tested recently 
_ (44). It appears from the results that both the Belgo-American and ECCS 
interaction formulas behave satisfactorily. The Belgo-American approach formula 
is approximately 3% more conservative than the ECCS formula; the Chen-Atsuta 
‘approach is generally on the safe side, but with a greater proportion of unsafe 
results than the aforementioned formula. It is shown, in addition, that most 
current equivalent moment formulas lead to a very marginal difference, a 
a slight advantage tothe AISC/ECCS formula, se 
‘The simulation tool first available for buckling of centrally loaded columns _ 
e no longer valid for biaxially loaded columns, due to torsional effects. Therefore, — a 


a new type of finite. element is needed to achieve this “goal. Such a a work is 


_ Tubular Columns. —Although interesting progress has been made in this field, 


| 
| 


this topic will not be reviewed here because a a separate discussion devoted — 
to ints presented by J. Rondel. 
Laterally Unsupported Beams.—Indeed, many tests are available on the lateral 
buckling of beams, and it has been clearly shown (45) that the behaviors of _ 
hot-rolled sections and welded built- “up ‘Sections are significantly different. It 
is obvious that residual stresses, initial out- -of-straightness, and cross- -sectional a 
_ distortions have a substantial influence on the ultimate strength by lateral buckling. * 
However, in the writers’ opinion, lateral buckling depends more on the more. 
or + less precise manner in which the actual conditions of support and load 


effects. while the support conditions remain rather too ideal. seat, alee 


__ The author states that the danger of lateral buckling is greater during erection 


of the structure, when it is not yet braced, than in service conditions, which _ 
: benefit from stabilizing secondary elements. There is, in addition, a range of 


problems that Tequires a attention | regarding the level of applied loads—problems _ 


at attachments located on the compressed bore) 
3 Stability of Frames.—The writers wish to emphasize the fact that, in spite 
7 of valuable information drawn from past research, the Gesigner needs a tool 
will | eliminate the need for sophisticated calculations, taking into 


frames, with respect to the requirement of a second order theory, and to provide — _ 
7 simplified second order methods of design. In Europe, the trend is to simplify 
on the basis of two classes of values for the sway angle } at each floor and 
two classes of values of coefficient e= aV N/N., , resulting in four approaches 
of increasing complexity, from the first order theory with no P-A and no oc 
effects to the full second-order theory with full account of both P-A and e« 
effects. Such a proposal, already partially introduced in German specifications ‘ 
(50), should be published in the first edition of Eurocode3. 
_ A practical method of second order design for sway frames has been suggested 7 
by Vandepitte.. It results from a modification of the horizontal equilibrium 
equations used in the regular slope deflection method. Such a ‘procedure leads 
to a noniterative system of linear equations and can be generalized for the _ 
_ Plate and Box Girders. —Plate and box girders are two fields in which significant — 
progress has been made during the past decade. In determining the ultimate 


Strength of transversally stiffened plate girders, ‘mathematical, models are now 


cr? 


available which provide accuracy. Unfortunately, when the web is stiffened, 
the design of longitudinal stiffeners is not sufficiently solved to allow the use 
of these ultimate models for designing longitudinally stiffened webs. The writers 
feel that stiffeners designed according to the linear theory of buckling (by 
‘their relative rigidity by a factor of the order of 4 to 


that | they remain rigid up to collapse) provide a a check for the buckling of oe 
plate subpanels and that, in this case, the postcritical strength reserve can be 
= appraised on the basis of a were oes numerical investigation performed 


Technical TC buckling) of ECCS is ) preparing | 


j 
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of in the area of stability problems with 
construction. Three other books in initial stages of development deal respectively : 

1. Stability behavior of components. 


Regarding fabrication tolerances and their effect on . ultimate strength, a a S. 0.A. 


; ‘mentioning (Tolerances in Steel Plated Structures. IABSE Surveys, S-14/80.) 4 
Composite Columns.—A model code has been prepared by an international | 
joint Committee ([ABSE-FIP-CEB-ECCS) (47). It deals with the design and hg 
- construction of composite structures and structural parts, and was designed © 
= use as a common basis or reference for national and international codes 


or specifications. The design method for composite columns presented in this 
model code is based on an ultimate load design philosophy. It is very broad 
i scope and covers encased columns and concrete filled circular and rectangular 
“steel tubes, which are loaded with any combination of end loads and moments. 
_ The code’s basic assumptions are that there is full interaction between concrete — ? 
-. steel up to the point of collapse, that proper account must be taken of 
the steel and concrete stress-strain curves as given in the ECCS and CEB 
~Recommendations dealing with steel and concrete respectively, and that allowance z 
must be made for imperfections whi ec 

assessing the strength of axially loaded base steel ‘columns. On the. basis of 
_ last assumption a new concept of column slenderness is introduced which 

leads to the same expression as that used in the European buckling curves 

and enables | these « curves to be used as the bas basic design curves for composite 


‘results on - several types of composite columns and agreement is shown to be 4 
_ Shells.—The influence of the amplitude and distribution of imperfections on 
; the ultimate strength capacity of shells is well known. Studies dealing with 
_ the connection between imperfections and ultimate loads would be welcome. 
Longitudinally stiffened cylinders axially or eccentrically compressed is an area . 
of practical interest in this field. The least investigated field concerns cylinders — 
42. Plumier, A., “Relation between the Welding Parameters and the Ultimate Buckling : 
Load of Columns,”’ thesis presented to the University of Liége, Belgium, in 1979, 
in partial fulfillment of the requirements for the degree of Doctor of Philosophy. gaily 
-- Rondal, J., and Maquoi, R., ‘‘The Buckling of Steel Columns,’ Chambre Syndicale 
des Fabricants de Tubes d’ Acier, Paris, Notice 1091,1980. 
44. Anslijn, R., and Massonnet, Ch., “New Tests on Steel I Beam-columns in Mild 


bg “7 Steel Subjected to Thrust and Biaxial Bending,” Festschrift Otto Jungbluth—60 Jahre, 
Keller, ed., Technische Hochschule, Darmstadt, 1978. 
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45. Fukumoto, ae Kubo, M., “An Experimental Review of Lateral Buckling of s 
- Beams and Girders,”’ Presented at the May 17-19, 1977, International Colloquium 
on Stability Structures under Static and Dynamic Loads, Washington, D.C. 0 cory 
6. Harding, J. E., and Dowling, P. J., ‘“‘The Basis of the Proposed New Design Rules 
“a for the Strength of Web Plates and Other Panels Subject to Complex Edge Loading,” 
ei) Stability Problems in Engineering Structures and pane Compenenny, | Appl. | Science P Publish. } 
Ltd., Barking, England, 1979, pp. 355-376. 
47. “‘Draft Model Code for Composite Structures,” Recommendations of the CEB- ECCS. 
_ FIP-IABSE Joint Committee on Composite Structures, The Construction Press, London — 
- 48. “The Design of Composite Steel Concrete Structures,’’ Second International Coll 
quium on Stability,” Introductory Report, Liege, Belgium, April, 1977. 
5 49. “*Mixed Construction,” Structural Design of Steel Buildings, Volume SB, Chapter 


50. Vogel, U., and Lindner, J. “Commentary on DIN 18800, Teil 2 (Gelbdruck)—Stabilitats- 
= _ falle im Stahlbau,”’ Knicken von Staben und Stabwerken, in ——~ DASt Berichte © 
aus Forchung und Entwicklung, Nov., 1981, hoe 


of plow the wishes to call attention to considerable 

in this field , the results of _ which were published after the author’s most — 

_ comprehensive survey of research needs was first presented in April, 1979. 
This work was done under a two-year research contract sponsored by the 


_ engineering knowledge i in this field, in the United States and Europe, are presented 
in the final report 0 on the ‘ ‘Proposed Design Specifications for Steel Box Girder ; 
The report contains a detailed discussion and proposed tentative engineering 
“solutions to the questions listed by the author (box girder webs, diaphragms, © 
_ Shear lag, fabrication tolerances) and other problems, such as the strength of 
7 compression flanges, interaction between the webs and the flanges in the cases 
_ where webs are designed with utilization of the ultimate (tension field) strength, 
design of the transverse, and the horizontal web stiffeners 


a a _ Criteria proposed for evaluation of the strength of stiffened compression ny 


q based on the “column approach, with consideration for the effects of 


Out- |-of-straightness and residual stresses. 
Of all components of a box girder the design ow the webs presents the most 
- difficult problems. Design recommendations given (51) are based on maximum — 
utilization of the elastic beam strength of the web subpanels (requiring rigid 
-_siffeners) and partial utilization of the tension field strength (corresponding ~ 
to “true Basler” solution) . The solutions | developed for the box girder) webs _ 
a also applicable, with some necessary modifications, to Oe HE of plate 
_ The proposed rules are necessarily conservative. Their liberalization and 


Wolchuk, | 4.0 E York, 
“Partner, olcl uk, , and | Mayrbeurl, Cor ons. Engrs., New Yor 
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and more fully answered by continuing research. 
Suggestions for future research on box girders are given in report: (42) pp. 
5-16 16 where they hay have be been su: summarized as follows: 
a) Parametric studies of applications of the proposed design rules for 
webs, web stiffeners, unstiffened and stiffened flanges in compression 
are very desirable. Such studies will permit evaluation of the range of 
gs and stiffener sizes, determine which of the various parameters oe 
principal importance in practical design, and may indicate 
improve and simplify the proposed rules. ne 
It is generally that many questions pertaining to the design 


tively solved and remain to be anowered by continuing research. Further E 
theoreti and experimental work on these problems i is much needed ded. 
contains an extensive bibliography of the subject. 


in (53). 


“Proposed Design Specifications ter ‘Steet Girder Bridges,” 
TS-80-205, Wolchuk and Mayrbaurl, Consuiting Engineers, Federal Highway Admin- 
® Wolchuk, R., ., “‘Proposed Specifications for Steel Box Girder Bridges,’’ Journal of 


the Structural Division, ASCE, Vol. 106, No. ST12, Proc. Paper 15946, Dec., 1980, 


33. Wolchuk, R., “Design Rules for Steel Box Girder Bridges,”” IABSE Proceedings, 


‘mei author investigates ss the ¢ capacity of shells shells to resist membrane forces and 
= procedures for the design of the associated membrane reinforcement. 
However, as this discussion will indicate, the writer believes that the capacity 
determination and design procedures of the paper are not in accord with the ‘ 
_ strength criteria and the strength design method of the ACI 318-77 Code (2). 
_ The strength design method is summarized in section 9.1.1 of ACI 318- 7 


“January, 1981, by Ajaya K. Gupta (Proc. Paper 15975). 
“Cw. Engrg., Specialist, Burns and Roe, Woodbury, N.Y. 
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d of the proposed solutions is discussed concisely + ; 
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as structural me: members shall be designed to have re design 
_ Strengths at all sections at roa equal to the required strengths calculated for | 7 
factored loads and forces. . . .”” The strength design method is based on — 
satisfaction of equilibrium and strength (yield) requirements and thereby complies _ 
_ with the lower bound theorem of Limit Analysis (16). It then follows from - 
Limit Analysis that the structure has sufficient capacity to sustain the applied — 
factored loads. Determination of the critical failure mode or cracking pattern — 
=u which is taken up in the paper is not involved in the strength design method 


_orinthe lower boundtheorem. = 
7 2 The ‘Strength ci criteria | of ‘ngs 318- “17 relative to the combined action of tension 


the reinforcement required for ‘each ‘type of stress be. explicitly. 7 
‘The shear strength provided by shear reinforcement is stated in section 11.5.6. 2 a 
when shear strength is governed by diagonal tension and in section 11.7.3 when 
ow Using the terminology of ACI 318-77 and notation in line with the paper, % 
N re 2, 83.8 * an and N * are the design strengths of the shell relative to membrane 
forces and N,, N,, and N ,, (equal to N,,) are the required strengths or internal 
forces which are in equilibrium with the applied factored loads. By strength 
design, the design equations are st) 
The design equations s for the are obtained by 
using the provisions of the Code for tensile strength and shear strength. After — 
- Setting N *, = N¥, and summing the reinforcement requirements in each direction 
for tension and. ‘shear, , the following sets of equations are obtained for the — 
- total reinforcement in each direction. Case 1: Shear strength controlled by shes shear 


if Case 2: ‘Shear strength controlled by diagonal tension. 


= oV, unit width) — wing 
‘The application of the design method as described above can be 
ac accomplished in any ny selected coordinate system with respect to the equilibrium. 
~ analysis and the provision of design strength by orthogonal reinforcement. a 
_ Investigation of the adequacy of the design relative to other coordinate systems 
_ which, the author indicates, is a necessary part of the design process is not 
required. This follows from the lower bound limit theorem which, in alternate — 
form (17), states that if an equilibrium dist distribution of stress can ote cients 


(ie., @,, the direction of the yield line) and the magnitude of the collapse 
& for a given concrete slab. Its application requires that the design strength 


be applied, as in references (6) and (18), to determine the critical failure mode 


_(N%,N*}) and the internal forces due to a loading pattern (N,,N,,N,,) be 
_ known. The internal forces due to the collapse load are represented in terms 
_ of a load parameter k by kN, kN,, andkN.,. The parameters to be evaluated | 

_ by the principle are © and k. Accordingly, it is apparent that the principle — 
of minimum resistance has no direct application to the design problem of 

- determining the amount of orthogonal reinforcement required to sustain a given _ 


7 .. factored load and that the principle is not pertinent to the Strength Design 
sociated yield criterion 6 are based on the principle of minimum 
resistance. However, Eqs. 5a and 5b (of the original article) are obtained directly - 
_ from Egs. 2, 3, and 4 and do not utilize the minimum principle Eq. | at all. 
. fact, it can be readily shown that the critical value of © obtained by using 
the principle via Eqs. .& and 4 bears no relationship to the values of © 
__¢. In view of the preceding it appears that Eqs. 5a and 5b have been adopted _ 
_as “‘the design equations for any assumed value of ©”’ on the basis of satisfying i 
- Egs. 2, 3, and 4. In this formulation the strength criterion is represented by ; 
Eq. 4. Review of this equation shows that, in contrast with aforementioned 7) 
ACI Code provisions, no explicit strength contribution due to the shear strength 


of reinforced concrete is included in the strength criterion. The resulting deficiency 

in the paper’s strength criterion is apparent when the design Eqs. 5a and 5b 
are applied for values of © = 90° and © = 0°. The design equations would 

then require infinite N * and infinite N * respectively and cannot be satisfied. _ 
 d. The statement by the author, leading to Eq. 7, that a compressive force 3 
parallel to the crack is necessary in order to maintain internal equilibrium, is’ 

- not concurred with. The internal forces acting on the shell membrane elements 
which have been determined by prior analysis to be in equilibrium with each 
other and with the applied factored loads are generically the forces N,, N,, 
and N,,. The actual internal state of stress relative to any direction follows 


from the values of N,, N,, and N,, and can be pictured via the associated 
Mohr circle of stress. The case of N, and N, both positive, has been pictured ; 

7. in Fig. 11. For the relative magnitudes used, it is apparent that no compression _ : 
: in any direction exists. The author considers the case in which the orthogonal — 


4 

_balances the applied load and is everywhere less than or equal to the design + ; 
- strength (yield stress), the structure will not fail. Physically, it follows because — 
the membrane element cannot fail if the design strength relative to two aa“ 
- _ Comments on the development of procedures in the paper are furnished below ; 
_ taking into account the preceding discussion of Code methods and criteria. 
a. The principle of minimum resistance, which is utilized by the author, can_ 
| 
if 
| 
i 
q 
> | 
ii 
. wrcinforcement is at the design strength levels, N* and N*. However, the ; 
sociated stress state is, by definition, N,, N,, and N,,, so that the absence &§ 
a e. The state of strain is investigated by the author assuming continuity of 
| 


In this segerd , the additional Eq. 50 below must be ntinted, besides: Eqs. 
z and 13, to. define the relationship between eve, 


Also, for case of N and N, at capacity strees, ws, paper incorrectly 
_ postulates that «, must be compressive. The invalidity of this assumption has 
been discussed i in preceding paragraph d. In addition, inconsistency i in the meaning — a 
of the angle © occurs in the “paper. By definition, via Eqs. 12, 13, and | 


© is the angle between the direction of the largest med strain and the 


ten , 
wa. 


SHEAR STRAIN 


NORMAL 


12.—Mohr Circle of Strain for Maximum «, 


sions relative to the iproperly omits any consideration 
t 
! 
= 


udivection. However _ it is stated that the minimum size cracks (or 
Strain) occur when the reinforcement and the principal strain are in the ‘same_ 
‘direction, i.e., with = 0°, Fig. 2 shows minimum strain when 45°. 

_ The unrealistic results which are obtained with the author’s methods are 
a apparent from Fig. 2 which pictures Eqs. 14 and 15. For the case of 6 = 
_ the orthogonal reinforcement lies in the principal directions so that €, — 
= c, and €, = «, as is evident from the standard strain Eqs. 12 and 13. This - 
contrasts 3 with the value of ¢€, from Fig. 2 and Eq. 14 which indicate infinite — 

= The author ‘reaches the conclusion that rapidly increasing strains result as 


the angle 9 between the reinforcement and the principal strain direction increases 
from 0° to 45°. The validity of this conclusion is now considered using the 
“Mohr circle of strain as in Fig. 12. Towards this end, we seek the largest 
ratio, R=€, /«, with «, where ¢€, is the tensile strain at yield. For any 
values of € 4 1g ‘ this will occur with Yep at maximum value. rd 


and are such as to preclude rete strains beyond the elastic range. On this - 
basis, a maximum shear strain y,, , of about 0.8¢, can = derived. The largest 

ratio R = /e, then occurs for €, = €, and y,, 0. and this case 

is pictured in Fig. 12. As is onions | from the figure and Eq. 50, the largest 
ratio R=1 4 and © = 45°. With other values of the orthogonal strains €,, €,, 
the ratio of the principal tensile strain €, to the tensile strain in the reinforcement 

€, will be between I. 0 and |. 


16. Wood, R. H., Plastic and Elastic ot Slabs and Plates, Ronald 
or 7. Plastic Design in Steel, Manual No. 41, 2nd od. by a Joint Committee of the yee 
Research Council and the American Society of Civil Engineers, ASCE, 1971, p. 10. — 
18. Jain, S. C., and Kennedy, J. B., ‘Yield Criterion for Reinforced Concrete Slabs,” 
Journal of the Structural Division, ASCE, Vol. 100, No. vs, Mar., 1974, pp. 631- as. 
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must be legible at the reduced size. Photographs should be submitted as glossy prints. Explanations a 


must be — in text rather than within the — 


o Silineis: cited in text should be arranged in alphabetical order in an appendix at oo ; 
of the paper, or preceding the Appendix.—Notation, as an Appendix.—References. 


10. A list of key words and an information retrieval abstract of 175 words should be provided 


= 


of 40 words must accompany the 


12. A set of 
set of conclusions must end the paper. 


13. Dual units, i.e., U.S. Customary followed by SI (International Sued units in pa entheses, 


14. A practical applications section should be included also, if appropriate. 
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